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Baryon number violation is discussed in gauge unified orbifold models of type II string theory 
with intersecting Dirichlet branes. We consider setups of Z)6-branes which extend along the flat 
Minkowski space-time directions and wrap around 3-cycles of the internal 6-d manifold. Our study 
is motivated by the enhancement effect of low energy amplitudes anticipated for M-theory and 
type // string theory models with matter modes localized at points of the internal manifold. The 
conformal field theory formalism is used to evaluate the open string amplitudes at tree level. We 
study the single baryon number violating processes of dimension 6 and 5, involving four quarks and 
leptons and in supersymmetry models, two pairs of matter fermions and superpartner sfermions. 
The higher order processes associated with the baryon number violating operators of dimension 7 
and 9 are also examined, but in a qualitative way. We discuss the low energy representation of string 
theory amplitudes in terms of infinite series of poles associated to exchange of string Regge resonance 
and compactification modes. The comparison of string amplitudes with the equivalent field theory 
amplitudes is first studied in the large compactification radius limit. Proceeding next to the finite 
compactification radius case, we present a numerical study of the ratio of string to field theory 
amplitudes based on semi-realistic gauge unified non-supersymmetric and supersymmetric models 
employing the Z3 and Z2 x Z2 orbifolds. We find a moderate enhancement of string amplitudes 
which becomes manifest in the regime where the gauge symmetry breaking mass parameter exceeds 
the compactification mass parameter, corresponding to a gauge unification in a seven dimensional 
space-time. 

PACS numbers: 12.10.Dm,11.25.Mj 



I. INTRODUCTION 

By suggesting the possibility that the Standard Model (SM) gauge and gravitational interactions unify in a higher 
dimensional space-time, string theory [1, 2] has opened up novel perspectives for the gauge unification proposal [3]. 
Based on the approach of Calabi-Yau (CY) manifold compactification, the string theory applications focused initially 
on the heterotic string in the regimes of weak coupling [4] and strong coupling [5-8] . These studies were soon followed 
by explicit constructions using orbifold compactification [9, 10] and free fermion [11] models of the heterotic string, 
which were later extended to orbifold compactification models of the type // string theories with single _Dp-branes [12] 
or the Dp/Dp' , [p' ^ p] branes-inside-branes type backgrounds [13]. 

The chief distinctive features of string gauge unification reside in the discrete gauge symmetry breaking scheme by 

Wilson flux lines [14] and the heavy string threshold effects [15]. Regarding, however, the issue of matter instability 

caused by baryon number violation, it is fair to say that no specific stringy effect has emerged from the earlier studies 

using the weakly coupled heterotic and type // string theories. A different situation seems to take place in the 11-d M- 

theory supersymmetric compactification on 7-d internal manifolds Xi of G2 holonomy [16], as discussed by Friedman 

and Witten [17] . The non-Abelian gauge symmetries in these models are supported on 7-d sub-manifolds, B = M4 x Q, 

loci for R^/r orbifold type singularities in the directions transverse to B, while the chiral massless fermions are 

supported at isolated singularities of the 3-d sub-manifolds, Q. This causes the existence of a natural regime where 

the grand unification occurs in a 7-d dimensional space-time with localized matter fermions. In comparison to the 

nucleon decay amplitudes of the equivalent unified gauge field theories, the string amplitudes are enhanced by a power 
_ 1 

a-^^ of the unified gauge coupling constant, which reflects on the short distance singularities from summing over the 
momentum modes propagating in the sub-manifold Q. As to the size of the enhancement effect, however, no definite 
conclusions could be drawn because of the poor understanding of M-theory perturbation theory, not mentioning the 
greater complexity of G2 holonomy manifolds [18]. 

Fortunately, it is possible to examine the enhancement effect in a controlled way by considering the weak coupling 
dual models based on type I la string theory orientifold compactification on M4 x Xq with D6-branes wrapped around 
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intersecting three-eyclcs of the internal Calabi-Yau complex threefold, Xg. Using a toy model realizing SU{5) gauge 
unification, Klebanov and Witten [19] showed that the four fermion string amplitude for localized matter modes in 
the gauge group representations, 10- 10^ • 10- 10^, featured a power dependence on the unified gauge coupling constant 

2/3 

of same form as in the M-theory model, namely, A cx Oj^'' . In the most favorable case where the fermion modes sit 
at coincident intersection points of the D-branes, the enhancement due to the gauge coupling constant dependence 
was found to be offset by certain constant factors which resulted in a rather modest net effect. 

Our goal in the present work is to document the enhancement of string amplitudes anticipated in M-theory [17] 
by developing further semiquantitative calculations in models with intersecting D-branes. The initiating discussion 
in Ref. [19] made use of the large compactification radius limit in which predictions are largely insensitive to the 
Wilson line mechanism responsible for the unified gauge symmetry breaking. We extend this study in three main 
directions. First, we consider semi-realistic orbifold-orientifold models realizing the 5'L''(5) gauge unification with a 
chiral spectrum of massless matter modes that closely reproduces the Standard Model spectrum. Second, we evaluate 
the nucleon decay four point amplitudes in the two independent configurations of gauge group matter representations, 
10 • 10^ • 10 • 10^ and 10 • 5 • 10^ • 5^. which contribute to the proton two-body decay channels with emission of left 
and right helicity positrons, p ^ tt^ + ^. This allows us to quantitatively assess the M-theory prediction for the 

ordering of partial decay rates, r{p ^ tt^ + ej) >> T{p — > 7r° -|- e^). Third, we discuss the string amplitudes at finite 
compactification radii, in order to weigh the importance of the string momentum and winding excitations relative 
to the string Regge resonances. Studying the interplay between the compactification and unified symmetry breaking 
mass scale parameters, and Mx, proves crucial in assessing the size of the enhancement effect. 

The contributions to nucleon decay processes from physics at high energy scales are conveniently represented by 
non-renormalizable local operators in the quark and lepton fields [20, 21] which violate the baryon and lepton numbers, 
B, L. In unified gauge theories, the exchange of massive gauge bosons with leptoquark quantum numbers induces in 
the effective action, [B + L) violating, {B — L) conserving operators of dimension T> = 6. In supersymmetry models, 
the exchange of massive color triplet matter higgsino like modes also induces dangerous operators in the quark and 
lepton superfields of dimension V = 5. Higher dimension operators initiating the exotic nucleon decay processes can 
possibly arise from mass scales significantly lower than the gauge interactions unification scale, Mx ~ 3. x 10^^ GeV. 
Of special interest are the B — L violating operators of dimension V > 7, and the double baryon number violating 
operators of dimension T> >9, responsible for nucleon- antinucleon oscillation. 

The 2-d conformal quantum field theory [22] provides a powerful approach for calculating the on-shell string theory 
amplitudes. The construction of string amplitudes [23] from vacuum correlators of open string vertex operators 
inserted on the world sheet disk boundary is readily applied to the 4-d amplitudes describing the nucleon decay 
processes. The tree level contributions to the dimension D = 6, 5 operators involve four fermions, ^/j*, or two pairs of 
fermions and bosons, ip'^cjp. For the intersecting brane models with matter modes localized in the internal directions, 
the calculations are made non-trivial by the need to account for the twisted like boundary conditions of the world sheet 
fields. Fortunately, the energy source approach, which was initially invented [24-28] and subsequently developed [29- 
31] in the context of closed string twisted sectors of orbifolds, can be readily extended to the localized open string 
sectors of intersecting brane models. Following previous studies devoted to the discussion of open string modes with 
mixed ND boundary conditions in Dp/Dp'-hr&ne models [32-35], the implementation of this approach for intersecting 
brane models has been clarified in discussions by Cvetic and Papadimitriou [36] , Abel and Owen [37] , Klebanov and 
Witten [19], Jones and Tye [38], Liist et al., [39] and Antoniadis and Tuckmantel [40]. A comprehensive review of 
the first quantization and conformal field theory formalisms for the open string string sectors of D-brane models is 
currently under preparation [41]. We should also mention here the studies by Billo et al., [42], Bertolini et al., [43], 
and Russo and Sciuto, [44] which develop the alternative approach based on the operator formalism. 

An intense activity has been devoted in recent years to the discussion type // string theory compactification using 
intersecting Z?-brane backgrounds [45-47]. To be complete, we should mention the parallel development on T-dual 
£)-brane models using magnetised backgrounds [48, 49] and on the magnetic field deformations of the heterotic 
string [50]. Important efforts towards building satisfactory models have been spent in works by Aldazabal et al., [51- 
56], Kokorelis [57-59], Blumenhagen et al., [60 67], and Cvetic et al., [68-74]. A useful review is presented in Ref. [75]. 
To develop our applications in the present work, we consider semi-realistic models already discussed in the literature. 
First, we should note that the minimal toroidal orientifold models developed by Cremades et al., [53-55] are of little 
use to us in the present work, because these realize direct compactifications to the Standard Model (SM), the minimal 
supersymmetric Standard Model (MSSM) [53-55] , or the Pati-Salam type unified model [58] , which all accommodate 
a built-in U{1)b+l global symmetry which guarantees B and L number conservation. A better answer to our needs is 
provided by the orbifold models [60, 61, 66-72, 76-85] owing to their richer structure and more constrained character. 
We have selected the two classes of non-supersymmetric and supersymmetric gauge unified models realizing a minimal 
type SU{5) gauge unification. The first class relates to the works of Blumenhagen et al., [66] and Ellis et al., [78] using 
the Zs orbifold, and the second class to the work of Cvetic et al., [68] using the Z2 x Z2 orbifold. It is important to 
emphasize at this point that the calculation of string amplitudes relies on features of the non-chiral mass spectrum and 
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the wave functions of low- lying modes which are not directly addressed in these works. For instance, the application to 
Z3 orbifold models rely on data involving the real wrapping numbers rather than the effective ones which are defined 
by summing over the orbifold group orbits. 

The discussion in the present paper is organized into five sections. In Section II, we review the first quantization and 
conformal field theory formalisms for type // orientifold models with intersecting D6-branes. The open string sectors 
are discussed first for tori and next for orbifolds. A review of the energy source approach for calculating the c;orrclators 
of coordinate twist field operators is provided in Appendix A. In Section III, we discuss the calculation of string 
amplitudes for nucleon decay processes in the gauge unified models. The world sheet disk amplitudes are considered 
in succession for the processes involving four massless fermions, two pairs of massless fermions and bosons, four massless 
fermions and a scalar, and six massless fermions. These are described by operators of dimension, T> = 6 and 5, 2? = 7, 
and D = 9, obeying the selection rules, AB = AL = -1, AB = -AL = -1 and AS = -2, AL = 0. In Section IV, we 
discuss the relations linking the low energy gauge and gravitational interaction coupling constants to the fundamental 
string theory parameters gs and mg , and to the mass parameters representing the average compactification radius of 
wrapped three-cycles, r = 1/Mc, and the unified gauge symmetry breaking, Mx, arising as an infrared mass cutoff. 
We consider in turn two distinct regularization procedures of the string amplitudes. The first uses the subtraction 
prescription replacing the massless pole terms by massive ones, and the second uses the displacement prescription 
splitting the unified Z?-brane into separated stacks. To illustrate the dependence of four fermion string amplitudes on 
the iJ-branes intersection angles, an initial numerical application is considered within the large compactification radius 
limit. However, the main thrust of the present work bears on the study of baryon number violating string amplitudes 
at finite compactification radius. The results illustrating the enhancement of nucleon decay string amplitudes are 
presented in Section V. We consider first the non-supersymmetric SU{5) unified models of the Z3 orbifold due to 
Blumenhagen et al., [66] and Ellis et al., [78] and next the supersymmetric SU{5) unified model of the Z2 x Z2 orbifold 
due to Cvetic et al., [68]. In Section VI, we summarize our main conclusions. For completeness, we provide a brief 
review of the Z3 orbifold models in Appendix B and a brief review of baryon number violating processes for gauge 
unified theories in Appendix C. 



II. REVIEW OF TYPE // STRING ORIENTIFOLDS IN INTERSECTING i36-BRANE BACKGROUNDS 



We present in this section a brief review of the open string sector of type // string theory orientifolds with inter- 
secting _D6-branes. The first-quantized and vertex operators formalisms are discussed for toroidal compactification in 
Subsections II A and II B and for orbifold compactification in Subsection II C. 



A. Toroidal compactification 

1. World sheet field theory 

The Neveu-Schwarz-Ramond type Ila string theory deals with 2-d world sheet Riemann surfaces on which live the 
coordinate and Majorana-Weyl 2-d spinor fields parameterizing the 10-d target space-time with signature ( — h • • • +). 
For a fiat space-time, one uses the diagonal representation of the metric tensor, —7700 = Vti/j, = Vmm = 1; associated to 
the orthogonal field basis, , , [M = (/i, m) = 0, 1 • • • , 9, = 0, 1, 2, 3; m = 4, 5, 6, 7, 8, 9] with the complexified 
basis of coordinate fields, X^, X^, [A = 0, 1, • • • , 4] defined by 

X^ = [X^; X™] = [A:^; X^,X\ [A={a,0) = (1,4,0), /= 1,2,3] 



V2 ' ' ' ' ^/2 



V2 V2 



(11.1) 



with similar linear combinations for the complexified basis of spinor fields, tp^, We consider the orientifold 

toroidal compactification on M4 x r^/(f27?.(— 1)^^), with factorisable 6-d tori, = nj=i ^/ > symmetric under the 
product of parities associated with the world sheet orientation twist, O, the antiholomorphic space reflection, TZ, and 
the left sector space-time fermion number, F^. For dcfinitcness, we choose the orientifold refiection about the real 
axes of the three complex planes, TZ = TZy : X^ X^ . The Tf tori may also be parameterized by the pairs of hatted 
real lattice coordinates, X^ = {X(, X^), [a = 1, 2] of periodicities 27r each, X^ — X^ + 27r. Except in special instances 
where the dependence on a' is made explicit, we generally use units for the string length scale such that, a' = 1. 
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Wc focus our discussion on the open string perturbation theory tree level with the world sheet surface given by the 
disk. The complex plane and strip parameterizations of the disk are described by the variables, z € C_|_, z € C_, 
and a E [0, tt], t E [—00, +00], related as, z = e*"" = 6**^+*, z = e"*"" = e~*'^+*, with the derivatives replaced as, 
dcr ^ i{d — B), dt ^ —{d + B). The string equations of motion are solved by writing the coordinate and spinor fields 
in terms of independent holomorphic and antiholomorphic (left and right moving) functions 

X'^izrz) = X''{z)+X^iz), = (V'^(^), i^'^iz))^, (II.2) 

living in the upper and lower halves of the complex plane. The complexified spinor fields are also cquivalcntly 
represented by the complex boson fields, H^{z), H^{z),[A = (a,0) = (1,2,3,4,0)] through the exponential map, 

(z),V'-^(z)] ~ e±^-^''(^) and [Vi-4(z), ^'^(z)] ~ e^'"^^^\ 
The coordinate fields obey one of the two possible boundary conditions at the open string end points, a = (0,7r): 
Neumann or free type, d„X^ = 0, and Dirichlet or fixed type, dtX^ = 0. The Neumann and Dirichlet boundary 
conditions (abbreviated as N and D) for the coordinate and the spinor fields read in the complex plane formalism 
as, {d - B)X'^ = 0, (V-^ T V-^) = 0, and {d + B)X^ = 0, (V'^ ± V'^) = 0, with the upper and lower signs 
corresponding to the R and NS sectors. In the doubling trick representation of the open string sector, one joins together 
the holomorphic and antiholomorphic coordinate and spinor fields into single holomorphic fields, X^{z), tp^{z), 
living in the full complex plane C, by means of boundary conditions along the real axis, x = ^(z) = ^(z) G R. 
These are implemented by having the holomorphic fields coincide with the left moving fields X^{z), tp^ {z) on 
C+ and identifying them with the right moving fields, up to a sign change for D directions, on C_. Specifically, 
X^{z) = ±X^{z), ifj^ [z) = ±-(/'*^(z), [z G C_] where the upper and lower signs are in correspondence with N 
and D boundary conditions. The holomorphic superconformal generators are defined in the lower half complex plane 
as, T{z) = T{z), Tp{z) = Tf{z)- The boundary conditions on the coordinate and spinor fields exactly combine 
to leave invariant the linear combination of supersymmetry generators, Tp(z) + Tp{z), corresponding to the 2-d 
supersymmetry transformation, 5X = i{eip — e^), Sip = edX, 6ip = —edX, with e = — e. For completeness, we 
note that in the presence of constant fluxes of NSNS and magnetic two- form fields, Bmn and Fmn, the doubling 
prescription for N directions [86-88] with the constant metric tensor, Gmn, takes the matrix notation form, X^{z) = 
{D-^)mnX''{z), [D = -G-^ + 2{G + J^)-\ J^MN = BMN + 2TTa'FMN, zeC-]. 



2. Geometric data 

We consider general non-orthogonal (oblique) 2-d tori symmetric under the reflection TZy. For definiteness, we 
choose the upwards tilted Tf tori with lattice generated by the pairs of complex cycles or vielbein vectors, e{ = 

r(e'^^~" )/sina^, 62 = iri, where denote the opening angles and rf, the projections of the cycles radii along 
the real and imaginary axes of the complex planes. The diagonal complex structure and Kahler moduli parameters 
of the Tf tori are then expressed as 

U' = Ul + iUi = -4 = 6^ + i^, 

J.I 

T' = tI + iTi = + iVDet& = + ir[ri, [4 cot a' = -h'] (II.3) 

^2 

where the 0.1^ symmetry restricts the 2-d tori tilt parameters to the two discrete values, 6^ = 0, \. The hat symbol 
on V is meant to remind us that this parameter idcntifles with the NSNS fleld VEV or flux in the factorizcd T-dual 
picture [89]. The parameters V should not hopefully be confused with the continuous angle parameter, 6^ ^ 6^ + 1, 
in the Kahler moduli . The complexified orthogonal and lattice bases for the Tj coordinates, X/2 and XI2, are 
related by the formulas 




(11.4) 



with similar formulas holding for spinor fields. 
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The fixed sub-manifolds under TZ define the orientifold 06-hyperplanes which are sources for the dual pair of closed 
string sector RR form fields, C\, Ci. With the choice TZ ~ TZy, the 06-planes extend along the three flat space 
directions of M4 and wrap around the three one-cycles along the real axes of Tf. The need to neutralize the 
net RR charges present inside the internal compact manifold Xq is what motivates introducing iPe-branes similarly 
extended along M4 and wrapped around three-cycles of T^. Recall that the -Dp-branes arise as soliton solutions of the 
type // string equations of motion for the closed string massless NSNS (supergravity) modes associated to the metric 
tensor, Gmn, and dilaton, and for the massless RR dual antisymmetric form fields, Cp+i ~ Cr-p- The linkage to 
open strings is realized by the characteristic property of the -Dp-branes to serve as boundaries or topological defect 
sub-manifolds, immersed in the 10-d space-time, which support the open string end points. Since the RR charges enter 
as central terms in the supersymmetry algebra, the supersymmetric Dp-hianes {p = 0, 2, 4, 6 for Ila and p = 1,3,5 
for lib) preserving a fraction 1 /2 of the 32 supersymmetry charges in the bulk, satisfy a Bogomolnyi type bound on 
their mass which guarantees them stability against decay. 

The {p — 3)-cycles 11^ of wrapped by Op-planes and Dp-branes solve the string equations of motion as 
equivalence classes for closed sub-manifolds modulo boundaries, hence as elements of the homology vector space, 
[n^] e Hp-3{X, R). The dual relationship with the cohomology vector space, Hp~^{X, R), generated by the equiva- 
lence classes of closed differential forms modulo exact forms, is used to define the cycles wrapping numbers, or electric 
and magnetic RR charges, in terms of volume integrals. The cycles intersection numbers are defined in terms of 
topologically invariant integrals obeying the antisymmetry property, [11^] • [Hi,] = — [n^] • [11^] = The homology 

basis of fundamental one-cycles [a^] and [b'^] for Tf, have the intersection numbers, [d^] ■ [b'^] = —[u] ■ [a"'] — 5^'^ ■ To 
avoid filling the internal space, one commonly restricts to the integer quantized homology classes, Hp-3{X, Z). For 
convenience, one also usually limits consideration to the subset of factorizable three-cycles, 11^ = Y[i=i n^, products 

of one-cycles of the Tf tori. In the lattice bases generated by the dual bases of one-cycles, (a^, 6^), along the Tf tori 
vielbein vectors, (ef , 62), these cycles are parameterized by the three pairs of integer quantized wrapping numbers, 
{nj^.mjj. The orthogonal bases representations for the factorisable three-cycles and their orientifold mirror images, 

= 11/=! n^J,, n^' = n/=i n^^') denoted by (n^, m^), {nj^,,ifij^,), are related to the lattice bases representations as 
= nj,[a'] + m'^[b'] = nj,[a'] + mj,[b'], 

[nl,] ^ [a'] + ml, [b'] = ni[a'] - m^[6^] = <[a^' 

The invariant volume of three-cycles, Vn^ , and the topologically conserved number of intersection points of three- 
cycle pairs, 7^1/, are given by 

1=1 

3 

J^. = [n^] ■ [n.] = H {nj^mi - nimj,). {11.6) 
7=1 

The one-cycles in Tf can also be described in terms of the rotation angles relative to the real axis one-cycle, (1,0), 
wrapped by the 06-planes 

^j, = Trej, = arctan(^) = arctan(^), [x' ^ jfj = ^l- (n.7) 

We use conventions in which the angles between Dp-hranes and the 06-plane, 9^ = ^ and the -Dpa/-Dpfe-branes, 
^ab — H ~^a^ ^^^y inside the range, 9^ G [—1,1], with the positive angles associated with counterclockwise rotations. 
Our determination of the D-brane-orientifold angle is related to the angle determination given by the inverse-tangent 

function, 9^ G [— ^, ^1, by the identification: 9^ — 9^ for nf, > and 9j, = -I- 9i for nl < 0. Furthermore, 

our determination of the interbrane angles is related to the angle differences, 9^i^ — ^1 ~ ^a' identification: 

9ab = 9ab, if \9ab\ < 1, Or 9ab = "2 + 9ab, H \9ab\ > 1- 

The Dirac-Born-Infeld world volume action gives the mass of single i'p-branes as the product of the tension 

parameter Tp by the wrapped cycle volume, = TpT/^^ = Tp YYi^i using the definitions of Tp given in Eq. (IV. 2) 
below. This suggests that the construction of energetically stable configurations of £)-branes should consider the 



+ (-mi+2U[ni)m. (II.5) 
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cycles of minimal volume. For the CY complex rf-folds, equipped with the metric tensor Gj j preserving the complex 
structure j/ = i<5/j, there exist two types of volume minimizing sub-manifolds, which correspond to the sets of two- 
cycles and d-cycles [90] calibrated by the Kahler and holomorphic volume forms, J = iGjjdX^ AdX"^ and O^^O) ^o,d, 
respectively. Because of the relation linking these forms, ^J''' = (— 1) ^ i'^^dfi A ^o,dj and the reality condition on 
the manifold volume, Vx ^ ^ Jx '^'^ ^ ^^'^ holomorphic form arise as the one-parameter family, e^'^^^x)-, parameterized 
by the angle </?. One then defines [91] the d-cycles 11^ calibrated with respect to e^'^ild.o as the Lagrangian sub- 
manifolds (with a vanishing restriction of the Kahler (1, 1) form, J\n.^ = 0) on which the holomorphic rf-form obeys 
the reality condition, 3(e*T23.o)|n = 0. The calibrated or special Lagrangian (sLag) sub-manifolds are defined by 
embedding maps which obey first order differential equations expressing the preservation of supersymmetry charges. 
These cycles have the property that their volume integral, Vn^ = /n 5^(e"''f2d,o), is minimized among the elements 
belonging to the same homology class [11^]. From the action of the antiholomorphic reflection on the covariantly 
constant forms, TZ : J ^ — J7, e^^'^fio.d, it follows that the orientifold Op-hypcrplanes, as fixed point loci of 

7?., must wrap the sLag cycles. In order to construct a supersymmetry preserving open string sector, one must then 
consider setups of Dp-branes which wrap the sLag cycles. For the factorizable tori, where J = idX^ A dX^ and 
^^3,0 = dX^ A dX"^ A dX^, the sLag cycles intersect at angles describing ^[/(S) rotations, in such a way such that 
the brane-orientifold intersection angles, 0^, defined by Eq. (II. 7), or the interbrane intersection angles 0^^, = ^ ^/d 
obey conditions of form, ±.Q^ = mod 2. For completeness, we note that the supersymmetry conditions for the 
M-theory intersecting branes have been discussed in Ref. [92]. 

The RR charge cancellation means the absence of RR tadpole divergences in the open strings one-loop vacuum 
amplitude. This condition suffices to guarantee that the world brane effective field theory is anomaly free. For 
general setups consisting of K stacks of parallel A?'^ £>6^-branes and their orientifold mirrors, {D&^ + D6^'), [(/U, /x') = 
1, - ■ ■ , K = a, 6, • • •] the RR tadpole cancellation condition requires that the sum over all the wrapped three-cycles 

belongs to the trivial homology class, ^^(1^^*] + [n^']) + Qo6[no6] = 0. Here, the Op-plane charge is given 

by, Qop = T2^~'*.fp, where fp = 2^~p counts the number of distinct Op-planes and the =F sign is correlated with 
the (orthogonal or symplectic gauge group) orientifold projection condition on the Chan-Paton (CP) matrices and 
with the sign of the Op-plane tension parameter, rop = =F2^~^Tp, as we discuss in the next paragraph. For toroidal 
orientifolds, using the relation, Qo6[no6] = T32nj[a^], along with the decomposition in Eq. (II. 5), translates the 
RR tadpole cancellation condition into the four equations for the wrapping numbers 

^ N^yy^ ^ 16 = 0, ^ N^nj,m-;,rh^ = 0, [I J K] (II.8) 

where, to avoid double counting, one must exclude the orientifold image branes from the above summations over 
brane stacks. For definiteness, we develop the following discussion in the case of orientifolds with negative tension 
and RR charge, Qop < 0, corresponding to the SO group type projection. 

The open string sectors are associated with the distinct pairs of _Dp-branes supporting the end points. In orientifolds, 
the diagonal and non-diagonal sectors include the pairs, (a, a), (a, a') and (a, 6), (a, 6'), with the equivalence relations 
between mirror sectors, (a, a) ~ {a', a'), {a,b) ~ {b',a'), {a,b') ~ {b,a'). The gauge or Chan-Paton (CP) factors 
are time independent wave functions, described for the diagonal and non-diagonal sectors by (2A^a x 2Afa) and 

{2Na + 2Nf,) X {2Na + 2Nb) matrices, (A^''''')ij and (A'-^'*'-')^, with the labels i, j nmning over members of the 

Dp-brane stacks and the labels A running over components of the gauge group representations. The matrices x[^'°''' 

decompose into 4 blocks of size Na x Na and the matrices A^**' into 8 diagonal blocks of size, NaX Na, N^x Nb, and 
8 non-diagonal blocks of size, Na x A^5. Since the modes in the conjugate sectors, (6, a) ^ (a, b^ , have opposite signs 
intersection numbers, lab = —ha, opposite helicities and conjugate gauge group representations, these are assigned 
the Hermitean conjugate matrices, A^**'"^ = A^"'*'-*^. We omit writing henceforth the upper suffix labels on the CP 
matrices specifying the sectors. The normalization and closure sum for the CP matrices are described in consistent 
conventions as 

Trace(AAA^) = '^Trace{0iXA)Trace{02XA) = Trace{0i02), (II.9) 

A 

where the summation extends over the complete set of states A in the gauge group representation. For the unitary 
group, U{Na), the completeness sum over the adjoint group representation uses the identity, '^Ai-^A)ij{XA)ki = SuSjk- 
The orientifold symmetry fllZ is embedded in the gauge group space of a Dp-brane setup through a unitary twist 
matrix, jmz, by imposing the projection condition, A^ = — TotcA^Tq^- convenient construct for 707^ is given 
by the direct product, jqh = ®fj,"/mi,fi- The anomaly cancellation constraint commonly imposes the tracelessness 
condition, Trace{'jcnz) = 0, along with the symmetry conditions, Jq-;^ = ±7^7^, in correspondence with the SO 
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and U Sp type projections. In the special case of a D6a-brane stack overlapping the 06-planc, hence coinciding 
with the mirror image I?6a'-brane, the gauge symmetry in the SO and U Sp projections is enhanced to the rank Ng, 
orthogonal or symplcctic groups, SO{2Na), USp{2Na), in correspondence with the negative and positive signs of the 
06-plane RR charge^ Qoe and tension parameter. To detail the construction of CP factors, wc consider, for the sake 

of illustration, the (a, a) + (a, a') sector. The SO type projection matrix, 'ymi,a ~ (^j^ ^ ' yields the SO{2Na) 

group adjoint representation, X^^'^^^ = | "^^t ) > with the restriction to ai = 02 = yielding the U{Na) group 

adjoint representation. We here use conventions where the symbols m, s and a for block entries designate generic, 
symmetric and antisymmetric matrices, respectively. The antisymmetric representations of U{na) (in the SO type 

projection) are realized in the {a, a') sector by the 2Na x 2Na matrix solution, A^'^^l ~ (o o)' ^'^^ conjugate, 

^[■A-lt Tijg bifundamental modes in the non-diagonal (a, 6) sectors are realized by the {2Na + 2Nb) x {2Na + 2Nb) 
matrix solutions 

A«-^-.i = (» ^),[B=(» ?).^ = (:f. Sy (n.io) 

obtained by setting successively the Na x TVf, block entries, a, (3, 7, S, to non-vanishing values. We recommend 
Refs. [22, 93, 94] for a further discussion of ZJ-branes. 

Proceeding now to the gauge group, we consider first the case of a single stack of Na coincident DSa-branes and its 
mirror DSo'-brane stack wrapped around three-cycles at generic angles 6^ and 9^, relative to the 06-plane, hence not 
overlapping the 06-planes. The massless states in the diagonal sectors, (a, a) ^ (a', a'), include the gauge bosons of 
the gauge group U{Na), along with adjoint representation matter modes. For the pair of intersecting D6a/l?6(,-brane 
stacks and their mirror image D6a' / D6h'-branes carrying the gauge symmetry group, U (Na) x U (Nf,), the non-diagonal 
sectors consist of conjugate pairs, [{a,b) + {b,a)] ^ [(b' ,a') + (a' ,b')] and [(a, &') + (&', a)] ^ [(6, a') + (a', 6)], with localized 
(improperly named twisted) states carrying the bi- fundamental representations, Iab{Na, N(,) and Iab'{Na, N(,), with 
multiplicities given by the wrapped cycles intersection numbers 

= ]J{niml - nlrhi), /ah' = ]J(n^m^, - n^,m^) = J| -(n^m^ + n^m^). (H.ll) 
/ / / 

The sectors {a, a') have a total number of intersection points, laa' = Ylii~'^i^i'>T^i)^ of which the I^'^J^' ~ 74 °°^ 

ill "a 

points, symmetric under the reflection TZ, give rise to modes carrying (in the SO type projection) the antisymmetric 

representation of the gauge group U{Na), X = —X^, while the remaining modes split into pairs of modes carrying 

the symmetric and antisymmetric representations with the same multiplicity, -/^^^^'^^ = |-^aa'(l — Yl^ n' ^' "^^^ 
multiplicities of the symmetric and antisymmetric representations are thus given by, /^^^ = ^/aa'(l ± "rr~T)- 

We follow the familiar description of fermion and boson modes in terms of the basis of left chirality states, (/l, /£), 
where the right chirality states are obtained by applying the complex conjugation operator exchanging particles with 
antiparticles, ^ f1 = (/£)^, //j ^ Jl = Note that the correspondence relations for the electrowcak SU{2) group 
doublets include extra signs, with, for instance, the quark doublet fields given by, Jl = {ul, dL), fn = {ur, da), /£ = 
{d'l, —u'l). For the spectrum of modes with the left-right chiral asymmetries, labfL and /cd/E, the presence of Aab 
and Acd mirror vector pairs results in the non-chiral spectrum, (lab + Aab).fL + ^abfu, {led + Ac(j)./£ + ^cdfR- Going 
from positive to negative intersection numbers entails changing the sign of the chirality (helicity for massless fermions) 
and conjugating the gauge group representations. For instance, the massless fermions with negative multiplicities 
lab, led: would rcfcr to right chirality fermions (or left chirality antifermions) carrying the conjugate bi- fundamental 
representations, \Iab\{Na, Nb), \Icd\{Nc, Nd)- 



3. First quantization formalism 

We only discuss here the non-diagonal open string sectors {a,b). The coordinate and spinor field components 
along the flat M4 space-time dimensions obey the N conditions, daX^^ = 0, ip^^ ^ i(;f^ = 0, at both end points, 
a = (0, tt), where the upper and lower signs apply to the R and NS sectors. For DGa/i'Ofj-brane pairs wrapped at 
the angles (/"^ 5 = ttO^ ^ in Tf, the rotated complexified coordinate components, e^*'^".''X^, e^*'^" '' (-0 T split into 
real and imaginary parts, longitudinal and transverse to the branes, hence obeying the N and D boundary conditions, 
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5^3?(e"*'^°.''X^) = and dt'^{e~''^-.''X^) = 0. The corresponding conditions for the rotated spinor field components 

read, JR(e~*'^"''' (V'^ T V'^)) = and ^{e~^'^'''^{ip^ ±ip^)) = 0, with the upper and lower signs referring to the R and NS 
sectors. In terms of the complex plane z, z variables, the N and D boundary conditions along the real axis for the 
pair of ^ rotated £>-branes read in full as 



{d + d)[e-''<''^''' BX' - e''<'i»dX'] = 0, e-''''i»{i)^ ± i>') - e'<''{'4)^ ± V^^) = 0, (11.12) 



where the labels [a, 6] correspond to the open string end points, a = [0, tt]. For convenience, we extend the notation 

for the intcrbranc angles to, = {S^abT^'ab)^ "^itb the understanding that = in our present discussions. In 
the light cone gauge of the 2-d world sheet superconformal field theory, the open string states include the quantized 
oscillator modes described by the number operators, Nx, N^, and the zero modes described by the momentum and 
winding quantum niunlxTS, p^,^, ,s^;^ G Z, and by the intcrbranc transverse distances, Y^j^f^, for the complex directions 
along which the braiics arc parallel. The string oscillation frequencies along the complex dimensions have integral 
modings shifted by the I^-brane angles, rij. = n" ± ^^"j, [n° E Z, a = 1,2, 3, 4]. In the boson representation of spinor 
fields, the oscillator number operators, N^a, are replaced by the i?" fields momentum vectors, = (r-^,r^), [/ = 
1,2,3] corresponding to the weight vectors for the Lorentz group SO{8) ~ Spin{8) Cartan torus lattice. The GSO 
projection for the world sheet fermion number parity symmetry, (—1)^, restricts the weight vectors to the sub-lattice, 
r" e Z + ^ + [X]a=i ^" "= 2^ + 1] where the boson and fermion (NS and R) sector modes with = ^, 0, are assigned 
vector and spinor type weight vectors. The fourth component of the 5*0(8) weight vectors, = 0, ±1,---, 
describes the chirality (or helicity quantum number for massless fermions) for the SO {2) little group of the flat space- 
time Lorentz group. The remaining three components, r^, describe the helicity quantum numbers for the three 50(2)7 
subgroup factors of 50(8). The quantized string mass shell condition for the (a, b) sector is expressed by the general 
formula for the string squared mass spectrum 



2' 

7=1 a=l a=l 7=1 



(11.13) 



where we continue using the conventional range for the brane intersection angles, — 1 < ^^(, < +1. The last two terms 

in the squared mass, a' A4^f^, describe the string zero energy vacuum contributions, while the first term, Mj^^^, given 
explicitly in the second line entry, includes the contribution from the momentum and winding string modes and from 
the transverse separation distance F/^ of the D6a/ Ddb-hianes along the Tf where they are parallel. The latter point 
is reminded by the symbol (5^^ which is non- vanishing whenever the displaced i36a/D666-branes are parallel along 
some complex plane / so that O^i^ = 0. 

We now explicitly describe the low lying string modes. The solutions for massless spin-half fermion modes select the 
unique conjugate pair of spinor weights, r° = (r^, r^) = ±(— i, — i, — i, i), with the overall ± signs corresponding to 
the two possible spatial helicitics. The abbreviated notation for the spinor weights illustrated by, ±(— i, — i, — i, i) = 

±( , +), will be adopted for convenience. The scalar modes of smallest squared mass select the four solutions for 

the vector weights, r" = ±(—1,0,0,0), ±(— 1, — 1, — 1, 0), with the underline symbol standing for the three distinct 
permutations of the entries. The resulting four solutions enter with the squared masses 



li-Ktl + + \olb\), k+\olb\ - Kbl + \olb\ 



^ (11.14) 



The lowest lying vector boson mode arises from the vector weight, ±(0, 0, 0, 1), with squared mass, M^^ = ^\0i^\. For 
completeness, we note that the towers of so-called Regge resonance gonion modes [51] of scalar and vector boson types 

correspond to the oscillator excited states, V"! '''kV'-r+Q:-ra+)™^ |0 >ns, with mass squared, M^^ = {'m^T^)\^ib\, [m^ S 
Z]. We recommend Rcf. [95, 96] for a further discussion of the mass spectrum in intersecting brane models. 

In parallel with the closed string geometric moduli, there arise open string sector moduli which correspond to order 
parameters of the world brane gauge field theory associated with the D-branes positions and orientations. Thus, 
the transverse coordinates of a D6a-brane stack are moduli fields in the adjoint representation of the U{Na) gauge 
theory which parameterize its Coulomb branch deformation. The recombination of a pair of intersecting branes into 
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a single brane, a + 5 — > e, or the rcconnection of two branes, a + 6 — > c + d, are described in terms of the moduli fields 
in bi-fundamental representations of the (a, 6) sector which parameterize the Higgs branch of the U(Na) x U{Ni,) 
gauge theory. The brane splitting fixes the VEVs of open string moduli while the brane recombination redefines the 
VEVs of open string localized moduli needed to avoid the vacuum instability from tachyon modes, in analogy with 
the Higgs gauge symmetry breaking mechanism. The splitting and recombination processes are accompanied by mass 
generation mechanisms which decouple pairs of fermion modes in vector and chiral representations. Representative 
examples of these deformations are the Higgs mechanisms for the unified and the electroweak gauge symmetries. 
The consistent description of D-brane recombination using non-factorisable cycles [53, 73, 97] does indeed lead to 
a reduction of the wrapped cycles volume and of the fermion spectrum chiral asymmetry, in agreement with the 
Higgs mechanism. In spite of the poor information on string non-perturbative dynamics, interesting results have been 
established concerning the existence of bound states for Dp/D{p + 4)-brane pairs and for Dp/D{p + 2)-brane pairs in 
backgrounds involving NSNS or magnetic field fluxes [33, 98, 99] or the T-dual backgrounds of Dpa/ Dpi,-hranc pairs 
wrapping intersecting cycles [100]. We also note that the recombination process can be partially formulated in the 
context of branes realized as gauge theory solitons [95, 96, 101]. 



4- Conformal field theory formalism 

The conformal field theory provides a powerful approach to calculate the on-shell string S-matrix in perturbation 
theory. The open string amplitudes are obtained by integrating the vacuum correlation functions of the modes vertex 
operators inserted on the world sheet boundary. We focus here on the tree level amplitudes of the (a, h) non-diagonal 
sectors of the D-hvajoe pair, D6a/D6i„ intersecting at the angles, (/>^j — (pl — cj)^^. With the field doubling prescription, 
the world sheet field propagators are simply given by 

(^(zi)(^(z2) >= - ln(zi2), < H^{z^W{z2) >= -5^^ ln(zi2), (n.l5) 

where Zi2 = Zi — Z2- Since the coordinate and spinor field components of M4 obey the N boundary conditions, 
[d — d)X^ = 0, one can formally replace the Minkowski space coordinate field components along the complex plane 
real axis boundary as, [X'^{z) + X'^{z)] 2X'^{x). The insertion of the open string mode (a.b) at the real axis 
point, Xi = 3^(2;^), modifies the boundary conditions on the right hand half axis, x > Xi, in such a way that the two 
orthogonal linear combinations, associated with the real and imaginary parts of the rotated complexified coordinate 
fields, obey the N and D boundary conditions: ^(e'^'^'-od^X^) = 0, ^{e' ip ^i)) = Q, and ^(e'^'^'-odtX) = 
0, 9(6"*''^°'' (V'''^ ± tp^)) = 0, where the upper and lower signs refer to the R and NS sectors. The left and right half 
lines, X G [—00, Xi] and x G [xi, 00], are mapped to the D6a- and -DSfc-branes with the boundary conditions given by 
Eq. (11.12). Since only the interbrane angle really matters, the boundary conditions on the coordinate and spinor 
field combinations along Tf can be expressed by the same formulas as in Eq. (11.12) with 0^ (, ^ (j^ib = 4'b~ 0o5 along 
the half line x £ [xi,oo], and ^ — > along the half line x G [— 00, a;^]. Taking the sum and difference of the two 
relations yields the equivalent form of the boundary conditions 

dX^ - e^^'^-f = 0, dX^ - e^''t'-^dX^ = 0, 

ij,! ip e^^'I'-^'i)' = 0, i)^ ± e^'^'-oi)^ = 0. (n.i6) 

Note that our sign convention for the brane angles is opposite to that used in Refs. [22, 36] and that we differ from 
Ref. [36] in certain relative signs. 

We now discuss the covariant conformal gauge formalism of the world sheet theory. Each open string state of the 
non-diagonal sector, C G {a,b), is assigned a primary vertex operator of ghost charge q and unit conformal weight, 

V(^'^^^ j^j(zi, fcc, Ac), with kc denoting the incoming four momentum and Ac the gauge wave function factor. The 
building blocks in constructing the vertex operators are the coordinate fields, X^'"^{z), their derivatives, dX^'^{z), 
and exponential maps, £^''c-x{z,z). ^.j^g spinor fields, ip^'^{z) = e'^'*''*(^); the superconformal ghost scalar field ip{z) 
exponential map, e^'^'^^ of ghost charge q; the spin and twist operators for spinor fields along the flat space-time 
and internal space directions, Sr<^{z) = e*'' ^^'^-^ and .s!j_gj(z), [r" = {r^,r^), = (r-'^, r^, r^)]; and the twist 

operators for coordinate fields along the internal space directions, a±gi{z). The weight vectors, = (r^,r°), denote 
the momentum vectors of the complex scalar fields, H^{z) = [H°'{z), H^{z)], belonging to the Cartan torus lattice 

of the Spin{10) Lorentz group. The twist and spin operator factors are needed to produce the requisite branch point 
singularities at the modes insertion points. These operators create the ground states of the twisted sectors upon acting 



10 



on the SL{2, R) invariant ground states of the NS and R sectors. For the low lying non-diagonal sector modes with 
excited coordinate oscillator states along the internal space directions, alongside with the ground state twist field, 
(T'^giiz), one needs to introduce the excited twist field operators, T±gi{z), f±gi{z). The spinor field ground state and 
low lying excited twist field twist operators, s^gi{z), f^giiz), i'^gi{z), are explicitly realized by the free field vertex 
operators, 

(Z) = e^^(0'+^">"'(^\ = e±^i9'+r' + l)HJi.)^ ^ ^ ^±,(,^+,^_l)ff,(,) ^ ^^^^^^ 

labeled by the angle 6^ and the 50(6) Lorentz group weight vector, . The GSO projection for the world sheet 
fermion number parity symmetry, (—1)^, correlates the weight vectors for the flat space-time 5*0(1,3) ~ 50(4) ~ 
SO{2) X SO{2) helicity, r" = (r'',r°), with those of the internal 50(6) helicity, r' = (r^r^, ?-'^). For the R sector 
fermions, this requires the number of — ^ entries in the five-component spinor weight vectors, = (r^,r"), to have 
a fixed parity (odd in our conventions). The left and right hehcity (chirality) fermions are thus described by the 
spin operators, S^-a = e"' ^» . with weights: = r° = ±| and = — r° = respectively. Note that the 

50(10) group weight vectors, — {r^ ,r"-) = (r'',r°), reduce in the light cone gauge to the 50(8) group weight 
vectors, r°' = {r\r*). The same description applies to the (a, a) diagonal open string sectors upon introducing the 
spin fields and the spinor twist fields, s'" (2;) = e"' ^ with vanishing angles. To develop a unified formalism 
for both the diagonal and non-diagonal sectors, we adopt the self-explanatory notation for the twist operators, 
s^^0a{z) = e^*'^ '^^^ encompassing the case, 0^ = 0. 

Unlike the spinor field twist operators, the coordinate field twist operators do not have a free field representation. 
An implicit definition can still be obtained by specifying the leading branch point singularities in the operator product 
expansions of these operators with the primary operators constructed from the coordinate fields 

dX\zi)aei{z2) ~ z[l ^^^9/(22) H , dX\zi)a_e'{z2) ~ "^^^-9^(22) H , 

dX'izi)cT_gl{z2) ^ Z^fi_gjiz2) + ■■■, dX ' {zi)a +gl {Z2) ^ Z^f f+gj iz2) + ■■■, (11.18) 

where T±e, f±0 are the excited state twist field operators introduced earlier, and the dots denote contributions 
operators which are regular in the limit 012 = zi — Z2 ^ 0. We have used here the abbreviated notation for the brane 
angles, 6*^^ — > ±6*^, [9^ G [0, 1]] with the sign made explicit in such a way that the results for the negative and positive 
brane intersection angles, , are related by the substitution, 9^ 1 — 9^ . For completeness, we also quote the 
operator product expansions for the spinor twist field operators in terms of the excited twist operators introduced 
above, 

i^'{z^)slliZ2) ~ zil'+'^'kl\{z2) + ■■■, i^'{Zl)s':!,,{Z2) ~ zg' k^!,, {Z2) + ■■■, 

i^'izi)sL0r{Z2) ~ Z^2^''+'^'n:l,riZ2) + ■■■, (z^) S^'r {Z2) ~ z;2^'' ^'^^ % {Z2) + ■■■ (11.19) 

It is of interest to note that the singular dependence on the brane angles cancels out in the operator product expansions 
of the coordinate and spinor twist fields with the energy-momentum and supersymmetry generators, 

T{z) = -GMNi^dx^dx'' + iv'^a^^), 

a 2 

Tf{z) = i\ ^GMNi>^dX^{z) = iJ^{ip''dX^ + iP^dX^ + ijj^dX^). (11.20) 
\ a \ a 

For instance, TF(zi)(sl,agi){z2) ^ i\[^{z':[^~'' iliTgi + zf^tlifei). 

The following formulas are of use in evaluating the conformal weights of various operator factors, 

hia±e') = ^O'il - e'), h{T±gi) = ^e'iS - O'), h{f±e') = \{l - e'){2 + 6'). (11.21) 

The mass shell condition for a mode of mass squared, Af^, is then determined by requiring that the total conformal 

weight of the mode C vertex operator, Vc{z) = Vc{z; kc, Ac), amounts to unity, 1 = h{Vc) = + ■ ■ ■ = —M^ H . 

The vertex operators take different forms depending on the superconformal ghost charge [22], q G Z + ly + ^, 
carried by the scalar ghost field exponential, e^^'^-', with v = ^,0 in the NS and R sectors. The canonical pic- 
tures (unintegrated with respect to the superspace variable, 9) involve the superconformal scalar ghost field fac- 
tors, V^~^\z) ^ e^'^^^''0^^^\z) and V(~^)(z) ^ e~^2— O^^^) (z), whereas the integrated (with respect to 9) ver- 
tex operators, of higher superconformal ghost charges, are obtained by acting on the canonical operators with the 
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picture changing operator, G_i = e'^'^^^Tp{z) + • • •, where the dots refer to ghost field terms. The isomorphic 
representations of the vertex operators of increasing ghost charges are obtained by the stepwise incrementation, 

= \mi^_^V{z,w)V'-'i\w) = hm^^^Tir(2;)e'^(^)y(9)(M;) H . Since the vacuum for the world sheet surface 

of genus g carries the defect ghost charge, (2g — 2), in order to conserve the ghost charge in the vacuum correlator 
involving np and ub fermion and boson vertex operators carrying the natural charges, — |, —1, one must apply the 
picture changing operator (PCO) on the number of vertex operator factors, Npco — riB + ^ + 2g — 2. For instance, 
the four point vacuum correlators on the disk surface require, Npco = + ^ — 2, while those on the annulus 
surface require, Npco = np + 

We are now ready to complete the construction of vertex operators. The matter and gauge boson modes are 
described, in the diagonal sectors of parallel Z>6-branes, by the following vertex operators in the canonical and once- 
derived ghost pictures, with charges q = —1,0 for bosons and q= —5, +5 for fermions, 

• V^-;^\z) = Ac/e-'^V^e''=-^, V^°^{z) = i^Xci[dX' - ia'{k ■ ^)V'^]e*'=-^, 

^ ^ \ a' 

^ 1 

. Vi-'\z) = XA^e-'^ea^{k)re"'-'', V^\z) = i\[^XA^Mk)[dX^ - ia'{k ■ i^WV'^''' , 



(11.22) 



where, A^/, A^a, denote the CP factors, k the 4-d momenta, u(k), e^(fc) the Dirac spinor and polarization vector 
wave functions for spin i, 1 particles, and the suffix labels v, s in r^^ are used to remind ourselves that boson and 
fermion modes carry vector and spinor 50(10) group weight vectors. In the non-diagonal ('twisted') sectors, (a, 5), the 
boson and fermion mode vertex operators in the canonical and once-derived ghost pictures are given by the following 
formulas: 



J i^j 



e""-^ , (11.23) 



(-) 

where the dots in V^^^ (z) represent 0{k) terms of complicated form that we shall not need in the sequel. We 
recommend Refs. [102, 103] for further discussions of the vertex operator construction. 

The processes of interest to us in this work involve four massless fermions belonging to two same or distinct pairs of 
conjugate modes, f — f, 9 = 9' and f ^ f , 9 ^ 9'. In the Polyakov functional integral formalism for the string world 
sheet, the n-point open string tree amplitudes are represented by the disk surface punctured by n points a;, inserted 
on the boundary with ghost charge obeying the selection rule, J27=i 1i ~ Since the unpunctured disk surface has 
no moduli, the integration over the moduli space consists of integrals over the ordered real variables, Xi G R, summed 
over their cyclically inequivalent permutations and divided by the Mobius symmetry group, SL{2,R), generated by 
the homography transformations of the disk boundary. Following the familiar Faddeev-Popov procedure of gauge 
fixing and division by the volume of the conformal Killing vectors (CKV) group, one can write the four-point tree 
string amplitude, A4 = A{f{ki)p{k2)f'{ks)f ^(^4)) as the integral of the vertex operators vacuum correlator, 

(11.24) 



where we follow the familiar convention in which all of the particle quantum numbers are incoming. The elements of 
the permutation group quotient, a G S4/C4, consist of the three pairs of direct and reverse orientation permutations. 
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and the integrations are carried over the ordered sequences of the Xi. The invariance under the SL{2, R) subgroup of 
the conformal group is used to fix three of the insertion points, say, at the values, Xi = 0, x^ = 1, x^ = X ^ oo, with 
the free variable varying inside the interval, X2 ^ x E [— oo, oo], so as to cover the three pairs of cyclically incquivalent 
permutations, and Vckg = '^/X'^. We have labelled the open string vertex operators in Eq. (11.24) by the pairs 
of associated branes, such that the disk surface is mapped in the Tf complex planes on closed four-polygons whose 
sides are parameterized by the linear combinations of N coordinates tracing the equations for the adjacent branes, 
D, A, B, C. This map is illustrated in Figure 8 of Appendix A. With the world sheet boundary represented by the 
complex plane real axis, the reference ordering of insertion points for the trivial permutation, cr = 1, determines the 
four segments, (— oo,a;i), {xi,X2), (0:2,2:3), {x3,X4 +00), on which the orthogonal linear combinations of internal 
coordinate fields, ^{e~^'^'^>>X^) and 3(e~*'^»!'X'^), obey N and D boundary conditions, with c/)^;, denoting the fixed 
interbranes angles at the corresponding insertion points, as displayed in Eq. (A. 10). The four point string amplitude 
may thus be represented by the sum of three reduced amplitudes, 

E/ 
CJ"G04/ C4 

xX^J^ °°dx< I>i7^(x.jf|7^(x.jfl7\a;,3)vi7\a;,j > 
= [^(1234) + (2 ^ 4)] + [2 ^ 3] + [1 ^ 2] 

= [^(1234) + ^(4321)] + [^(1324) + A(4231)] + [A(1342) + A(2431)], (11.25) 

where we have denoted, —6 — 6a — 0d = Sa — Ob, —d' = 0c — &b = &c ^ ^d, and introduced the hat symbol to denote 
the vertex operators with the CP matrix factor removed, V^'^^Xi) = V^'^^Xi)Xi. The factor X"^ from the gauge fixing 
cancels out with the X-dependent contributions from the correlator. With the incoming flat space-time four-momenta 
denoted by fc,, obeying the conservation law, fci + fc2 + fcs + fc4 = 0, one can express the Lorentz invariant Mandelstam 
kinematic variables as, s = — (fci + fc2)^, t = — (fc2 + ^3)^, u = — (fci + ^3)^. A compact representation of the reduced 
amplitudes may be obtained by considering the definition of the correlator with the dependence on the kinematic 
invariants extracted out, 

<V^'s^\mU''\^)vl.~e,^\m^^^^^ (11.26) 

while rewriting the second and third reduced amplitudes, ^4(1324) and A(1342), in terms of the first reduced amplitude, 
A(1234), through the change of integration variables, x € [l,oo] ^ x' = € [0, 1] and x € [— oo,0] x' = G 
[0, 1]. These steps lead to the compact representation of the disk level string amplitude 

^(1234) = [rrace(AiA2A3A4) + Trace(A4A3A2Ai)] / cixx-^-^il - x)-*-^Ci234{x) 

L Jo 

/■I 1 

-F[Trace(AiA3A2A4) + Trace( A4A2 A3 Ai)] / dxx-*-\l - a;)-"Ci324(-; ) 

+ [Trace{X2XiX3X4)+Trace{X4X3XiX2)] [ dxx'" {1 - x)-''-^Ci342{^^^—^) ■ 

Jo ^ 



(11.27) 



B. String amplitudes from world sheet correlators 

We discuss here some practical details of use in evaluating the open string amplitudes for the configuration of 

non-diagonal sector modes involved in Eq. (11.24) for the amplitude -4(/i/2/3/4^). Since the ordering of adjacent 
DG-branes is determined by that of the vertex operator insertion points, Xi, we deduce by simple inspection that only 
the direct and reverse orientation permutation terms for the first reduced amplitude, A(1234) + A(4321), is allowed, 
while the other two pairs of reduced amplitudes are forbidden. Only for symmetric configurations involving subsets 
of identical Z)-branes, do exceptions to this rule occur. 

The correlators receive contributions from three sources. There are first the quantum or oscillator terms coming from 
the Wick pair contractions of free field operators, which arc determined by the world sheet field propagators. The sec- 
ond source is associated with the CP factors which arc grouped inside traces of ordered products. The third source is as- 
sociated with the classical action factor in the functional integral which accounts for the string momentum and winding 
zero modes for the coordinate field components along the compact directions. The heaviest calculational task resides in 
the coordinate twist field correlator factor. The correlation function, Zi{xi) =< a_gi{xi)agi{x2)<7_g'i{x3)ag'i{x4) >, 
is evaluated by making use of the stress energy source approach initiated by Dixon et al., [25] and Bcrshadsky and 
Radul [24]. One expresses the constraints from operator product expansions, holomorphy and boundary conditions 
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on the two correlators, g{z,'w;xi), h{z,w;Xi), obtained from Z^{xi) by inserting the bilocal operators, dzX^ d^X^ 
and d^X^ dyjX^ . The resulting formula for (xi) comprises two factors including the contributions from quantum 
(oscillator) and classical (zero mode) terms, Z^{xi) = ^^^(a^OSc; ■^c;(*^Oi where the classical summation is over the 
lattice generated by the closed 4-polygons with sides along the branes, A, B, C, D. We have found it useful to provide 
in Appendix A a comprehensive discussion of the correlators of open string modes involving distinct angles, 0^0', 
since this application has not been addressed in great detail in the literature. Our presentation there closely parallels 
that of Biirwick et al., [29] for the closed string orbifolds. 

Two important constraints follow upon requiring that the world sheet boundary is embedded on closed polygons in 
the Tj planes. For the coordinate twist field correlator, < <J0iji{xi)(Tg^j^{x2)o'e^j^{x7i)ag^j^{x4) >, the closed four- 
polygons have edges along the N directions of the iP-branes, with vertices and angles Oi € [0, 1] identified to the 
intersection points and angles ff and Oi of the adjacent branes. We use here the index a to label the intersection points 
and the notational convention for the angles, 9i = [9i, l—9i] for ±(?i. The first condition expresses the obvious geometric 
property of the angles, Ylt=i H ~ 2- The second condition is related to the consistent configuration for the intersections 
of the various branes pairs. Following the initial discussion for 3-point couplings by Cremades et al., [54], a general 
comprehensive discussion of this problem was provided by Higaki et al., [104], whose presentation is closely followed 
here. We start by observing that each pair of branes a = A, B intersect at Iab = nArriB — nsmA points lying along 
the branes A, B with coordinates, Xa = ^° ^° , labelled by the integers Ua, ks G [0, 1, • • • , Iab — 1], such that each 
intersection point is associated with a unique choice for the pair of integers kA, ks- In the case of branes intersecting 
at the origin, solving the complex linear equation, Xa — X^-, \Xfx — ^al'a ~1~ Qa^i ~^ Pa^2i Lot — ^a^i H~ ^a^2) ^ 
-R; (Pa)9a) S Z] yields the explicit representation for the integer parameters, ka = UaPAB — "maQAB, \pab = 
PA — Pb, Qab = qA — Qb]- Since the intersection points are defined modulo the grand lattice, Aab, generated 
by La, Lb, they form equivalence classes defined modulo the addition of vectors of A^^. A convenient way to 
characterize these Iab classes is in terms of the shift vectors, wab = Xa — Xb = ^^tZb' ~ ^Iab ' s-ssociated to the 
choices of integers kA, kg appropriate to the various intersection points. Since the vectors wab belong to the torus 
lattice A, generated by the cycles ei, 62, and are defined modulo Aab, they arise as the independent elements of the 
lattice coset. A/ Aab- One can also interpret the shift vectors as the A lattice translations which bring the intersection 
points on branes A, B in coincidence, or equivalently, as the segments linking the open string end points located on 
the branes A, B. For the 4-point correlator with the configuration of adjacent branes, ABCD, the condition that 
the 4-polygon closes, may now be expressed by the selection rule involving the shift vectors associated to the four 
adjacent brane pairs, wda + wab +wbc + wcd = modulo A. While the Iab independent classes of shift vectors are 
in one-to-one correspondence with the intersection points, they do not specify the coordinates of these points which 
must be calciilatcd independently. Higaki et al., [104] have given a simple useful procedure to explicitly evaluate 
the shift vectors. One starts by testing whether the winding numbers of the A, B brane pair along the two lattice 
cycles are relative primes, by considering their greatest common divisors (gcd) defined by, gcd{nA,nB) = Nab and 
gcd{m,A, ms) = Mab- The independent set of shift vectors is then given by wab = Pab£2, [pab = 0, 1, • • • , Iab — 1] 
if Nab = 1, or by wab = QabGi, [qab = 0,1,- ■■, Iab - 1] if Mab = 1- Otherwise, for Nab 7^ 1, Mab ^ 1, 
the independent set of shift vectors can be chosen as, wab = QABei + PABG2, [pab = 0,1, ■■■,Mab — 1, (Iab = 
0, 1, • • • , Ij^g — 1]. The above rules readily generalize to the case of n-point correlator, < W^^i Ugi (^iAi+i) f°i i^i) > 
, [^ra+i = Al] where the requirement that the embedding n— polygons, A1A2 ■ • ■ An, close in each Tf, is expressed by 
the selection rules on the angles and the shift vectors [104], 

n 

J2sl = {n-2), u;i,A.+w'iA3 + --- + w'A„A. =0. (11.28) 

These results hold irrespective of whether the branes Ai intersect at a common point, chosen above as the origin 
of the coordinate system. Finally, we observe that there exist a close formal similarity with the shift vectors and 
fixed points, Wh, introduced in T"^ jZjq orbifolds with lattice A symmetric under the point group rotations, 

d^, [/i = 0, 1, • • • , A^] by using the definition, (1 — ■(?'')(/" + A) = . The shift vectors described by the lattice 
translations which bring the corresponding fixed points in coincidence with themselves after applying the rotation 1?'*, 
arise here as the representative elements of the lattice coset, A/(l — i9'')A. However, it is important to realize that for 
the open strings in intersecting brane models, in contrast to the closed strings in orbifold models, the selection rules 
have nothing to do with the point and space group symmetries of the torus lattice. 

For the four-point string amplitude in the equal angle case, 9 = 9', the embedding 4-polygon is a parallelogram, 
so that the selection rule takes the simple geometric form, /2 — /i + ./'4 — /a = 0, in terms of the intersection points 
/i, • • • ,/4. Since the intersection points are naturally associated with the generation (flavor) quantum numbers of 
matter modes, one might wonder whether generation non-diagonal four fermion processes may be allowed at the tree 
level, only subject to suppression from the classical action factor. However, the combined constraints from gauge 
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symmetries and the above tree level selection rules on angles and intersection points, are seen to conserve flavor 
and hence disallow the flavor changing neutral current processes. Assuming for the sake of illustration that the 
intersection points label the quark flavors, fi{q), one indeed finds that the AS" = 1, 2 strangeness changing processes, 
s\d2d\d4 and s\d2s\d4, require the relations, /2(d) — fi{s) = and f2{d) — + .f4{d) — /3(s) = 0, which cannot 

be satisfied unless the intersection points for d, s quarks are coincident. From these observations it follows that 
the matter fermions trilinear effective Lagrangian couplings with the non-localized massless or massive boson modes 
are necessarily flavor diagonal. The quark and lepton flavor symmetries are broken only by the fermions Yukawa 
couplings with the electroweak Higgs bosons with the flavor mixing arising in the familiar way through the fermion 
mass generation mechanism. 

The 4-d space-time structure of amplitudes is strongly restricted by the symmetry constraints. The GSO projection 
correlates the helicities r^(i) along the internal and non-compact space directions (odd number of — | for fermions), 
as already noted, while the selection rules from the 50(10) Lorentz symmetry group imposes the Ha momentum 
conservation, rA{i) = 0, summed over the n modes of the correlator. The Hj momentum conservation, following 

from the symmetry under the SO{6) space group, imposes the conditions on the branes intersection angles, Y^"^i Of = 
0, [—1 < 9( < 1, I = 1,2,3] which identify with the previously quoted selection rule. These conditions often suffice 
to uniquely determine the Lorentz group covariant structure of matrix elements with respect to the Dirac spinors. 
For the configuration, fi{—6), /2(^), h{—0'), fl{0'), the restrictions on the spinor weight solutions for the massless 
localized fermions entail that only the configurations involving pairs of conjugate modes with same or distinct angles, 
±9^ are allowed, so that only the reduced amplitude >1(1234) is non vanishing. Since the fermions in the two pairs of 
conjugate states with opposite space-time chiralities require setting the weight vectors as 

-rHl) = i ,+-), r^(2) = ( ,— ), 

-r^(3) = ( , -+), r^(4) = ( , ++), (11.29) 

the Dirac spinors can only be contracted via the 10-d vectorial coupling, (uj CT mU2){uJ CT'^'^ U4) , which reduces in 
4-d to the matrix element with vector contraction, {ui^ij,U2){u3j^U4). This unique structure, up to Fierz reordering 
of Dirac spinors, is antisymmetric under all of pair permutations of the (commuting c-number) spinor factors, as 
it should be. Note that the scalar coupling of Dirac spinors would appear upon considering configurations mixing 
the fermion and antifermion modes, / and Z*^. Translating between different structure of the Dirac spinors matrix 
elements is conveniently performed by making use of the 4-d Fierz-Michel identities, given for the c-number Dirac 
spinors by 

^'il7^^*2L = W27?,7^"lfi, (Mlff7''wiff)(w3ff7M"4i/) = -{u\Hl'^UiH){U5HlfiU2H), [H = L,R] 
{Uil1^U4l){u3R1ij,U4r) = 2{uilU4r){u3RU2l) ■ (11.30) 



C. Orbifold compactification 

The covering space formalism of orbifold compactification is developed by including all the states produced by the 
orbifold group action prior to projecting on the physical states invariant under the orbifold symmetry. We restrict 
consideration to the subset of Abelian orbifolds, T^/Zpf, with the cyclic groups generated by the order TV unitary 
matrices, 9 e SU{3), [Q^ = 1] yielding Af = 2 supersymmetry in the closed string sector. The complexified bases 
of coordinate and spinor fields of the symmetric 6-d factorisable tori, — fli Tf, transform by the diagonal unitary 
matrix transformations 

[e» = diag (e^^'^s"', e^'^^s"', e^'^^s"'), 5 = 0, 1, • • • , TV - 1] (11.31) 

where the generator Q is represented by the twist vector, v = [v^], satisfying the conditions, X^/V^ = 0, A^v^ 
i I mod N. For the compactification on x T^/{Z]\r + ZnUTZ), with the orientifold point symmetry group including 
the elements, flTZQ^, [5 = 0, • • • , TV — 1] one must require that the generator QTZ acts crystallographically on the 6-d 
torus T^. This introduces conditions on the torus moduli which transform certain continuous vacuum degeneracies 
into discrete ones [60, 62, 105]. Thus, for tori, the reflection QTZ has only two inequivalent actions up to coordinate 
rescalings. The first corresponds to the diagonal reflection about one of the two torus cycles, say, e\ (Case A) and 
the second to the reflection about the diagonal sum of cycles, say, ei + 62 (Case B). An equivalent action to Case 
B corresponds to the diagonal coordinate refiection about the single cycle, ei, followed by a complex rotation, 
g»7 x^. Explicit solutions for the allowed A, B lattices have been obtained in the various Abelian orbifolds [79, 83]. 
Extensions to non- factorisable tori [81] and to exceptional cycles in orbifolds and smooth Calabi-Yau manifolds [80] 
have also been discussed in the literature. 
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The invariancc under the orientifold group, Z^mZ, produces N distinct orientifold planes, 06g, defined as the 
fixed loci of mZQs^ [g = 0,1, - ■ ■ ,N - 1]. Interpreting the operator identity, Q's/'^nes/'^ = mZQ^ , as a similarity 
transformation by the generator 6~^/^, shows that the 06g-planes trace in Tf a sequence of A'' lines related by the 
half-rotation angles, 0^2^ with 06g = &~'^06o. The orbifold symmetry also imposes conditions on the individual 
I?6-brane stacks and the open string sectors. For the case of pairs of brane stacks /x, v passing by orbifold fixed points, 
the invariance under 6^, ClTZQ^ is ensured by requiring that the CP gauge factors of open string sectors, (/x, v), realize 
projective representations for the gauge embedding unitary matrices, 703 , 7n-R0s , obeying the projection conditions 

Aa = Tes./xAATei,!.; = -lanes, iyXA'yane9,fi^ (11.32) 

holding for fj, = v or fx i^, where one must allow for mode dependent complex phase factors determined by the 
quantum numbers of the modes. The RR tadpole cancellation conditions generally admit the simple solution involving 
traceless twist gauge embedding matrices for the orbifold group elements, Trace(7e,^) = 0. 

For the case of branes intersecting the orientifold planes at generic angles, 6^, hence not traversing the orbifold 
fixed points, the rotations 0^ act non-trivially on brane stacks, so that the CP factors are not constrained. To 
ensure that the ZJ-brane setup is orbifold group invariant, one must include for each stack of D6^-branes its — 1 
images D6,j,^^ mider the rotations 6^, and similarly for the mirror images, DQ^'^, [5 = 0, 1, • • • , iV — 1]. Thus, the 
intersecting DG^^-brane stack wrapped around the three-cycles [11^] is made symmetric under the 0^ identification of 
by introducing the image _D6^g-brane stack wrapped around the image three-cycles, Iles/i- Each Z)6^-brane stack 
at generic angles is then described by the equivalence class (orbit) of — 1 branes, rotated images of the reference 
brane, Hg = 0^/i, accompanied by the orbit of N rotated mirror image branes, Q^jj,' = flTZQ^n. The non-diagonal 
open string sectors are described by the orbits, (/i, Vg) ^ {v'g, /i'), (/i, v'^) ^ {vg, /i'). Note that the reverse orientation 
pairs, (/Lt, Vg) ~ {vg,ii)\ are related by conjugation, as discussed previously, and that being identical, the modes (/x, Vg) 
and {fikji'g+k) need not be included simultaneously. The Zn + Zj^QTZ group elements are represented on the 2-d 
vector space of wrapping numbers by the matrix transformations 

^e^f""") = ("""A, rnze^ ^e'^f'^A = f^^Y (n.33) 

such that the column vectors of lattice cycles, (g^), transform by the transposed matrices, 0^ and WZQ^ . The 3-cycle 
volume is a function of the cycles equivalence classes given by the product of three one-cycle lengths, 

V{Q,) = n(20, [Lf = {niri my,) • • ("^Z^,)] (11.34) 

with g^ denoting the 2x2 matrix for the metric tensor in the lattice coordinate basis. 

The diagonal orbit for the Na mirror stacks, (-D6[a], -D6[a']) = {D&ag^ D6a' ), generates the gauge symmetry group, 
U{Na), as the diagonal subgroup of the 2N direct product of group factors. The sectors ([a], [b]) = {agi,hg^) for the 
pair of Na, stacks of £>6Q/D66-branes, carrying the gauge group U{Na) x U{Nb), consist of N distinct subsectors, 
{a,bg), labelled by the relative rotations, g = g2gi^- The A'' distinct mirror subsectors, {a',bg), are similarly defined. 
The off-diagonal open string sectors, {[a], [b]) ~ (a, bg) and {[a], [b']) ~ (o, b'g) carry bifundamental representations for 
U{Na) X U{Ni,) whose multiplicities must be combined algebraically. Recalling that opposite signs are assigned to 
the complex conjugate modes of opposite helicities, one can express the net chiral multiplicities of the bifundamental 
modes as 

N-1 N-1 

3=0 3=0 

[lab, = - nimi), U'^ = n -{nirhl + n^^m^)] (11.35) 

where the summations extend over the A^ distinct subsectors in a given equivalence class belonging to the same gauge 
group representations. The non-chiral spectrum may thus be expressed as, {lab + ^ab){Na,N})) + Aab{Na, Nj,) and 
{lab' + ^ab'){Na, ^b) + ^ab'i^a, ^b), with the model dependent integer numbers of vector pairs denoted by, Aab and 

Aab'. 

The diagonal sectors, ([a], [a]) and ([a], [a']), include the A^ distinct subsectors, {a,ag), carrying the adjoint repre- 
sentation Adja of U{Na), and the A^ distinct subsectors, {a,a'g), carrying the antisymmetric and symmetric tensor 
representations, Aa, Sa, of U{Na). The net chiral multiplicities for these modes are given by 

Adja : = ^ E^-«' s„] : = E ± ttV)' ^"-^6) 

g g L LI ag 
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where the multiplicity for the adjoint (real) representation modes is halved in order not to double count the equivalent 
charge conjugate pairs, (a, Og) and (a, a^v-g)- Since the rotation angles in the supersymmetric type orbifolds 
are given by SU{3) unitary matrices, the adjoint representation modes, ([a], [a]), are localized at branes with angles 
obeying the supersymmetric cycle conditions, J2i ^i,b ~ oYlii = mod(l). Hence, they form chiral supermultiplets 
of A/" = 1 supersymmetry. By contrast, for branes intersecting at generic angles, none of the non-diagonal sector modes 
form chiral supermultiplets. We also note that the extra vector mode pairs, Aab, in bifundamental representations are 
expected to decouple through the tree level Yukawa couplings, fatflh^a, flbf«-bKb, involving the singlet components 
of the adjoint representation scalar modes, Ka and Kb of U{Na) and U{Nb). 



III. TREE LEVEL STRING AMPLITUDES FOR BARYON NUMBER VIOLATING PROCESSES 



We discuss in the present section the string amplitudes for the baryon number violating tree level processes taking 
place in the gauge unified models with intersecting branes. For the familiar two-body nucleon decay channels into 
meson-lepton pairs, the dominant contributions arise from the subproccsses involving four matter fermion fields 
and, in supersymmetry models, from the subprocesses involving two pairs of matter fermions and sfermions, which 
must be subsequently dressed by one-loop gaugino or higgsino exchange mechanisms. The low energy limit of these 
amplitudes is represented by baryon number violating local operators of dimension 6 and 5, obeying the selection 
rules, AB = AL = — 1. Other exotic nucleon decay channels can also arise from higher order subprocesses involving 
either four fermions interacting with a gauge or scalar boson or six fermions [20] . These are represented by dimension 
7 and 9 operators obeying the selection rules, AB = —AL = — 1 and AB = — 2, AL = 0. We shall present here a 
detailed treatment for the former processes but only a qualitative treatment for the latter. 



A. Four fermion processes 

L General structure of amplitude 

The string amplitude for processes involving two distinct conjugate pairs of incoming massless fermions. A' = 
■^'{hflfafJ)^ obtained from the general formula in Eq. (11.24) in the simplified form 

A' = {T,+T,)X^l\x < Vll^^^^mi/^^^^ixW^^^ >, 

[Ti = Tr(AiA2A3A4)° % = Tr{XiX3X2Xi)] (III.l) 

involving only the reduced amplitude A(1234) and its reverse orientation counterpart, since the other two permutations 
refer to forbidden target space embeddings. We find it convenient in the present work to introduce the primed 
amplitude, A', obtained by removing the space-time momentum conservation factor. A' = A/[i{2Tr)^5^^\ki + k2 + 
kz + ki)]. Making use of the useful formulas 

< . . . >= [x{l - x)]-^/\ < e*'=-^(-0 . . . e^k,.x(.,) >^ ^-.(^ _ ^yt^ 
«^ • • •<*) < So^.ixi) ■ ■ ■ So^Jxi) >= x-^l - x)-Huiru3){u2^^u^), 

< s_e.{x,) . . . s^,,s{x,) >= x-(^^-i)'(l - x)-(^^-i)(^''-i), 
Z^{x) = Zl^{x)Zli{x) =< u_ei{x{) ■ ■■cT^f)>,{xi) >, 

[Zl^{x) = C,x-'^''^^-''\l - x)-'^''^^-''')-'^'''^^-'')ip{x)] (III.2) 

one finds that the dependence on intersection angles in the power exponents of x and {I — x) cancels out upon 
combining the various correlator factors. The coordinate twist field correlator consists of quantum and classical 
partition function factors, Z^^{x) and Z^i(x), which are evaluated in Appendix A. The combined contributions from 
the trace over CP factors, the Wick contractions of superconformal ghost fields and of spinor twist fields and the 
coordinate twist field correlator lead to the following final formula for the string amplitude 



A' I dx 







5irix-^-i(l-a;)~*-i/-^(x)^e-^'cJ'(^) 

cl 



-52T2ar-*-i(l - x)-'-'l'-i {x) J2 e" 



j(2) 
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= C(Ti + T2)5i / dxx-'-\l - a;)-*-i [| (2sin(7r6l^)7^ ^ (x) ^ e-^"' 
•^0 7 ^ cl 

[Si = (mi7^U2)(m37p"4), ^2 = (wi7''w4)(w37m'"2), 

^(^) =n^^(^)' ^Ka;) = ''"^^^'^ (i32Gig2 + i3iG2gi), 

7 

C" = C]j2sin(7r6'^), C = 2wgs] 



(III.3) 



where the coefficients, Bi,B2, and the functions, Gi{x),G2{x),Hi{x),H2{x) and Ii{x), are defined by the formulas 
quoted in Eq. (A. 16), while the classical action factors, S^j\x), S^^\x), arc defined in Eqs. (A. 22) of Appendix A, 
and will be discussed in more detail in the next subsection. The equality of the direct and reverse orientation reduced 
amplitudes, which is used in obtaining the second form of the amplitude with the dependence on CP factors factored 
out, follows in a non-trivial way from the selection rules on the branes intersection points and angles, the relation 
between the Dirac spinor matrix elements. Si = —S2, differing by the permutation 2 <-^- 4, and the transformation 
properties of the function //(x) and of the classical action under the a; — > 1 — a; change of integration variable, 
Ij{x) — > 7/(1 — x), sl^\x) ^ 2;). We shall keep track in intermediate results of the direct and reverse 

orientation reduced amplitudes, A(1234) + ^4.(4321), despite the fact that these arc equal in orientifold models. 

The overall normalization factor C = 2ngs in Eq. (III. 3) is determined from the factorization of the low energy 
amplitude with respect to the massless gauge boson s- and ^-channel exchange poles. These are asssociated with the 
contributions from the regions a; — » and a:: — » 1, as will be discussed in detail below. For a shortcut derivation at this 
stage, we consider the large radius limit, r ^ 00, where the classical action factor reduces to unity, e^"^":' 1, and the 
a; — > limit of the amplitude reproduces the s-channel massless gauge boson exchange pole of the 7-d gauge theory 
on the D6 world brane with gauge coupling constant, = g%Q = 2TTgsa' {2TT\/a')^ . Matching the leading massless 
pole term in the string amplitude A'^f to the massless gauge boson exchange term in the field theory amplitude, A'j^^, 
suitably transformed by applying the closure formula for the trace over four CP factors, 



A',t^C{Ti+T2)SnTi [ dxx-'-\ln{l/x))-i 
Jo 

= A'f, = gUTi+T2)Si J (03^ 

2 

yields the previously quoted result [19], C = ^^^^^^ = 27: g sa' . 



(III.4) 



2. Classical action factor 

_q Area 

The tree level contributions to the classical partition function from zero mocks , Zci{x) = er ''' ~ '^^ci , 
are represented by the sum over the embedding of the disk on the lattice of four-polygons weighted by the expo- 
nential of the classical action which identifies with the polygons area. These correspond formally to the holomorphic 
and antiholomorphic instanton embeddings of the disk on . The classical action term in the reduced amplitude 
^(1234) contains a factor for each given by a double series sum over the lattice of the large 2-d tori with cycles 
given by the Z)6A/756B-brane segments, (i^, L^^, images of the disc intervals, (a;i,a;2), (a;2,a;3), as illustrated by 
Figure 8 in Appendix A. 

Using the complex number notation for the plane coordinates, we denote the equations for the branes A, B as, 
Xa = ^aLai Xb = S,bLb, Ha; is G -R]; the position of intersection points as, /i ~ Ar(a;i), /2 = X{x2), /s = X{x^); 
the straight line distances between them as, = f2 ~ fi^ ^23 = fs ~ f2': ^^'^ the winding numbers around the large 
tori as, va, vb- For notational simplicity, we suppress here the index / of the Tf complex planes. While superfluous, 
the suffices A, B on and (5^ are retained for the sake of mnemonics. Upon circling the cycles Ci, C2 surromiding 
the segments (a;i,a;2), (a;2,a;3) the classical part of the coordinate fields along the branes A, B transform by elements 
in the lattice {La, Lb) shifted by the distance of intersection points. The extended boundary conditions are expressed 
by the monodromies 

V2Ac,X = 2TTVA = 27r(l - e''^'''^){6t2+ PaLa), [5^2 = fi " f^, La = haBi + mACa, Pa e Z] 
V2Ac,X = 2nvB = 2tt{1 - e'''''%{5g + pbLb), [Sg = fi - fi, Lb = ubBi + msCs, Pb & Z] 

(III.5) 
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using the definition in Eq. (A. 17), where the integers pA, Pb label the winding numbers of classical solutions and the 
factors 2TT represent the Tf tori periodicities. The factors depending explicitly on the open string sector intersection 
angles, Oda = —0, reflect our use of closed contours, Ci = C'l — C'[ and C2 — C2 — C2, composed of the mirror contours 
around {xi,X2) + {x2,xi) and {x2,X3) + {x3,X2), while noting that the coordinates along the lower and upper paths 
are related by the complex rotation of angle, 2Tr6nA = — Stt^. 

A convenient parameterization for the D-brane equation is obtained by introducing the longitudinal (L) and trans- 
verse {t) vector directions with respect to the wrapped cycle, (n, m) , 



L = Li + iL2 = nei + me2 = (m + r'n)e2 = (— + im), 
£1 

62 



m = m + nT[, t' = — = t[ + ir^ = -U, \L\ = (m^ + nVa^)^ = {L*L)i 



t = t^+it^ = -^(m + mr^), [t- L = = ^{t*L) = 0, P = ^{t*t) = 1] (III.6) 

where we continue using the case of up-tilted torus, while setting the length scale to unity, 62 = i, for simplicity. 

The c-numbers L and t may also be represented geometrically by the 2-d orthogonal vectors L and t with Cartesian 
components given by the real and imaginary parts of L and t. The segments joining intersection points, (5^, 
decompose into the longitudinal components, efg^^' ^23^b, and transverse components, rf^j = Vri^A, = vEz^b- 
The equations for the branes A, B are then represented as, = Ca-^a: Ka E R, a = A, B] those of the wrapped 
cycles as, Xa = PaLa, while the segments joining the brane intersection points are decomposed in longitudinal and 
transverse directions as, Sa = f-aLa + da = ^aLa + flatai [^a, "Ha S ii] so that the grand tori cycles have squared 
length given by, (27r)^[(pc( + + In the present notations, the monodromy conditions in Eq.(III.5) can 

be rewritten as, \f2h.cX = 2-kv = 27r(l — e^"'')[(p + e)L + 77^], [p G Z; e,?? G i?] and the classical action may be 
expressed by the quadratic form 

Sci = Vrx\vA\^ + F22IVBI' + ^^{Vxivav'b) 

= Vi, [{pA + 42f\LA? + 4'] + Vi^[{pB + el)\LB\'' + d^i] 

+ 2^[Vi2{{VA + 4)iA + + 4)^5 + 4)*], (ni.7) 

where Vij, [V-j = VijY\j{2sm-K9^), (ij) = [11,22,12]] are defined in Eqs. (A. 23) in terms of the Hypergeometric 
functions of the integration variable x with parameters depending on the angles 9, 6'. Useful limiting formulas for 
these functions are quoted in Eqs. (A. 24) and (A. 26). The extra factors 2sin(7r(?^) in the primed coefficients V^j, arise 
as a result of expressing the global monodromies in terms of closed rather than open contours. The low energy field 
theory limit of the classical partition function factor is determined by the end point regions of the x-integral. These 
contribute an infinite series of pole terms associated to the exchange of massive string excitations, with infinite singular 
terms with respect to the kinematical variables s, t, u, occurring whenever massless poles are exchanged in the relevant 
channels. Upon including the classical action contributions, the classical partition function must be transformed by 
the Poisson resummation formula prior to the analytic continuation of the .r-integral, in order to ensure the instanton 
series convergence. Considering, for instance, the unequal angle case, ^ 9', where the massless gauge boson pole 
occurs only in the s-channel via the x~^ singularity, then the use of the limiting formulas in Eqs. (A. 24) and (A. 26) 
shows that V22 is the only function in the classical action, Eq. (III. 7), which vanishes in the limit a; ^ 0, hence 
indicating the need to perform the Poisson formula resummation on the pb series. The resulting modified low energy 
representation of the classical partition function factor in the case, dA 7^ 0, rfs = 0, reads as 



Scl = ^-lipA+eAf+dl]\LA\Wli-{pB+eBf\LB\^V^2-2{pA+€A){PB+eBmLAL*BVl2) 



cl PA,PB&Z 



.{pA+eAr\LA\\Vl,-^-^) 



^C-KiMaP( ^ ^ f,-Vi^\LB\^maBA))\ v^',l4l%-WB3*(«BA)g- 

Pb 'Pa 

[aBA = eB + {pa + ca) , K,- = Vij JJ 2 sin(7r^0] • (ULS) 

The analogous Poisson resummation at the end point, x 1, performed on the sum over pA, is obtained from the 
above formula by substituting, A <-> Vn V22- 

Wc consider at this point the brane stack parallel splitting process which realizes the unified gauge symmetry 
breaking. Since the resulting massive gauge boson arises from the open strings stretched between distant pairs of 
brane substacks, its mass Mx is related to the minimal interbrane transverse distance, d^, by the familiar term 
in the string mass squared spectrum, cx d\. For our configuration DABC of branes, the relationship can be 
derived from the s-channel mass spectrum of the open string sector [B, D) by examining the contribution to the 
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classical action in Eq. (III. 7) from the term Vu. Identifying the leading term in the limit .t ^ of the classical factor 
a;-integrand as, e" 



-Sci ^ ^ leads to the result for the (-B, D) sector gauge boson squared mass 



lim e-^=' =e-l''-l'^i'i =ne" 



-\d'^\^ sin^{-we')ln{5i/x) 



(III.9) 



where Sj = 5{9^ ,6 ^) is the auxiliary angle dependent parameter defined by Eq.(A.24). The consistent implemen- 
tation of the broken gauge symmetry by the brane displacement requires including by hand in the string amplitude 
the extra normalization factor, n/(^^)''™ ■> needed to cancel the prefactor of Recall that the over- 

all constant normalization was previously determined by matching the gauge coupling constant in the large radius 
limit after factoring the classical partition function out of the a;-integral. With the same normalization prescription 
based on the identification of the gauge coupling constant, the partition function must then include the extra factor. 



Special configuration with same pairs of conjugate intersection angles 



We now specialize the results of the previous subsection to the simpler case involving two fermion pairs of equal 
angles, 9^ = 9 \ which is realized by the brane configuration with D = B, C = A. The derivation is straightforward 
provided that due care is taken in dealing with the limit, 9^ ^ 9 ^ . The combined contributions to the four fermion 
string amplitude from quantum and classical terms yields the formula 
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(III.IO) 



where we have used the result, I{x) — > 2F{x)F{l — x), and the abbreviated notation for the Hypergeometric function, 
F{x) = F{9, 1 — 9;l]x). The equality of the direct and reverse orientation terms enclosed inside the brackets is 
established by using the change of integration variables, x ^ (1 — x). The classical contributions consist of two 
multiplicative factors associated to the winding and momentum states with respect to the large 2-d tori generated by 
the lattice vectors. La, Lb- The squared mass spectrum of the {A, A) and {B,B) open string sectors are deduced 
by examining the end point regions x and a; ^ 1 of the x-integral which select the i-channel and s-channel 
poles. Before showing this explicitly, we rewrite the lattice summations over wrapped cycles in a compact form by 
introducing the Jacobi theta function with moduli parameter r, defined by the familiar series representation 



(III.ll) 



The resulting formula for the string amplitude reads 



A' = 2'Kgs / dx 
Jo 
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The duality transformation formula for the theta function, i?[q](t) = (— ir) 1). accomplishes the same task 

as the Poisson resummation formula. At x ^ 0, the theta function factors with argument, ta(x) ioo, are safe, while 
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those with argument, tb{x) — ^ iO, are unsafe, hence requiring the use of a duahty transformation to avoid the singular 
behavior from the factor, F~'^/'^(l — x), as needed to interpret the field theory limit in terms of an infinite series of 
s-channel poles. For eg 7^ 0, the same argument with ta{x) and tb{x) interchanged leads to a series of t-ehannel 
poles. The following two representations of the string amplitude, obtained by applying the duality transformations 
on tb{x) and ta{x), achieve the a;-integral convergence at small x and small 1 — x, respectively, 

nl 
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Substituting now the a;-integrand in Eq. (III. 10) 
expansion of the amplitude, 



](tb) . 

by its leading term in the limit x 
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gives the low energy 
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where we have refrained from writing the suffix / on pa,b, La,b and TA,Bix). It is interesting to note that if the 
reverse orientation term were evaluated after performing the change of integration variable, x ^ {1 — x), one would 
obtain the equivalent representation involving only the s-channel poles. 



A' ~ 27rg, 
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The existence of distinct representations of string amplitudes as dual infinite series of poles is a familiar consequence 
of the world sheet duality symmetry. This makes the comparison with the field theory limit appear rather subtle, as 
will be discussed in the next section. The squared mass spectrum of states in the s-channel, M\ g, include momentum 
and winding modes from the {B, B) and [A, A) sectors. The pole positions reproduce the squared mass spectrum 
for the momentum modes along the N directions longitudinal to brane _B, and winding modes along the D directions 
transverse to brane A. The spectrum in the t-channel ^, is similar with A, B interchanged. The structure of the 
compactification mass spectrum is formally equivalent to that for rotated branes parallel along some torus, with 
the roles of the torus cycles e(, 62 played here by the two brane sides, L^, Lb- The momentum modes are associated 
with the cycle La and the winding modes with the transversally projected distance between the l)ranos A and B. The 

above string squared mass spectrum conforms with the familiar formula, = ^ ^ ^ for open 

strings stretched between parallel Dl-branes wrapped around the torus generated by cycles of length L and shape 
angle a. We recommend Ref. [106] for further discussion of the string mass spectrum in intersecting brane models. 



B. Processes with two pairs of fermion and scalar superpartner modes 

In the supersymmetric unified models, alongside with the B, L number violating contributions of D term type to 
four fermion subprocesses exchanging colored gauge bosons, F term type contributions can occur from tree level sub- 
processes exchanging colored higgsino modes between two pairs of massless matter fermions and sfermions, A{tp^p(j)(f)). 
The dominant chiral F term operators of dimension 5 are of form, [QQQL]p and [U'^D'^U'''E''']f- We study here the 
string theory predictions for supersymmetric models by focusing on the configurations with two conjugate pairs of 
localized open string modes with intersection angles 9 and 6'. The tree level contributions to the D and F term 
operators can be identified by considering in turn the two transition amplitudes on the disc surface, 

Av{4i^24>3U'd = 1 < v[y\x,)vt^^\x2)v[-e}^{x3)V^\xi) >, 

J VCKG 
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As{4h'<l''H'^) = j < v[-'^\x,)vt'\x2)V^_y\x^)V^\xi) >, (III.16) 

where we have signaUed the vector and scalar character of the two couphngs by the sufBx labels, V, S. The correlator 
involves same inputs for the vertex operators as those introduced in the preceding subsection. 

Let us start with the amplitude Ay- The massless fermion modes and the low- lying scalar modes arc assigned 

the S'O(IO) weight vectors, r^(l) = (+ + +,-+), r^(2) = ( , ), r^(3) = (111,00), r^{4) = (-1 - 1 - 

1,00)+, (000, 10), implying the squared masses, Af| = M| = 1 - ^{\9[\ + lO'^l + |6'(,|). Note that r^{A) refers to 
the weight vector shifted by picture changing. Since the massless fermions have opposite space-time helicities, the 
only allowed Dirac spinor matrix element is the Lorentz vectorial coupling, ui7^fc4/iU2. This structure can also be 
inferred by making use of the operator identity, u"u2 Sa{xi)Sa{x2)'4'iJ.{x4) = Cr^U2 — > ■^ui'ynU2. Since the 

only relevant term in Eq. (11.23) for the picture changed vertex operator Vg!^\x4) is that involving the spinor term 
{k ■ ip)4'^, the same correlator factor, Z{x), with the ground state twist operators only, appears as in the four fermion 
amplitude. The resulting chirality diagonal string amplitude given by 

/•I / 

X / dxx-''-\l-x)-'-^+^^i'^^'-^'h-^x)e-^c]\ (III.17) 
Jo 

has indeed the Dirac spinor structure expected from a gauge boson and gaugino exchange amplitude. For massless 
scalars, setting, i J2i — \ ^ ^ = !> consistently with the assumed supersymmetry of the model, reduces A'y to 
an amplitude of similar form to the four fermion amplitude. 

We discuss next the amplitude, ^5, in the low energy limit, ki 0. One expects this to contribute to F term 
operators since the incoming massless fermions have the same helicity. In the picture changed operator, V0?\4), we 
need only retain the terms associated with the spinor fields along the internal space directions, tlJ'^{z). The following 
choice for the 50(10) spinor and vector weights of the states is then practically forced on us, 

r^(l) = r^(3) = (+ + +, T±), r^(2) = ( -1,0,0 ; 0,0), 

r^(4) = (-1, -1,-1; 0, 0) + ( 1,0,0 ; 0, 0), (III.18) 

where the constant shift in the spinor weight, ^^"^(4), arises from the picture changing, tp'^{z)dX'^{z)a0'{x4) ~ (2 — 

Xi)~^ Tgi J {xi)e'^^''^^-\ and the imdcrlincs refer to the possible permutations of the entries in correspondence with the 
choice of the index J. We make henceforth the definite choice, J = 3, which corresponds to using the vector weights 
for scalar modes, r^(2) = (00 - 1,00), f^(4) = (-1 - 1 - 1,00) + (001,00) = (-1 - 10,00), with squared masses 
determined by the condition of unit total conformal weights, 

= h{Ve{2)) - 1 = l(|^i| + _ |^3|) _ 

= h{Ve'{A)) -1 = 1- i(|(?'i| + \6'^\ + |^'3|) - Ml. (III.19) 

The coordinate twist field correlator includes two factors, Zi{xi) =< a_Qi{\)crgi{2)a_gii{i)aQii{'^) >, [/ = 1,2] 
identical to the previously studied correlator factor, and one new factor involving a single excited state twist field, 
Zj=z{xi) =< cT-e ' {^)'^eJ (2)CT_g' j (3)fg'j (4) >. For the present, we introduce the following shorthand notation for the 
latter four point correlator with a single excited twist field, 

Zj{x) =< a_gj{xi)agj{x2)(J_g'jix3)fg>j{x4) > 

= fj{xi) < a^gJ {xi)agj {x2)(T_g' J {X3)ag' J {X4) > . (III. 20) 

Upon performing the familiar gauge fixing choice of insertion points, xi = 0, X2 = x, X3 = 1, X4 = X ^ 00, 
extracting out the appropriate factors of X by writing, Zj{x)/Zj{x) fj{x), and using the known correlator factors 

Ez + _ 0'(e''-i) + (0''-i)^ (jjj^21) 

one finds the string amplitude in the form 



-J-/-, 

< e 2 e ^^""'e 2 >^ x 2 (1 - 

(< SSSS > \l) < SSSt > |/=3 = 



7=1,2 



A's = iCs{Mki)u2{k2)){Ti+T2) 



22 



/' 
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-ar)-*-^+^[(^^+^=-^^'-(^^+^^-^3)l/3(a;) J] (^/7^(x)e-^«'). (III.22) 



The expHcit dependence on the intersection angles in the exponent of (1 — x) cancels out in the case of interest 
involving supersymmetric cycles. 

To reach a concrete final result, we modify here our initial choice of the initial states so as to deal with 
the more tractable correlator involving a conjugate pair of excited coordinate twist operators rather than a sin- 
gle one as in the above case. This choice is motivated by the fact that the correlator of interest, Zj{xi) =< 
(y-0j{xi)Tgj{x2)(T_g'.T{x3)Tg'j{x4) >, Can bc morc readily accessed within the formalism set up in Appendix A. 
In fact, the latter correlator would arise if one assigned to the scalar mode (/>2 in -4s of Eq. (III. 16) the choice of vertex 
operator 

[r„^(2) = (00 - 1,00), M| = \{\e^\ + l^^l + \e^\)] (III.23) 

where we have displayed the weight vector of the mode ^2 and the formula of its squared mass. For the remaining 
modes, 03^- <i>'i, in -^s, we continue using the same inputs as above with the choice of complex plane, J = 3. 
The excited state twist field correlator, Zj{xi), can be evaluated in terms of the ground state correlator, Zj[xi), by 
considering the representation for the ratio of these functions 

J, _ Z'j{xi) ^ < a-ejixi)Tgj{x2)a_e'j{x3)fg'j{x4) > 

Zj{Xi) ~ < (7_gj{xi)a0j{X2)cr_g'j{x3)crg'j{X4) > 

= -1 hm {z-X2)^-^'{w-X4/'gj{z,w), (III.24) 

where gj{z, w) is the same correlator denoted as g{z, w) in Appendix A. We next perform the same choice of SL{2, R) 
gauge group fixing of the Xi variables and extract out the appropriate X factors from the two partition functions, 
Zj ~ X^^ Zj ~ X^^^~^ \ A simple calculation yields the following result for the ratio of correlation functions 

f^ix) = C'x'-'\l - xy-'''d,lnls{x), (III.25) 

where the calculable normalization factor C" will be left unspecified. Combining this with the familiar results for the 
correlator factors, 

i 

< ssst >\j= a;-(^'-i)(^'-i)(l - a;)-(^'-i)(^''-i), (III.26) 
yields the final formula for the low energy string amplitude 

X dxx-"(l-x)-*+^S/(^'-^'''^(ln/3(x))]J(/7*(x)e-^'''(")). (III.27) 

This result exhibits the chiral structure expected from the exchange of a fermion mode. At this point, we remark 
that since Ml = O(TOg) is finite, the process at hand is energetically forbidden so that the interest of the present 
calculation is academic at best. One can suppress M| by choosing vanishing small angles, 9^ , but at the price of 
dealing with delocalized modes "01 and (/)2. 

Before closing this discussion, we briefly indicate one possible route to evaluate the partition function Zj in 
Eq. (III. 20). For this one may use the limiting representation, 

Zj{xi) = lim {X5 - x^f' < dX-^{x5)a_sj{l)asj{2)a_g,j{3)as'j{4) >, (III.28) 

and evaluate the resulting five point correlator by a similar method to that used by Frolich et al. , [34] for the open string 
modes with mixed ND boundary conditions. Alternatively, one could consider the bilocal correlators, g{z, w), h{z, w), 
obtained by inserting the quadratic products, dzX^ d-ujX'' , d^X"^ d^X^ in the correlator for Zj{xi), and apply the 
energy source approach reviewed in Appendix A by writing the general representation on the functions consistent 
with the constraints. We leave to a later work the feasible task of implementing these calculations. 
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C. Higher order processes with baryon and lepton number non conservation 

In spite of the stronger suppression of baryon number violating processes initiated by the dangerous operators of 
dimension V > 7, the study of these contributions in grand unified theories is motivated by the need to test variant 
gauge unification schemes involving lower mass scales and different selection rules on B, L non conservation [20, 21]. 
For orientation, we provide a brief overview of the baryon number violating processes from higher dimension operators 
in Appendix C. In the present subsection, we present a qualitative discussion of the string amplitudes in intersecting 
brane models associated with the V = 7 local operators coupling three quarks with single lepton and Higgs boson 
and the V = 9 local operators coupling six quarks. Our presentation here will remain at a general level without 
commitment to any specific model. 



1. Five point amplitudes 

With hindsight from the general structure ofD = 7 operators, we consider the tree level five point amplitude involv- 
ing the localized open string modes of four matter Dirac fermions and a single scalar boson, ^5 = ^(V'lV'IV'sV'l'/'s)) 
defined by 

-^5 = E /n < ^-7^(^i)^r'^(^2)K!;3^\x3)<-^\a;4)C(x5) > . (III.29) 

perms i ^(^^G 

For an acceptable embedding of the disk onto the internal torus, one must require that the modes intersection 
angles obey the condition, —^3 + ^4 + ^5 = 0. In the low energy limit of interest, ki — > 0, only the terms in the picture 
changed operator, Vq^\x^), depending on the internal space directions, contribute. The following unique choice of 
weight vectors for the massless modes, consistent with iJA-momentum conservation and the GSO projection (odd 
number of — | entries), must be assigned to the vertex operators in Eqs. (III. 29), 



-fil) = -r-i?,) = ( ,±T), r-(2) = r-(4) = ( ,±±), 

(5) = (- 100, 00) + (100, 00) = (000, 00) . (III.30) 



The string amplitude can be expressed in abbreviated form as 
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VcKG 
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Ti = Tr(Ai • • • As), T2 = Tr(A5 • • • Ai), =< s_gi{xi)sei{x2)s_gi(x3)sgi{x4)( *'''^^''''^''\] >, 

UJ 3 4 \SgI=2.3{X5)J 



T ( ^6'=^{x^^ \ 

^,f^=<a_ei{xx)uei{x2)o_Qi{x3)aQi[xM = >]. (III.31) 

5 

With the Mobius group gauge fixing choice, [xi, • • • ,0:5] = [0,x, y, \,X 00], the correlator factors for the ghost, 
spin and spinor twist field operators are evaluated by means of the familiar rules 

< e~^^ • • • e~^^^ >= [xy{y - x){l - x)(l - y)]~^^^, 

(u"' • ■■ul") < Saii^l) ■ ■■Sa4,{x4) >= 2/^^(1 - x)^^ {ui{ki)Yu2{k2)){u3{k3)jfiUi(k4)), 

= ^l= a;-(«'-5)\(«'-5)(e3-i)(y _ a;)-(^'-5)(^l-i)(i _ a;)(^'-5)(«i-5)(i _ 

(III.32) 

We shall not attempt here an exact evaluation of the five point coordinate twist field correlators, ^^(jy because of 
the significant labor involved in this task. Nevertheless, the n point correlators of the coordinate twist fields are 
expected to have the general structure [24, 26], < n"=i '^Oii^i) >= C<Tlli:^j{^i ~ 2j)~ = (i-^')(i-^j)Det~ 2 (VF), where 
W denotes the period matrix, whose (n — 2) x (n — 2) entries give the period integrals of the (n — 2) independent 
holomorphic dififerential forms over the (n — 2) independent cycles circling the pairs of insertion points in the cut 
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complex plane. This result motivates us in introducing the following definitions, obtained by including only the 
explicit pair contraction factors with the dependence on angles determined by the conformal symmetry 

x(i - y)-m(^-ol)+ela-ei)) ( -^iC^' y) ) ^ (ni.33) 

V-^2,3(a;,y)/ 

where !Fi{x,y), ^2.3{x,y) include the determinants of the period matrices for the correlators, $f and Combining 
now the various contributions to the five point string correlator, we deduce the final form of the amplitude 

A'^ = Cz{Tl+T2){ul{kl)-^^U2{k2)){m{k:^)^^Ui{ki)) 

X dx dyx~'^{y - xy^il- y)~^PiJ^2:F3e~^''' + perms, (III.34) 
Jo Jx 

where we have explicitly displayed the reduced amplitude associated to the reference cyclic ordering of the insertion 
points, Xi < ■■■ < xc,. Upon associating the five modes in V-'iV-'IV-'s'^l'/'s to the sequence of open string sectors 
{D, A){A, D){D, B)(B,C){C, D), one sees that only the reduced amplitude with the trivial permutation (12345) is 
non-vanishing if the branes D, A, B, C are all distinct. In the case B = A, the reduced amplitude with the permutation 
(14532) would be added, and in the case B = A, C = D, the reduced amplitudes with the permutations (14532) and 
(12534) would be added. 



2. Six point amplitudes 

We turn next to the string amplitude for the six quark subprocesses initiating the AB = —2,AL = processes 
of N — N oscillation and two nucleon disintegration. The discussion developed by Kostelecky et al., [102, 103] is 
followed to some extent, since we only focus here on the low energy limit of six fermion tree level string amplitudes 
involving three conjugate pairs of localized fermion modes. The assignment of CP gauge factors and flavor and color 
quantun numbers will remain implicit, without reference to any specific model. The relevant six point amplitude, 
Af, = A{ijjlip2'^jliJ-'4:'4'l'>p6) , admits the following representation in terms of the correlator with vertex operators inserted 
on the disk boundary 

t6 



■^6 = E / < K^7\a.r)yJ-^^(a=2)ld7\x3)Fjr^\a.4)l/i;i^(a;5)Fj+^\a='') >, 

(III.35) 



perms 



where the massless fermion modes are assigned the unique choice of spinor weights, ( jTT) or (+ + +,=F±)- 

Without proceeding any further, it is easy to convince oneself that the /fA-momentum conservation forces the low 
energy limit of this amplitude to vanish. Nevertheless, whether the AB = — 2 processes are forbidden or not at tree 
level cannot be concluded until one has examined the possibility that VEV induced contributions might arise from 
higher order processes. Hindsight from the unified field theory models with left-right gauge symmetry [107], suggests 
considering the higher order baryon number violating operator, u'^ud'^dd'^dA^ , coupling six quarks with a scalar mode 
carrying the SU (3)c x SU{2) l x SU{2)fi gauge group quantum numbers of dileptons. A; ~ (1, 1, 3). The corresponding 
V = 10 operator would yield a finite contribution to the V — 9 operators once the electric charge neutral component 
of A' = T^A^ acquires a finite VEV. The six quarks in the local operators occur in three pairs, two pairs or a single 
pair of electroweak group singlets while the scalar boson multiplet A' must have at least one component which is a 
singlet under the Standard Model group. Guided by these observations, we consider the seven point string amplitude, 
Ar = ^(V'lV'IV'a^l^sV'e^T')) with localized open string modes inserted on the disk boundary, 

Ar= f < v:!7^(^i)^^^^(^2)vl7'(x3)vjr^\x4)yi;/^(x5)v:^;^^(x6)<°^(x7) >,(iii.36) 

subject to the condition on the intersection angles, ^7 = 6c, + 0^. while adhering to our convention that all the angles 
6i are positive. The assignment of spinor weights, consistent with the if^-momentum conservation and the expected 
structure of the Dirac spinor matrix element, allows for the unique choice 

r^(l) = (+ + +, T±), r^(2) = ( , ±±), r^(3) = (+ + +, ±t), r^(4) = ( , TT), 
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r^(5) = (+ + +, ±T), r^(6) = (+ + +, T±) + (00 - 1, 00) = (+ + -, T±), 

r^(7) = (-1 - 1 - 1, 00) + (001, 00) = (-1 - 10, 00), (III.37) 

where the scalar mode has squared mass, = 1 — 5(1^7] + l^rl + l^rl). The single surviving reduced string amplitude 
can be expressed in abbreviated form as 



A'j = CriTi + 7^)[ay(ui7^U2)(u37M"4)(u5W6) + as{uiU3){u2U4){u5Ue)] 

i_gI^{xQ)tgl{x7) 
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r, , , , 1 ^(^1) vi'^d) , , f t_gi{xe)tei{x7)\ 



X <a_0i{xi)--- { « \ > +perms, 

\(T_gl {X6)agI {X7)J 

[Ti = Trace{\i ■ ■ ■ A7),' T2 = TraceiXr ■ ■ ■ Ai)] (III.38) 

where the upper and lower entries for the twist field correlators refer to / = 3 and / = 1,2, and the Lorentz covariant 
structure is described by the calculable functions 05, ay which we shall not attempt to determine here. The Mobius 
group gauge fixing of the insertion variables is set as, [xi, 2:2, X3, X4, X5, xg, X7] = [0, w, u, y, x,l, X ^ 00]. In order to 
obtain finite contributions to At, the brane models must comply with restrictive constraints. Among these are the 
conditions that the scalar boson be a localized mode from a non-diagonal open string sector and that the Z?6-brane 
embedding of the world sheet boundary forms a 7-polygon which closes up to finite transverse separations. The 
vector like gauge unification schemes are out of the game because of the automatic B + L or B — L conservation 
present in these models. While the left-right symmetric gauge models are natural candidates, as illustrated by the 
brief review in Appendix C, we find that none of the intersecting brane models discussed in the literature [54, 57] 
can be extended so as to include finite contributions to the AB = — 2 seven point amplitude at tree level. Indeed, 
for the D-brane models with Pati-Salam gauge symmetry, SU{4)c x SU{2)]^ x SU{2)ij, no pcrturbativc open string 
mode can occur with the required gauge group representation, (10, 1,3), which includes the boson A'. The higher 
dimensional multiplets which arise in the realization presented in Ref. [57] carry rather the representation, (10, 1, 1). 
For the minimal left-right symmetric model with gauge group, SU{'i)c x SU{2)l x SU{2)fj, a localized dilepton scalar 
mode A' carrying the triplet representation under SU(2)ii — SU{2)a can possibly arise in orbifold models from the 
non-diagonal open string sectors of type (a, 9^a). 

We close at this point the discussion of string amplitudes for the exotic baryon number violating processes. Our 
admittedly unfinished presentation here is due to the non-trivial structure of the higher order correlators for the 
coordinate twist operators. We have attempted here to clarify the conditions leading to finite string amplitudes for 
five, six and seven localized matter modes of similar structure to the four fermion amplitudes. For the V = 9 baryon 
number violating operators, only the VEV induced contributions can possibly be present in the intersecting brane 
models. As in the effective field theory model context (cf. Appendix C), finite amplitudes may arise from dimension 
V = 10 operators involving massive scalar multiplets in higher dimensional representations of the unified gauge group 
having Standard Model singlet scalar components which acquire a VEV at some larger mass scales. 



IV. FIELD THEORY LIMIT OF STRING AMPLITUDES 



A. Effective low energy field theory 

Before going into the applications, it is important to identify the relevant fundamental parameters and determine 
their relation to the parameters of the world brane field theory. The underlying string theory is characterized by the 
Regge slope mass scale parameter, a' = l/m^, and the string coupling constant which is linked to the dilaton field 

VEV, Qs = e<*>. Alongside with the complex structure moduh, = + and the wrapped three-cycles winding 

'"2 

numbers, (n^,m^), we introduce here the two real Kahler moduli parameters, R and r = 1/Mc, representing in an 
average sense the length scale radii for the 6-d torus and the three-cycles wrapped by the D6-branes. The volumes of 
the internal manifold and three-cycles, Vx and Vq^, are described in terms of these parameters as 

Vx = {2nfR^ Vq^ ^ Ui^nLi) = r']l{2nCi) = {2nrfC,, 

1=1 I 

[£. ^l[ci,Ci = ^ = + ||)^ rhi = mi- <t//], (IV.l) 
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where Qa denotes the three-cycle wrapped by the D6a-brane. To establish contact with the low energy physics, we 
must also introduce an infrared cutoff mass parameter, Mx, whose meaning will be discussed at length shortly. It is 
useful to keep in sight the familiar formulas expressing the 10-d gauge and gravitational coupling constants of type II 
supergravity theory, g\(j and Kio, and the tension, RR charge and gauge coupling constant parameters of Dp-branes, 
Tp, jip and QDp, as a function of the string coupling constants, 

1^ o 7/2 '2 9ym 2kio 

Kio = — = Stt'/ a , gw = — — = —f=, 
9s gi Va' 

^ _ _ 1 _ \/7r 



^ 9s gs\/a' {2tt \f^)P K.(2^:^f^y-^' 
9-A,. - \9-nl = i(2.V^)Vp = (IV.2) 

The 4-d gravitational Newton coupling constant. Gat, and the gauge interactions coupling constants on the Dpo-branes 
world volume theories, ga, are given by dimensional reduction as 



167r _ Vo-kVx _ SmfVx 

^ = "^^ = 2^ " ~ WNI' 



where the enhanced gauge symmetry case for Dp-branes overlapping the Op-planes is obtained by replacing, g\ g\l^- 
For reference, we also display below the inverse relations expressing the fundamental string parameters as a function 
of the effective low energy parameters for the case of £>6-branes with p = 6, 




(27r)3 83/2(27r)3Ar' 
( Mpgl 



(IV.4) 



The dependence on the gauge group factors included in the parameters, ka = ^j"; with Vq denoting the volume of some 
reference three-cycle, allows one to write the string gauge unification relations as, kaUa — ax- The gauge interactions 
coupling constant, gx, may be traded with either gg or rUg, by using the relation, imsr)Ylj C^^ = gl^^ax^^^. For 
fixed ax, increasing m,,r, causes gs to increase. The ratio between the manifold to three-cycle radii, A = Vx/Vq = 
{R/r)^, controls the relative strength of the gauge and gravitational interactions. We note that the duality between 
compactifications of M-theory on G2 holonomy manifolds and of type II string theory on Calabi-Yau manifolds with 
intersecting _D-brancs suggests [17] that the values of A arc restricted by an upper bound of order unity, while the 
larger values of A are favored [108] in order to reduce the mass scale of the string theory axions, F^. However, we do 
not gain much in the present work in including the gravitational interactions input, since this trades Gn for the bulk 
radius parameter R which has no direct impact on the open string observables. One cannot simultaneously match 
Gat and ga to the observed values without imposing a wide disparity between the bulk and brane radius parameters, 
r/R « 1, hence requiring finely tuned shape moduli parameters, rl/r^ << 1. This circumstance rules out the TeV 
scale models with fS-branes, which require large extra dimensions in the bulk. Of course, the restriction to toroidal 
models with r R need not be in force in models with lower dimensional D5- or Z)4-brane backgrounds, where the 
existence of subtori unwrapped by D-brancs gives the ability to arbitrarily raise the bulk volume Vx relative to Vq. 
For the £)5-brane and Z)4-brane models, the elimination of Wg in favor of Mp would yield the formulas for the string 
coupling constant 



a^VQml alMpV^ g^VQnis aT M^i^V'J^ 

(2^)' 4V2ttV^ 277 V2(2^)i/3y_^/6 

which are seen to exhibit a variant dependence of gg on the gauge coupling constant, aa- 

Let us now discuss the mass cutoff which is inevitably needed to regularize the infrared divergences in string 

amplitudes. In unified gauge theories, this is naturally associated with the gauge symmetry breaking mass scale 

parameter, Mx- Since the wrapped cycle volume Vq^ fulfills the same role as Mx, it is convenient to introduce 

1 

the dimensionless parameter [19], L{Qa) = Vq^Mx- The compactification mass parameter, = ^ = , [C = 
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ri/("-i^''i^ + '^^^^2^)^] is linked to the unified symmetry breaking mass scales A4x through the string threshold 
corrections to the gauge coupling constants, Aq. The moduli dependent functions, Aq, are commonly defined by 
writing the one-loop renormalization group scale evolution for the running gauge coupling constants as 



Aa = ^277^7 =2 ba log( — ), fca = "f^ (IV.6) 

9a{Q) 9x Vq 



where ba denote the one-loop slope parameters of the gauge group factors, Ga, determined by the massless charged 
modes, and denotes the string theory effective unification mass, determined by matching onto the low energy field 
theory. We choose here to treat Mx as a free parameter without imposing any condition on its relative ordering with 
respect to the compactification mass, Mc, except for the request that both masses are boimded by the string mass 
scale, {Mx, M^) < nis- Loosely speaking, the condition Mx < M^ implies that the unified symmetry occurs in 4-d, 
while the condition Mx > M^ implies that that this occurs in 7-d. Prom the interval of variation assigned to Mc by 
the combined consistency conditions, < 1, mg/Mc > 1, yielding, (a^l^CDs < ^ < 1, we see that the possibility 
Mx > Mc is not excluded. 

To describe the ratio of the string theory amplitude to that of the equivalent field theory, we follow the procedure 
of Klebanov and Witten [19]. Ignoring momentarily the regularization of divergences, we write the string amplitude 
as, Ast = 27r5s/(6')TiS/(2m^), where the extra factor ^ accounts for the orientifold projection, the a;-integral I{9) 
is a calculable function of the dimensionless free parameters, m^r = ms/Mc, s = mg/Mx, and T, S contain the 
dependence on the gauge and Dirac spinor wave functions. The equivalent field theory amplitude is now assumed to 
include the same factors. Aft = 2-KaxTS /M\. In correspondence with the different prescriptions of identifying the 
effective parameters, we deduce the following three predictions for the ratio of string to field theory amplitudes. 



_^ Ast _ 1^9 si ml 
^^'f =Ait- l.axIMl'^^^ 

gl"'aa''^V^'^Ml MlmsVQ aaVg^MpMl 



m, (IV.7) 



where the successive entries inside brackets are obtained by eliminating gs alone, nis alone, or gs and nis together, 
in favor of the low energy gauge and gravitational coupling constants. It is instructive to compare these results with 
the ratio of M-theory to field theory amplitude [17] 



Am = Cm^^^J^^J^ =^ UM/f = ^ ^ ^CmL^/HQ)c^xK [Mn' = 2K?i(27r)«]. (IV.8) 



Getting an estimate for the coefficient Cm was the main motivation of Rcf. [19]. The appearance in the string theory 

2 

amplitude at fixed gs of the same power dependence on the 4-d gauge coupling constant, A oc a^, joins with the fact 
that in the large radius limit of the string and M-thcory models the gauge unification occurs in 7-d. A useful discussion 
of this issue is presented by Burikham [109]. That this property is specific to the D6-branc models appears clearly 
from the different dependence on the 4-d coupling constants of in Eq.(IV.5) for the £)5-brane and £)4-brane models. 
The Z)6-brane case is known to be unique in admitting a purely geometrical lifting to the M-theory models [75]. 



B. Regularization of string amplitudes 

The enhancement of string amplitudes reflected by the non-analytic dependence on the 4-d gauge coupling constant 
takes place naturally in the large compactification radius limit, r = l/M^ oo. The effect is maximized by considering 
fermions localized at vanishingly small distances relative to the compactification radius [19]. It is important to realize 
that what justifies considering directly the ratio of amplitudes in Eq.(IV.7) is the assumption that the string amplitude 
is self-regularized. However, that the large compactification radius limit of the string amplitude is only finite in the 
case of I^G-branes can be easily seen by rewriting the a;-integral truncated to the leading contribution at a; ^ as 

TT* / da;a;-'*-^(ln(l/a;))-* = / d^q^^, (IV.9) 
Jo J Q ^ s 

where the right hand side exhibits the sum over propagators for the plane wave modes along the L = 3 effectively fiat 
directions. The convergence at g — > would not hold for L < 3. (Of course, the factor Inl/x slows the convergence 

of the integral near the endpoint, x = 0, but the problem this would pose is only of technical order.) In the D5- and 
£)4-brane models with the x-integral at L = 2 or 1, the singularities at a; ^ are seen to be not integrable ones. 
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Owing to the slowly convergent .T-integrals, it is not clear whether the enhancement survives finite compactification 
radii. Indeed, as just noted, the integrable singularity at x = is specific to the case of _D6-branes where the infrared 
singularity at (f ^ is cancelled by the integration measure (fiq. At finite compactification radius, the massless gauge 
boson modes in the s— and channels are separated by a finite mass gap from the massive string modes, and the 
presence of infrared divergences from massless pole terms is then unavoidable. The pole terms from exchange of string 
modes can be separated out by means of the familiar analytic continuation method. While the non-trivial functional 
dependence of the a;-integrand from the coordinate twist field correlator precludes using a fully analytical procedure, 
it is still possible to consider a semi-numerical regularization of the x-integral where one removes by hand the massless 
s- and t-channel gauge boson poles originating from the contributions to the x-integral of form, dxx~''~^!Fo{x) and 

dx{l — x)~^~^J^i{x), where J^o{x), J^i{x) represent the limiting forms of the remaining factors in the integral at 
a; ^ and x ^ 1. 

We now discuss two distinct infrared regularizations that will be used in the sequel to obtain numerical predictions in 
the finite compactification radius case. In the first prescription, the regularized amplitude is constructed by subtracting 
out by hand the massless gauge boson pole terms in the relevant channels and substituting these by the corresponding 
massive pole contributions, 1/s l/{s — M'^). We thus write the string amplitude as, Ast = T^9sOi'[SiTiIi + S2T2I2], 
and substitute for the direct and reverse orientation integrals, Ji, I2, the subtraction regularized integrals with the 
massless poles removed, 



hie) ^ ir^ie) = h- /i,o - [/i,o = 

h{e) ^ 7r^(^) = 72 - 72,0 - 72,1, [72,0 



\Lb\{-s) 



72,1 = 



'a' 



\LA\{-ty 
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where the indices and 1 in 7j^o, 7j_i, [i = 1, 2] refer to the limits x — > and x ^ 1. In the case, = es = 0, where 
massless poles are present in both s- and f-channels, the regularized string amplitude is defined schematically as. 



A=[A- A\:o^o - ^Ix^i] +AI + A. 



\LB\{s-Miy 



Ca 



\LA\{t-Miy 



A\ 



Cb 

\Lb\ Jo 



[ dxx ^ \ A\x^i = -^ f rfx(l-x) 

Jo \^a\ Jo 



-t-i 
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where the constant coefficients Ca, Cb are determined from the limits at the end points of the x-integral. We now 
assume that the field theory amplitude has the same dependence on the Dirac spinor and gauge wave functions. Using, 
for simplicity, the approximate equalities, |Lyt| = \Lb\ = \L\, 7[^^ = 72^^ = I^'^s, and ax = eta, f^x = Mx, we can 
express the ratio of string to field amplitudes in the fixed gs prescription as 



// 



+S2T2{ir'\L\ 



+ 



-s + 



'a' 



l + -9iax^\L\-^M%r^o{e). 



X 



27raa 



5iTi 
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S2T2 
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In the fixed nXg prescription, the ratio would instead read as, 'R-sif — 1 + \M\m,s\L\V^^{G). For the case with non 
coincident intersection points, = 0, £5 7^ 0, the f-channel pole is absent, so that using the above prescription 
without the terms A\x^\ and A^, yields the ratio, T^-s// = 1 + M\ms\L\V^^ . An analogous result holds in the case, 
EA 7^ 0, es = 0. 

The second, perhaps more natural, regularization procedure is based on the description of the unified gauge sym- 
metry breaking by the deformed configuration obtained by splitting a finite distance apart the unified gauge theory 
7)6-brane stack into two stacks. The transverse displacement parameter is related to the symmetry breaking mass 
scale Mx by the leading contribution at x — !■ in the classical action factor, as displayed by Eq. (III. 9). This amounts 
to moving along the Coulomb branch of the gauge group moduli space of vacua and is T-dual to a Wilson flux line 
around the cycles normal to the brane stack. Denoting the resulting multiplicatively regularized x-integrals as 7['f', 
and using the same simplifying assumptions as above, leads to the following regularized ratio of string to field theory 
amplitudes in the fixed Qs prescription. 



Aft 



reg 



-S2T2il'') 



SiTi 
-s+M'i 



+ 



S2T2 

-t+M'i 
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The ratio in the fixed rus prescription would instead read, TZs/ / — ^M'^ms\L\i^^3 . Recall that the regularized integral, 
7''^^, now includes the brane induced form factor, n/(<^/)*^'" '"^^ ; needed to restore the correctly normalized 
gauge coupling constant. In the case with equal angles, the correction factor simplifies to (6)^^ . The numerical 
comparison of the above discussed regularization prescriptions will be presented in the next section. 

To conclude, we comment on the calculation of the ratio of amplitudes in orbifold models. The vertex operators of 
the {a,b) modes must now be summed over the orbifold group images, V)^a,b) ^ J2g ^{a,bg)- We choose not to include 
the normalization factor at this stage, since the amplitude normalization will be determined in this case by 

comparison with the gauge coupling constant. Other conventions would lead to the same end result. Upon combining 
the summations over the orbits in the various vertex operator factors, the compatibility conditions from the target 
space embedding on the brane boundaries leaves a single orbifold group summation of equal reduced amplitudes, which 
then introduces an overall factor N . One must now recall that promoting the torus to an orbifold entails a reduction 
of the torus volume by the group order factor N. The gauge and gravitational interaction parameters are expressed 

ytor 

via dimensional reduction by same formulas as for tori, with the three-cycle volume identified as, Vq'^ = -jf-- This 
change can be implemented by using either {ga^)orb = J*"'' and {Mp)orb = , with fixed string parameter, 

Qs, or {gs)orb = ^^"jj'"' = ^2(2ot)^jv^ ' '^^^h fixed gauge and gravitational coupling constants. Repeating the calculation 
of the string amplitude factorization on the gauge boson pole term, one arrives at the formula, gs = "(2^p°jv° > with 

the imderstanding that the three-cycle volume is that evaluated for the torus, Vq^^ = {2nr)'^ 11/ (-^D- thus 
summarize the conversion from torus to orbifold descriptions by the following schematic correspondence 

(At)tor = -^^TSIi9) =^ iAst)orb = ^^TSI{9). (IV.14) 

This shows that upon passing from tori to orbifolds, one may retain the same normalization constant, C = 2Trgs, 
while using the modified formula for the string coupling constant, gg = mlax\L\/N. 



C. Large compactification radius limit 

To prepare the ground for the calculation of four fermion amplitudes in semi-realistic models, we present the results 
obtained in the large compactification radius limit using the self-regularized x-integral. The ratio of string to field 
theory depends then on the branes intersection angles only. Lacking a simple analytic approximation for the a;-integral, 
we perform the quadrature numerically. The relevant integral 1(9) in the equal angles case, 9^ ~ 9 ^ , with the classical 
action factors omitted, is evaluated numerically for the toy orientifold model [19] realizing the SU (5) group unification 
with the mirror pair of I?6/Z)6'-brane stacks, allowing for a single massless 10 matter multiplct. One advantage of 
this simple model is in studying the anticipated enhancement effect in the amplitude, 10 • 10^ • 10 • 10^, as a function 
of the brane-orientifold angle only. The amplitude in Eq. (III. 10) can be expressed in terms of the x-integral 

where we used the identity, [x(l — x)]~^ = x~^ -|- (1 — x)]~^, and the symmetry of the integrand under a; ^ (1 — x). 
We pursue here the study in Ref [19] by evaluating 1(0) for generic three-cycles, not obeying the supersymmetric 
restriction on the angles, J^i = 0. As already noted, the infrared finitcness of the integral does not dispense 
us from taking care of the sensitive integrations at the end points of the integration interval. Reasonably accurate 
numerical values for this integral can be obtained by the simple procedure which consists in subtracting the leading 
term at a; ^ of the integrand and adding back its contribution to the integral, as illustrated in the following formula 

I{9) = [ dxT{x) = [ dx[J^{x) - J^oix)] + [ dxTo{x) = [I{e) - Io{9)] + Io{9), 
Jo Jo Jo 

[^(-) = m^l^^/^'' ^o(x) = lim ^(.), 

(IV. 16) 
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The same subtraction procedure can be used for the x = I end point, with the limiting function .^"0(2;) substituted by 
J^i{x) = lima;^i J^{x) = JFo(l — x). The subtracted integral displayed in the second entry of the above equation has a 
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simple representation in the two special cases involving three equal angles, 9^ = 9^ ^ 9^ = 9, or two unequal angles, 
9 = 9^ = 9"^ ^ 9^ , respectively. The resulting analytic formulas for these two cases are given by 



h{9,9,9) = 277^/2 (In (5)- 
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Io{0,0,9s) 



r3/2 



|ln(^3/^)l^ 
+e//(-A)[7r - 2arctan 



Oh (A) In 



1-(1 



InS \ h 
In 53 ' 



1 + 



-1 + 



In (5 
IniSs 



(IV. 17) 



In ^ positive and negative 



where 6jj denotes the Heaviside function (Qjj (y) = 1 for y > and for y < 0) with A ; 

in correspondence with the cases, S < (^3 and S > S3. 

The numerical results for the integral I{9) as a function of the single variable intersection angle are displayed in 
Figure 1. We note that the integral rapidly vanishes at the end points, corresponding to parallel or anti-parallel brane 
stacks, and is maximized at the midpoints, 9^ = ^, as expected from the symmetry under 9^ ^ {1 ~ 9^), at fixed 
/. Away from the end points, the integral vary slowly and monotonically inside the interval, I{9) ~ (6 — 12) for 
9 e [0, i]. Recall that for supersymmetric cycles satisfying 9^ = 2, the integral was found to vary inside the range, 
I {9) ~ (7 — 11). No significant differences thus arise in going from the case of supersymmetric cycles [19] to that of 
non-supersymmetric cycles. 



V. BARYON NUMBER VIOLATION IN GAUGED UNIFIED ORBIFOLD MODELS 



A. Z3 orbifold 5*17(5) unified model 



Wc here employ the formalism developed in Section III to examine the enhancement of string amplitudes for 
nucleon decay processes at finite compactification radius. We consider two type // string theory realizations of SU (5) 
gauge unification using £>6-branes. The first uses a non-supersymmetric Z3 orbifold model [66] and the second a 
supersymmetric Z2 x Z2 orbifold model [68]. To complement our discussion of the solutions for the non-chiral 
spectrum of the first model, we summarize in Appendix B the properties of models using the Z3 orbifold with £'6- 
branes, based on the work of Blumenhagen et al.,[66]. 



1. Low lying mass spectrum 

The SU{5) gauge unified Z3 orbifold model of Blumenhagen et al.,[66] uses the minimal setup of two D6-brane 
stacks Na = 5, Nc = l with effective wrapping numbers, (Y^, Za) = (3, i), (Yc, Zc) = (3, — ^), solving the RR tadpole 
cancellation condition, N^Z^^, = 2. This realizes the gauge group SU{5) x U(l)a x U{l)c with three chiral matter 
fermion generations, Fi, fi, z/f, localized at the brane intersections and several copies of non-chiral fermions. Hi, Hi 
and Ka^i, K^^i- The gauge group representation content and multiplicity of these modes is displayed in the table 
below, in correspondence with the open string sectors to which they belong. 



Mode 


F 




H + H 


Ka 


Brane 


{a', a) 


(a,c) (c,c') 


(a',c) 


{a, a) 


Irrep 


3(10,1)2,0 


3(5,l)_i,i 3(1,1)0,-2 


n(5, 1)1,1 +n(5,l)_i. 


_i 7„^* (24, 1)0,0 



A vertical bar has been inserted to separate the chiral from the vector representations. The complete spectrum 
for matter fermions is of the form, (3 + Af)F + ApF, (3 + Aj)/ + Aj/, (3 + A^.y + A^.V, where the model 
dependent integers Ap, Aj, A,yc count the multiplicities of mirror vector generations. Of the two Abelian gauge 
group factors, the linear combination with charge, Qan = ^Qa + Qo is anomalous while the unbroken orthogonal 
combination, Qfree = %" ~ is anomaly free. When the unified gauge group breaks down to the Standard Model 
group, the leftover unbroken Cartan subalgebra generator of SU{5) combines with the anomaly free gauge charge 
so as to yield the unbroken gauge symmetry, U{1)b-l- The anomaly free gauge symmetry U{l)free, which assigns 
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charges Qfree = —3, —2) to the representations (10, 5, 5), precisely identifies with the accidental symmetry of the 
minimal SU{5) unification related to the Y and B — L charges as, Qfree + |^ — Qb-l = 0. 

The scalar modes content of the model includes the ingredients needed to accomplish the Higgs mechanism breaking 
of the various gauge symmetries. The scalar field VEVs for the Ka^i singlet components of the adjoint representation 
modes of SU{5), are needed to break the unified SU{5)a gauge symmetry at the unification mass scale Mx; those 
for the modes are needed to break the U{1)b-l gauge symmetry at some low intermediate mass scale; and those 
for the modes {Hi, Hj) are needed to break the electroweak SU{2)l gauge symmetry at the Fermi mass scale. The 
presence of tachyon or low-lying scalar modes Ka^i with mass considerably lower than the string mass parameter, is 
facilitated by the fact that these modes arise in massless chiral supermultiplets. The U{5) singlet components of the 
Ka^i scalars have the ability to produce large Dirac masses for the various vector pairs of spin-half fermions. The tree 
level string couplings of the matter fermions, represented by the disk configuration in the amplitude I.a of Figure 2 
and the similar couplings for the higgsino and gaugino like fermions, Hi + Hi, Ka,i, produce the Yukawa trilinear 
interactions, 

LeFF = hfj^iFiF]Ka,l + h{./fi~f]Ka,l + hfj^iHiHjKa,l + h^„^iKa,mKa,nKa,l, (V.l) 

which are necessary to decouple the extra fermion modes. For this to occur, it is necessary that the fermion mass 
matrices, hfjj < Kaj >, h{j ^ < Ka^i >, hf^ i < Ka^i > and h^^ i < Ka,i >, have rank not smaller than the vector 
representation multiplicities, Ap, Ajr, n and I^'^^- 

The intermediate mass scale breaking of U (1) b-l can be accomplished by the localized scalars of the (c, c') sector, 
which enter with multiplicity 3. For a tachyon mode to be present, one must cancel the negative contribution to the 
string squared mass depending on the angles by the positive contribution depending on the branes c, c' distance. For 
this, one must require that the brane c does not pass through the origin of the coordinate system. While the Dirac and 
Majorana mass matrices of the left chirality neutrinos vanish, finite Majorana masses can occur for the right chirality 
neutrinos via the quartic coupling, i)'^ i)'^ v'^ v'^ , once the v'^ tachyon scalar raises a VEV. The corresponding correlator, 

< v[y\v'')V^°\v''^)v[y\v''^)V^^^\v'') >, represented by I.d in Fi gure 2, is expected from the structure of the 

^ -c ^2 

string amplitude in Eq. (III. 22) to give a contribution to the Majorana neutrino mass of form, mt,c . 

Although the orbifold brane setup of interest does not preserve any supercharges, one may still use the freedom 
in the moduli parameters to dynamically select a vacuum having low-lying or tachyon scalars with the appropriate 
gauge quantum numbers. For the locahzed modes in non-diagonal sectors, (a, 6), this possibility can be explored 
by considering the stability tetrahedron with edges traced in the space of the brane pair intersection angles, 9^^, by 
the four equations expressing the vanishing of the scalar modes minimal squared mass in Eq.(II.14). The preserved 
supercharges on the tetrahedron faces, edges and vertices generate the supersymmctries, J\f = 1,2,4. The faces 
separate the inside domain where all four squared masses are positive from the outside region where at least one 
squared mass is negative. Clearly, one should avoid scalar tachyons with the gauge quantum numbers of matter 
fermions, F, f, but would welcome those with Standard Model singlet components whose condensation can accomplish 
the Higgs mechanism symmetry breaking [52]. 

The unified SU (5)a gauge symmetry breaking down to the Standard Model corresponds to the deformation described 
by the brane splitting, a a. + b. Following Blumenhagen et al.,[66], we consider the setup consisting of the three 
stacks, Na = 3, Nb = 2, Nc = 1, with the effective wrapping numbers, (Ya, Za) — (Yb, Z^) = (3, ^), solving the RR 
tadpole cancellation condition. While the Coulomb branch deformation is restricted to mutually parallel D-branes 
in the orbifold equivalence classes, [a] and [h], one may also consider the deformation leading to non-parallel a and b 
brane stacks. Removing from the gauge symmetry group, U (3)a x U (2)b x U (l)c, the anomalous Abelian factor, Qan = 
3Qa + 2Qb + Qc, leaves the extended Standard Model model gauge group, SU(3) x SU(2) x U{1)y x U{1)b-l, [Y = 
— 5<3a + ^Qb, Qb-l = — g(Qa — SQfc -|- 3Qc)]- The localized fermions consist of three chiral generations of quarks, 
leptons and antineutrinos, {qi, u^, dl, U, e^, t'f), and several copies of adjoint representation modes, Ka^i, K^^i, 
with model dependent multiplicities, as displayed in Table I. We have not included in the table the vector pairs 
of fermion modes in the sectors (6',c) -I- (a',c) with gauge quantum numbers. Hi ^ [(l,2)o,i,i -|- (3, l)o,i,i], Hi ^ 
[(1, 2)o,-i,-i -I- (3, l)o,-i,-i], = 1, • • • tIt]. In spite of the vanishing chiral multiphcity, /^^ = 0, one expects that 
the massless spin half modes and the low-lying complex scalar modes in these sectors enter with a finite non-chiral 
multiplicity of the form, h'cg^ = —h'cg^ = n, [g\ ^ 52] characterized by the non-vanishing integers, n. The presence 
of tachyons among the low- lying electroweak doublet complex scalars in the modes, Hi+Hi, is needed to accomplish 
the electroweak gauge symmetry breaking, as indicated already. To relate the electroweak Higgs mechanism to the 
brane recombination process, b' + c ^ e, one would need a fine compensation in the tachyon mode mass squared 
between the negative contribution from the angle dependent terms, a'M"^ = min[i '^ji^lOl/J), 1 — ^ X^/ ^^'^ 
the positive contribution from the interbrane distance. In spite of the analogy with the models incorporating an 
M = 2 supersymmetry subsector [52] or with the supersymmetry-quiver models (first article in Ref. [53]), involving 
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D6-branes wrapping parallel one-cycles, the present model features specific differences, unlike the claim made by 
Axenides et al., [84]. In particular, the conjugate modes, carrying the representation n{Nh,Nc) + n{Ni,,Nc), are 
neither A/" = 2 hypermultiplets, since the cycles in the three planes are not parallel, nor Af = 1 chiral supermultiplets, 
since none of the wrapped cycles are supersymmetric. 

The candidate electroweak Higgs boson modes. Hi Hi, account only in part for the quark and lepton mass generation. 
The graphs I.b and I.c in Figure 2 represent the trilinear Yukawa couplings of down quarks and leptons with the 
electroweak Higgs boson, Qi,i aD'g^^^{H) c^v and Li,c{H)cb'gE^, ^. One sees that the trilinear Yukawa couplings exist only 

for the down-quarks and charged leptons, qd'^H and le'^H , while the corresponding ones for the up-quarks, qu'^H, are 
forbidden. Indeed, the scalar mode with the needed quantum numbers, H ^ (l,2)_3 _i,o, cannot be realized as a 
perturbative open string mode. The obstruction to generate masses for both the up and down quarks is a general 
feature of the minimal type models which realize the right handed triplet quarks as the antisymmetric representation 
of SU{3)c- The same problem also occurs in models using the Z2 x Z2 orbifold [68] and the x Z^ orbifold [85]. In 
the flipped SU{5) model, it is rather the down-quarks and leptons Yukawa couplings that vanish. One way to account 
for the top quark mass here is by appealing to the dynamical gauge symmetry breaking or to composite Higgs boson 
mechanisms. For completeness, wc note that to remedy the shortcoming of intersecting brane models in favoring a 
separable structure of the quark and lepton mass matrices with rank one, recent proposals consider using a hypercolor 
type scheme [110], loop corrections [111], or higher dimensional operators [112]. 

2. Unified gauge model and Standard Model realizations 

To each set of effective wrapping numbers {Y^, Z^) of a given lattice type, there arise a large number of realizations 
described by the wrapping numbers {nj^,iTijj) solving Eqs. (B.9). Each solution comes with a specific non-chiral 
spectrum and D-hva,ne intersection angles, constituting input data indispensable for evaluating string amplitudes. As 
found by Blumenhagen et al., [66] for the model with three brane stacks, {njj_,mjj, [/x = a,b,c] each pair 
admits some 36 solutions, hence yielding a total of 4 x 36^ realizations for the four inequivalent lattices. However, 
only a few solutions gave acceptable mass spectra for the scalar modes. Some of the natural constraints to impose 
on the physical solutions include the absence of tachyon scalar modes with the quark and lepton quantum numbers, 
and the presence of at least a single scalar tachyon among the electroweak singlet bosons u'^ ~ (c, c'), charged under 
U{1)b-l, and the electroweak doublet bosons H ~ {b' ,c). 

The scan over solutions for the wrapping numbers of the extended Standard Model with three brane stacks, is 
conveniently performed by means of a numerical computer program. We study the two cases involving parallel and 
non-parallel £)6a/£'6(,-brane stacks. We begin with the simpler case of parallel branes a and c, whose discussion 
overlaps that of the gauge unified models. After elimimating the solutions with tachyon scalars carrying quark and 
lepton quantum numbers, we could not find any solution satisfying the requirement that tachyon scalars be present for 
both i^^ and H, H modes. However, since several solutions exist with tachyon scalars in either the i>f or Hi, Hi modes, 
we consider the minimal retreat from our initial goal allowing for the solutions with at least one tachyon among the 
modes i^f. Restricting the numerical search to the interval of integer wrapping numbers, [—6, +6], we find 16 and 10 
inequivalent solutions for the T^/Z^ lattices AAA and BBB, but no solutions for the lattices AAB and ABB. The 
wrapping numbers for the solutions of the Z)6a/-C)6c-brane setup with parallel a, b stacks are quoted in Table II. The 
displayed information on the intersection numbers for members of orbifold orbits, la'ag , lacg , led ■ la'cg , [5 = 0,1,2] 
fully determines the non-chiral spectrum. There is no need here in quoting the intersection angles since these are 
easily determined from the displayed information. (Note that = —0^.) In the various realizations of the unified 
model, the number of mirror matter generations take values in the range, Aj? ~ [0 — 5], ~ [0 — 3], Ai,c ~ [0 — 5], 
the vector higgsino modes. Hi, Hi, [i = 1, • • • n] appear with the finite multiplicities, n = 2, 6, 7, and the gaugino like 
modes, K^^u = 1; " " ■ > I^'^^i x = a, b], enter with the uniform value of the multiplicity, J^'^j = I^"^^ = 7. Note that 
none of the solutions in Table II satisfies the property that the brane pairs, b', c, are parallel along a single complex 
plane. In the cases where tachyon modes for Hi + Hi appear, the squared mass of these modes sets at the typically 
large value, a'M'^{H) ~ —0.11, thus justifying our previous claim that no J\f = 2 subsectors arise in the present 
model. 

We next examine the Standard Model realizations with non-parallel stacks, (n^,m^) ^ (n^,m„). Requiring the 
absence of scalar tachyons with quark and lepton quantum numbers, yields some 50 solutions for lattice AAA, 40 
solutions for lattice BBB, and no solutions for the lattices AAB, ABB. The maximally constrained selection of 
solutions (Case I), requiring a nearly massless or tachyon scalar for the modes v'^ and H + H, retains 4 solutions for 
lattice BBB. The less constrained selection of solutions (Case II), requiring a nearly massless or tachyon scalar for 
the modes i/^, retains 5 solutions for lattice AAA and 10 solutions for lattice BBB. We display in Table III the 
solutions found for the wrapping numbers of the three non-parallel stacks realizing the extended Standard Model in 
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the partially and fully constrained searches of Cases II and I. 

We now discuss the conditions set on the effective string mass scale, m'^, by the low energy data for the SU{3) x 
SU{2) X U{1)y gauge coupling constants, Qa- The hypercharge embedding, Y = — |Qa + ^Qb, implies the relation, 
g-^ = Qy = + 92 I which can be compared with experiment by substituting the one-loop running coupling 

constants. 



47r 



2sin2 OwiQ) - 1 



= Blog^ + A, 



3a3(Q) «(Q) 
[B = ^63 + 62 - 61, A = ^Aa + A2- Ai] (V.2) 

subject to the assumption that we deal with a non-supersymmetric unified model with no new physics thresholds be- 
tween Q = niz and m^. We can estimate m'g by substituting the experimental values, a{mz) = 127.95, sin^ Owimz) = 
0.23113, as{mz) = 0.1172, [mz = 91.187 GeV] while using suitable inputs for the slope parameters ba of the massless 
modes below m',. (The slopes are calculated from the general formula, ba = \{—Ca{Rs) — ^CaiRp) + llCa(J?v)), with 
Ca denoting Dynkin indices for complex scalar, Weyl fermions and vector modes.) Since all the above listed solutions 
feature equal volumes for the a, &, c cycles, we expect to find, as in the conventional field theory type unification, 
poorly convergent running coupling constants indicative of a low unification mass scale. It is known, however, that 
the presence of massive charged vector modes can significantly improve the compatibility with data. Thus, assuming 
that a new threshold is present at some mass, < modifies the above quoted formula to 
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3a3(<3) ol{Q) 



{B + 6B) \og^=B log ^ + A, (V.3) 
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where the slope parameter correction, SB = -gSbs + 6b2 — Sbi, includes the contributions from the vector modes that 
decouple at niA- We now test this relation by assuming the Standard Model evolution in the interval, Q € [mz,mA\ 
with slope parameters, 

11 4F 22 4F Nh , 20 iV^ 44 

03 = 11 , bo = , oi = — F B = — , (V.4) 

3' 336' 9 6 3' ^ ' 

and including in the interval Q G [niA, m'^] the extra massless vector multiplets, (3 + A)(F + / + v") + A{F + f + 
+ /^*'/ra + if^'^Kb + nH + nS, with slope parameters, 

Sbs = -M^''^ - n, 5b2 = -24^* - n, Sb, = ^ SB = -2{I^^^ + I^'^) + ^. (V.5) 

Assigning the tentative values, = 7^* = n = [0,2,4,7] and ^ = 10~^, we find that the predicted string 

compactification scale sets at the values: m'^ ~ [5. x 10^^, 3.7 x 10^'', 2.7 x 10"^^, 5.3 x 10"^^] GeV. These results 
explicitly demonstrate how a mild hierarchy for the intermediate mass threshold can significantly pull the string mass 
scale down. 



3. Non-supersymmetric flipped SU{5) x U(l)fi model 

We here consider the flipped gauge unified model of the Z3 orbifold constructed by Ellis et al. , [78] (Model I in 
Table II of their paper) constructed with three ZJG-branc stacks Na — 5, Ni, = 1, A^c = 1 with effective wrapping 
numbers, {Ya,Za) = (3, 5), {Yb, Zb) = (0, -5), {Yc,Zc) = (3, -g), which realize the gauge group, SU{5) x U{l)a x 
U{l)b X U{l)c- The chiral spectrum of localized fermions is displayed in the table below in correspondence with the 
group representations and the assigned open string sectors. 



Mode 


F f 


H H S 


Sector 


(a', a) (a, c) + (a, b) (c', c) + (6', c) 


{a',b) (a',c) (6,c) 


Irrep 


3(10)2,0,0 2(5)_i,o,i 2(1)0,0,-2 


l(5)_i,_i,o 1(5)1,0,1 1(1)0,1,-1 




+l(5)_i,i,o +1(1)0,-1,-1 
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On side of the anomalous symmetry group, Qan = 5Qa + Qc, one finds the unbroken gauge symmetry group, 
SU{5) xU{l)fi xU{l)freej [Qfi = 5 (Qa " SQfo — SQc), Q free = Qb] with three chiral matter generations, {F+f+e'^) ~ 
(10 + 5 + 1), along with massless higgsino like modes, (H + H + S) ^ {5 + 5 + 1). The light or tachyon scalars with 
same quantum numbers as the higgsino modes are needed for the Higgs mechanism breaking of U{l)free and the 
electroweak SU{2)l gauge symmetries. The Yukawa coupling. SHH, decouple the higgsino modes once the scalar S 
mode raises a VEV. The absence of low-lying scalars carrying the representations 10, 10 means that one must resort 
to a higher dimensional Higgs mechanism to accomplish the unified gauge symmetry breaking. 

We have developed a computer aided scan of the solutions for the three branes wrapping numbers selecting the 
realizations with scalar tachyons absent for quark and Icpton modes but present for the H, H and S modes. These 
conditions select only the T^/Z^ lattices, ABB, BBB. The maximally constrained selection (Case I) with tachyon 
or nearly massless H,H,S scalar modes gives 6 solutions for ABB and no solutions for BBB. The least constrained 
selection (Case II), with tachyon or nearly massless scalar modes for S only, retains 35 solutions for ABB and 2 
solutions for BBB. We display in Table IV a representative sample of the solutions found for the lattice ABB in 
Case I and the lattice BBB in Case II. 



4- Enhancement effect in four fermion haryon number violating processes 

We now discuss the two-body nucleon decay amplitudes in the orbifold model realization of SU{b). The four 
fermion amplitudes displayed in Figure 2 by the graphs II. a and II. b refer to the configurations 10 • 10^ • 10 • 10^ and 
10 • 5''' • 5 • 10^, which contribute to the proton decay channels, p ^ 7r° + and p ^ 7r° + e^j. The corresponding 
amplitudes in the broken unified gauge symmetry version are displayed by the graphs III. a, Ill.b of Figure 2, with 
graph III.c referring to the neutrino emission decay channel. 

The numerical calculations are performed by setting the unified coupling constant at the value, ax = assuming 
that the wrapped cycles have the same radii, = r, and regarding the gauge unification and compactification mass 
scales, ,s = and nisr = j^. as free; paraincrters. For simplicity, we further set Mgut = Mx. The winding 
numbers and angles of the wrapped three-cycles are determined for each solution. The same numerical value for the 
three-cycles volume, \Ca\ = Yli{Li/r^) = \/7, is found in all the solutions. The low energy constraints allow then to 
express the string coupling constant by the numerical relation, 9., ~ 0.11 {m.srY /N . The weak coupling condition on 
the string coupling constant entails an upper bound on the wrapped three-cycle radius 

.. = ^^ = ^^<1 ^ m.r<(^)..2.Ar., (v.6) 

where N denotes the orbifold group order. Using these inputs in Eq. (IV. 7), yields the following approximate numerical 
formulas for the ratio of string to field theory amplitudes in the fixed Qg and mg cases 

C-^fim.r)H'-^), (V^7) 

ms 10 

where we have set the reference value of the regularized a;-integral at the approximate value obtained in the large 
radius, = 10, as follows from Figure 1, while ignoring momentarily its dependence on the parameters, Mx = mg/s 
and mgr. We note that the ratio has a fast power growth, proportional to {nigr)'^ in the fixed g,, and to {mgr)^ in 
the fixed nig case, if one discounts the suppression effect from the classical action factor in /''^s'. In comparison, the 
analysis of Klebanov and Witten [19] with the reference value of the Wilson line parameter set at, L{Q) = 8, predicts 

the ratio, TZg/f ~ 1.5 (^^)55l(f^)"^/^(^). The fact that this is a factor 4 smaller than the result in Eq. (V.7) 

reflects on the uncertainties in the input parameters. It is also interesting to compare with the ratio predicted by 
making use of the gravitational interactions input, TZg/f — 2.3 10~^( ^*^^ )( '^g^-' )((fM)("~L5^)' since this suggests that 

the compactifications with small wrapped cycles, hence close to the upper limit on = R/r oi order unity, would 
not support an enhancement effect. 

A brief aside on numerical issues is in order. For a good convergence of the series sums in the classical partition 
function, the summation labels should cover the range, \pa\i \pb\ < 6, throughout the interval of radii, m^r e [1,4], 
The larger is mgr, the slower is the convergence of the zero mode summations. The sigma model perturbation theory 
sets the lower bound, m^r > 1, while the string perturbation theory sets the upper bound, Qg « 0.11 {mgr)^/N < 1. 
Nevertheless, we shall consider small excursions at small and large radii for the sake of illustrating the rapid growth 
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of string amplitudes in these forbidden regions. To simplify the implementation of the displaced brane regularization, 
we shall also assume that the transverse displacement takes place in a single complex plane with fixed /, with the 

leading term in the classical action factor of Eq.(III.9) reading as, (Sj)'^^ e^^^'''^ . To illustrate the slow numerical 
convergence of the .T-integral, we show in Figure C a plot of the integrand in the subtraction regularization prescription 
using the definition 

2mx) ^ l.[A'J{CS,{T,+T,))] = [x{l - x)]-\j^^^;^^^)h-'-, (V.8) 

which generalizes the definition in Eq. (IV. 15) for I{9) = dxi{d;x). in the finite radius case. We see that after 
removing the end point singularities due to the massless poles, the integrand still picks up its leading contribution 

from narrow intervals close to the end points. This feature holds true at all the relevant values of m^r. The slow 
convergence of the x-integral is a critical slowing factor in the cases involving distinct pairs of intersections angles. 

Turning now to the predictions for the ratio of string to field theory amplitudes. Ti-s/f; we examine first how the 
results vary for different realizations of the unified Z3 model. Since the intersection angles in the various solutions 
follow repetitive patterns, we need to perform independent calculations only for a few cases. The ratio TZs/ / calculated 
with the parameter values, rUsr = 1, s = nis/Mx = 1, = = (0, 0, 0), is displayed in the following table for the 
representative sample of four solutions covering the complete set in Table II. 



Solution 




TZg/f (Pole Subtraction) 


TZg/f (Brane Displacement) 


AAA 








/ 


(0.879,0.333,0.667) 


3.15 


3.16 


II 


(0.545,0.667,0.667) 


5.93 


2.48 


III 


(0.667,0.667,0.545) 


5.66 


2.64 


V 


(0.333,0.667,0.879) 


5.93 


2.67 


BBB 








/ 


(0.454,0.333,0.333) 


7.54 


3.06 


II 


(0.333,0.667,0.121) 


7.02 


3.14 



We see that small variations in the angles can have a significant impact on the predictions. The sensitivity is 
stronger in the subtraction than in the displacement regularization scheme. The near equality of the results in the 
two prescriptions for solution / must be viewed as a coincidence. Although the angle and volume parameters are 
nearly the same for all solutions, there arc important variations in the wrapping numbers and hence in the fixed 
complex plane volume factors, L^, from one solution to the other. Indeed, we find different predictions for solutions 
with same angles but different wrapping numbers. Prom the formal relation between the a;-integrals in the subtraction 
and displacement regularization prescriptions, /'"'^^ = /'"'^^/(l + ^M'j^ms\L\V^), one infers that I"^"^^ < /""^s. In fact, 
we notice that the predictions for TZ^^ f is a factor 2 — 3 larger in the subtraction regularization scheme than in the 
displacement regularization scheme. Since the displaced brane prescription is more realistic, we conclude that the 
enhancement effect is of typical order of magnitude, TZg/f — 3. (Mx/Mc)^ . 

We now discuss the dependence on mgr based on predictions displayed in Pigures 4 and 5 for the ratio of string to 
field theory amplitudes in equal brane angle configurations. The plots show a rapid rise of the ratio with increasing 
mgr. In both the subtraction and displacement regularizations, starting from 0(1) values near mgr = 1, the ratio 
TZg/f rapidly increases to 0(10) for m^r ~ 2. Note that further increasing m^r would clash with the perturbative 
constraint on g,. The amplitudes ratio rapidly increases with increasing Mx at fixed r. It is larger for Mx > Mc, 
corresponding to the region in the plots, r > s, lying on the right hand side of the vertical line, r = s. The comparison 
between Pigures 4 and 5 shows that the ratio is smaller over the full interval of mgr by a factor of order 2 — 4 in the 
brane displacement regularization relative to the pole subtraction regularization. 

The ratio of string to field theory amplitudes is strongly reduced for intersection points separated by a finite 

longitudinal distance, es 7^ or 7^ 0. In the coupling, Waa'^ot'aPac^ca, the distance cb between the 5 and 10 

modes is generically non- vanishing, while = 0. This is expected from the exponential dependence of the classical 
factor. The results illustrate the level of suppression that may be expected for the flavor non-diagonal nuclcon decay 
processes involving second generation leptons and/or quarks. We see on the predicted ratios displayed in the right 
hand panels (b) of Figures 4 and 5, that a significant suppression occurs even for the moderate value, €b = {jq, jq, jq)- 
Comparing the results in the left and right hand panels (a) and (b) of Figures 4 and Figure 5, we note that the two 
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regularization schemes give a similar dependence on the parameters, with a stronger reduction for increasing TOsT 
taking place in the displaced brane regularization. Following an initial rise, the ratio undergoes a drastic reduction 
for radii, m^r > 2. The fact that near mgr ~ 1, the ratio starts out larger for ee 7^ than for eb — 0, comes about 
because of cancellation effects in the series sum over zero modes due to the nontrivial complex phases introduced by 
using the Poisson resummation formula. The latter factor is a counterpart of the overlap integral between the localized 
modes wave functions encountered in the context of the field theory orbifold models. Since the flavor changing nucleon 
decay processes are drastically reduced except in narrow intervals close to mgr = 1 or for Mx/mg > 1/3, one may 
expect a strong sensitivity to the quark and lepton flavors in these models. 

The ratio of amplitudes in the unequal angle case is plotted as a function of m^r in Figure 6. We see that although 
somewhat weaker, the enhancement effect is still of same order as in the equal angle case. The comparison between 
the left and right hand panels (a) and (6) again shows that a drastic reduction takes place in the non-diagonal 
configuration even with moderately distant intersection points. We note that for distant modes at large values of 
the s parameter, the ratio takes values below unity. That the string theory mechanism studied here might introduce 
a reduction rather than an enhancement joins with the conclusion reached by Acharya and Valandro [113] upon 
evaluating the contribution from the tower of Kaluza-Klein modes in a field theory version of the M-theory amplitude 
compactified on the solvable 7-d lens space manifolds, Q = /Zp. 

B. Supersymmetric SU{5) model of Z2 x Z2 orbifold 

1. Unified model and threshold corrections 

We consider here the supersymmetric gauge unified models with intersecting _D6-branes constructed by Cvetic 
et al., [68] for the Z2 x Z2 orbifold. To avoid repetition, we briefly recall that the point symmetry group consists 
of the two generators, 0, to, and allows the Hodge numbers are, /i(i.i) = 3, = 3 + 48 = 51. The minimal 

SU{5) model is obtained with three D6-brane stacks, Na^ = 10, A^^^ =6, Nf, = 16, with wrapping numbers, 
[(n^, m^)] = [(1, 1) (1,-1) (1, ^)], [{nl, m^)] = [(0, 1) (1, 0) (0, —1)], satisfying the RR tadpole cancellation conditions. 
Each brane stack is made to preserve some J\f = 1 supersymmetry by by tuning the ratio of radii, Xi, so as to satisfy 

the condition on the brane-orientifold intersection angles, J2i ^1 ~ w Cretan "^J^^ = mod 2, [xi = ^]. With the 

_D6-brancs assumed to pass by the orbifold fixed points, one must impose on the CP factors the projection conditions 
involving the gauge twist embedding matrices, 7^,/^, Jui,n alongside with the orientifold projection condition involving 
1(111,11 ■ We shall not further elaborate on the dependence with respect to the gauge factors since this should cancel 
out in the simple-minded definition we have adopted for the ratio of string to field theory amplitudes. The issue of 
the string amplitude normalization has been examined Refs. [73, 114]. 

After the orbifold and orientifold projections, the resulting gauge group is given by SU (5)ai x SU (3)a2 x U Sp{16)b x 
U{l)ai X f/(l)a2, where the symplectic gauge group factor arises from the D6b-hrane being parallel to the 06-plane 
fixed under QTZOu). The massless modes from the diagonal sectors include the adjoint representations, 3(24, 1, l)o,o + 
3(1,8,1)0.0 + 3(1, 1, 136)o,Oj and three massless adjoint representation A/" = 1 chiral supermultiplets, $1, $2, ^3, 
which descend from the initial A/" = 4 gauge multiplet. The non-diagonal sector includes four chiral generations of 
matter chiral supermultiplets, Fi, {fa,i + fa), [^ = 1, • " " 4; a = 1, 2, 3; a = 1, • • • , 16] where the indices a and a label 
the fundamental representations of U{3) and USp{16), along with Higgs boson chiral supermultiplets, C", C", as 
listed in the table below. 



Mode 


F /a 


fa 




Sector 


(ai,ai) (01,4) 


{ai,b) 


(02,4) {a2,b) 


Irrep 


4(10,1,1)2,0 4(5,3,1)1,1 


l(5,l,16)_i,o 


4(1,3,1)0,2 1(1,3,16)0,-1 



The intersection angles for the 10 and 5 modes are parameterized by, F ~ 10 : (^i-^a' ^ ^i' ~ ^01 = 
(— 2a^, 2a^, —2a^), / ~ 5 : 9^^^ = ^6 ~ ^ai — ~ ot^ ■lO? : —\ — ce^), in terms of the angle parameters obeying 
the supersymmetry condition, — + = mod 2, [ira^^^^'^^ = arctan(x/=(i,2))) ^ra^ = arctan^]. We only 
consider here the solution defined by setting ai, a2, at fixed values. 
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The unified SU{5) symmetry breaking to the Standard Model and the breaking of (7(3) and USp{16) are accom- 
plished through the finite VEVs of the adjoint representation scalar modes. The 5, 5 mirror fcrmion generations 
decouple through the mass terms induced by the trilinear superpotential, W = XifaifaC", provided one assumes that 
the tachyon scalars, present among the C^, give the 12 x 16 mass matrix, Aj < C" >, rank 12. The electroweak 
Higgs doublets arise from the massless linear combinations of fai, fa- Since the trilinear superpotential FifjH only 
is allowed by the gauge symmetries, while FiFjH is forbidden, one finds again that the Higgs mechanism generates a 
mass matrix only for down-quarks and leptons. 

We now discuss an indicative prediction for the unified mass scale, Mx , based on a rough calculation of the massive 
string threshold corrections to the gauge coupling constants. These have the additive structure, Aa = A-^^^ + A-^^^, 
in correspondence with the contributions from the string state subsectors of = 2, 1 supersymmetry. From the 
discussion in Ref. [115] for type // supersymmetric intersecting brane models, we borrow the following schematic 
formulas for these two components of the threshold corrections 

A^=2 = Y,biMTi\n{T')\^]+\og{4ne-^-V^)], 

b,I 

= E ba, log (i^) , [y/ = ^] (V.9) 

where 6^^ denote the slope parameters due to the (a, b) modes with branes parallel in a single Tj, -fE = 0.577... denotes 
the Euler constant and the correct use of the formula for A^=^ requires that the interbranc angles lie inside the range, 
\6^\ < 1. A natural definition for the effective; string unification scale, m'^, can be considered by absorbing inside the 
logarithmic term, InQ^/m^, the contributions from A-^^^ which are finite in the large radius limit, thus yielding, 

iTLciy similcirly tra,nsfGr tliG explicit T2 mod-uli dependent contributions to 

to the massless mode logarithmic term by introducing the effective mass scale, = 

which then naturally identifies with the compactification scale, ~ Mc- We now consider the parameterization 
of the string threshold corrections, A^ = — 6aA + kay, assuming that the moduli independent universal constant 
terms in A^, are absorbed into the string unified coupling constant, m^. The redefined mass and coupling constant 
unification parameters can then be expressed as 

Mx A Aa 1 1 y 

-JL = e^ =e-^, —^ — + (V.IO) 
< 9x 9x (47r)2 

To obtain a rough estimate of the effect of A/" = 1 sectors on the ratio Mx/Mc for the unified brane stack 
a, we consider the contribution to A^"^ from a single brane stack b for the choice of interbrane angles, = 
±2(— |, 1^ — A)^ [p = 1, . . . , 10]. With increasing p > 1, the numerical results for A^^^/bab start from the low values 
±5. X 10~* and increase rapidly up to the value, ±1.8, at p = 10. It is important that both positive and negative 
signs of Aa/bab, leading to Mx/Mc ratios smaller and larger than unity, can occur. The favored values lie in the 
range, Mx/Mc 5 — 2, although one may not exclude a pile up effect from several branes which would widen this 
interval. 

It is interesting to compare these results with those found in the M-theory compactification on a G2 holonomy 
manifold, X-j ~ Q x K3, using the lens manifold, Q — S^/Zq, with SU{5) gauge group broken by Wilson lines 
around Q. The predicted threshold corrections [17], A^ = '2,'}2i%TrR.{Q'j/) ~ lOkaT^ + |6o(7o — 7^), involve 
the analytic torsion index, 7^., in the representations of the subgroup U{1)y commuting with 5'[/(5). Defining 
Vq = (27rr)^ = (27r/Mc)^ and using the numerical estimate, L{Q) = e^"~^° = Aq siia^ {birw / q) ~ 10, yields the 

prediction for the mass ratio, = = "^n"''^^ — 



2. Nucleon decay amplitudes 

The numerical results for the ratio of string to field theory amplitudes in the subtraction regularization procedure 

are displayed in Figure 7. The panels (a) and (6) refer to the equal and unequal intersection brane angle cases, 
associated to the nucleon decay amplitudes in the configurations 10 • 10''' • 10 • 10''' and 10 • 5 • 10'^ • 5'^. The dependence 
on the compactification radius is shown for a discrete set of values of Mx at a fixed value of the free complex 
structure moduli, determined by the choice of ai, 0:2 • The results are very similar to those found previously for 
the Z3 orbifold. We thus conclude, in conformity with our previous conclusion from the study of x-integrals, that 
no essential difference arise for the ratio in the supersymmetric models. However, to be more conclusive in the case 
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of the dimension 2? = 5 amplitudes, one should improve the quantitative understanding of the correlators involving 
excited coordinate twist fields. 

We have used so far the information bearing only on the unified model by assuming a schematic representation 
of the gauge symmetry breaking. A fully consistent calculation of the four point baryon number violating processes 
requires specifying the brane setup at low energies. Explicit realizations of the deformed vacua with displaced brane 
stacks have been discussed in Ref. [68]. The brane displacement regularization is explicitly realized by splitting 
some D6-brane stack into separated quark and lepton substacks. One can then adequately suppress the baryon 
number violating couplings of 2? = 6, 5 by increasing the distance between the substacks. We consider here the 
three family extended Standard Model defined in Tables IV and V of Ref. [68] with the setup of six brane stacks, 
8D6a„ 8D6^„ 4D6b„ 2D6b„ 2D6^„ 4D6^„ yielding the gauge group, USp{8)a, x USp{2)^, x U{2)b, x USp{2)b, x 
[/■(4)ci x USp{4:)c^. The brane stack splitting, 8D6c, = 6D6c,^ + 2D6c,i, resulting in C/(4)ci, U{3)c,^ x U{l)c„, 
is used to suppress the baryon number violating processes. With the first gauge factor broken as, USp{8)ai — > 

U(l)s X [/(l)^/. the hypercharge is embedded in the linear combination, Y = — + ^a±5ai_ xhe open 

string sector assignment for the quarks and leptons, q ~ (^i,cig) + {bi,c[^), u'^ ~ (ai)Cig), d'^ ~ {ciq,ai), I ~ 
{c[i,bi) + (ci;,6i), e'^ ^ (ci;,ai), allows to identify the configurations of open string modes which contribute to the 
relevant V = 6 and 5 dangerous operators, 

Oe^ ~ ueu^'^d"^ ~ {ciq,b[){b[, cii){cui. ai){ai, cig), 
• qqql ~ udue ^ (hi, c[q){bi, ciq){bi, c[g){c[i,bi), 

{u^e^u^it) ~ (ai, cig)(cu, ai)(ai, ci,)(ci,, ai). (V.H) 

The above brane embeddings of the disk boundary show that while the non-supersymmetric operators, Oe^ and Oe^, 
and the supersymmetric operator, (w^e^u^d^), are allowed, the supersymmetric operator, qqql, is disallowed. Since 
the one-loop dressing of the operator u'^e'^if^d'^ involves the exchange of electroweak higgsinos, one concludes that the 
nucleon decay from X> = 5 operators should be adequately suppressed in the present model. 



VI. DISCUSSION AND CONCLUSIONS 



We have studied in this work the nucleon decay processes in semi-realistic string models of gauge unification. Our 
aim was to assess the stringy enhancement effect caused in the M-theory and type II string theory models by the 
localization of matter modes in the internal space manifold. The effect has maximal strength in the large radius limit. 
In accounting for the effect, the reduced power dependence of the string amplitude on the unified gauge coupling 
constant, A oc c?x^ , appears as a secondary manifestation in comparison to the strong peaking of the ar-integral at 
the end point regions. The leading contributions from these end point regions arise from the towers of momentum 
modes propagating along the internal space directions wrapped by the branes. 

The localization of open string modes at brane intersection points is formally analogous with that of closed string 
twisted modes at the orbifold fixed points. On practical groimds, an important difference is that while the distance 
between intersection points may be very small relative to the compactification distance scales, the fixed points of 
orbifold groups are generically a finite distance apart. Thus, the inevitable suppression from the classical partition 
function factor in closed string amplitudes may be minimized in brane models provided that the distances between 
intersection point are small compared to the wrapped three-cycle radius. 

The consistent discussion of string unified models requires introducing an infrared matching mass scale in addition 
to the string parameters Qs, tUs and r. This is identified here with the unified gauge symmetry breaking mass scale, 
Mx, which may stand for the open string moduli representing either the Wilson flux line or the unified brane splitting. 
Since we lack a fundamental understanding of how the compactification and D-brane processes stabilize the closed 
and open string moduli, we must regard and Mx as free parameters bounded by and differing by at most a 
few orders of magnitude from each other and from mg- The relative ordering of the parameters Mx and M^ has a 
crucial impact on the enhancement effect. While a definite relation between them is implied by the heavy threshold 
corrections to the gauge coupling constants, we lack a quantitative understanding of these effects to make a learned 
choice on the relative order of Mx and M^- Although our study of the size of Mx/Mc implied by the threshold 
corrections is not fully conclusive, this seems to indicate that the ordering of these scales can go both ways. 

Our numerical predictions confirm that some enhancement of the string amplitude may be present at finite radius. 
The consitcncy requirements select the intervals, Mx/tRs G [1/4 — 1] and m^r G [1 — 4]. We find a significant 
growth of the ratio with TOgr, which is maximal for coincident intersection points. For distant intersection points 
the suppression eflfect from the classical action factor limits the ratio to values of 0(1). If Mx < Mc, the smallest 
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contribution found with lowest Mx leads to ratios of order 2 at coincident intersection points and order 1 for distant 
intersection points. If Mx > Mc, the ratio attains 0(10). The same conclusion applies to the configurations of 
brane pairs with same and distinct intersection angles. Thus, the effect is not restricted to the branc setup with same 
intersection angles, as concluded from the study of Ref. [19], based on a brane setup which did not include the unequal 
angle configuration. 

The application to the supersymmetric model using the Z2 x Z2 orbifold indicates that a similar enhancement 
effect may be present in the nucleon decay processes from V = 5 operators. No essential diflFerence appears in 
supersymmetric models, as also indicated by the similar size of the x-integrals for the supersymmetric and non- 
supcrsymmetric configurations of intersecting angles. Our study has an indirect bearing on the nucleon decay processes 
described by the dimension T> = 7, 9, 10 operators, although we have not carried out explicit calculation of the string 
ampliudes in these cases. We found that restrictive selection rules are set at the tree level by the D-brane embedding. 
No finite contributions to these operators are found in the D-brane models discussed in the present work. 

In conclusion, our calculations confirm the presence of a small enhancement effect which is maximal for diagonal 
configurations with same intersection points, but not necessarily carrying the same group representations. No essential 
changes are observed upon passing to supersymmetric amplitudes associated with the dimension 5 operators. The 
string amplitudes feature a rapid growth in the regimes of large and small compactification radii relative to the string 
length parameter in which the string and world sheet perturbation theories are invalidated. The suppression effect in 
the case of distant intersection points is found to be substantial. The enhancement effect is also potentially present 
for the six fermion or higher order processes contributing to the AB = —2 amplitudes at tree level. However, the 
examination of proposed intersecting brane models [58] realizing Pati-Salam or related gauge unified models indicates 
that this source of baryon number violating is generically suppressed because of the absence of the higher dimensional 
matter modes needed to contribute to VEV induced AB = — 2 operators. 
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APPENDIX A: VACUUM CORRELATOR OF COORDINATE TWIST FIELDS 



We review in the present appendix the calculation of four point correlators for the open string coordinate fields 
obeying the twisted boundary conditions appropriate to intersecting branes wrapped around three-cycles. The re- 
striction to factorisable tori allows us to specialize the discussion to a single Tf torus, parameterized by the 
pair of complexified coordinate fields, X^{z), X^{z), with the relevant part of the stress energy generator given 
by, T^{z) = —-^dzX^ (z)dzX-' (z). For notational convenience, we restrict in the sequel to a single complex plane 
and omit the space component indices / by identifying the string coordinates as, X = X^ , X = X^ . We are 
interested in the correlators involving two distinct pairs of conjugate twist fields, (J±0i , with opposite sign angles, 
Z^{zi) =< (T_gi {zi)(Tgi {z2)(T_g' I {z^)(jgi I {Z4) >. Thc calculatious will be pursued along same lines as in the conformal 
field theory approach developed for the twisted sectors of heterotic string orbifold models by Dixon et al., [25] and 
Bershadsky and Radul [24]. Our presentation closely follows the work of Biirwick et al., [29]. 



a. Quantum partition function 



The energy source approach [25] exploits the observation that the world sheet stress energy tensor acquires a non 
vanishing VEV in the twisted sector vacuum created by the primary twist fields. One is thus motivated to consider 
the five point correlator, ZtIz: Zi] = <^(^)°"-''(^i)°'''(^^)°'-o' (^3)0-0/ (z4)> Q]-,^g^jj^g(j inserting T(z), along with the two 
six point correlators involving the insertion of bilocal quadratic products of the coordinate fields, 

, N / X < -^dzX{z)dy,X{w)a^e{zi)(Je{z2)(J-e'{z2,)ae'{zi) > 

g{z, w; Zi) = g{z, w) = s — — — -— , 

< a-e{zi)a0{z2)cr-0'[Z3)a0^[Z4) > 

h(- \ - < -^d^X{z)d^X{w)a-0{z-i)ae{z2)(T-e'{z3)a0'{z4) 

n[Z,W, Zi) = n{Z,W) — r r r r— . {■^■'^) 

< a-0{zi)a0{z2)a-0>{z3)a0'{z4,) > 

The correlators g{z, w), h{z, w) are mcromorphic functions of z, z and ui, whose singularities are fully determined 
by the operator product expansion of the primary twist fields, as given in Eqs. (11.18). The fact that the short distance 
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expansion in the limit z ^ w of g{z,w) exhibits a double pole term, allows one to relate this function to the correlator 

ZT{z;Zi) as 



—dX{z)dX{w) 



1 



{z — w)' 



+nz)+- 



g{z,w) 



1 



{z — w)' 



+ Zt{z; Zi) + ■ 



(A.2) 



where the dots represent the higher order terms in the expansion in powers of [z — w). Since no singularities are 
present in the limit z w, one concludes that the bilocal correlator h{z, w) has no pole terms and so must tend to a 
constant in this limit. Writing the short distance expansion of ZT{z;Zi) in the limit z ^ zj, for fixed index /, with 
the help of the familiar operator product expansion 



yields the differential equation relating the correlators Z{zi) and Z^iz; Zi) 

^ In Z{zi) = {z- zi)Zt{z; Zi) - 

ozi z — Zi 



(A.3) 



(A.4) 



The holomorphy properties of the functions g{z,w) and h{z,w), as the variables z, z and w approach the insertion 
points Zi, allows one to write the following general representations for these functions, 



g{z,w) = u)e,e'{z)uji-0,i-e'{w) 



P{z,w) 



■ A{zi) 



{z — w)"^ 

uj0,e'{z) = {z- zi)-\z - Z2f-^{z - Zi)-^'{z - Zif'^, 

W-i-0^l-0'{w) = {w- zi)^~'^{w - Z2)~^{W - zzf~'^{w - ZiY 

= {Z- Z^f-\Z - Z2)-\z - Zsf-\z - Zi)-^', 
2 



, h{z,w) = u)i-e,i-e'{z)LUi-0^i-e>{w)B{zi), 



P{z,w)= ' 



aoo + aoiz + aww + a2ow 



+{aiiw + a2iw^)z + (ao2 + Q!i2W + a22w'^)z'^ 



(A.5) 



where the coefficients A{zi), B{zi) (to be determined in the sequel via the global monodromy conditions) and the 
coefficients cijk of the second order polynomial P(z, w) arc hmctions of the insertion points Zi. The polynomial P{z, w) 
is determined by the requirement that one reproduces the pole structure, g{z,w) {z — w)~^, in the limit z ^ w. 
Expanding the various factors of g{z, w) in powers of {w — z) leads to the two sets of relations: 

• P{z, z) = [uje,e'{z)LOi-e,i-e'{z)]~'^ = ^{z - Zi) 

i 

=^ aoo = 21^2^324, "10 + aoi = -{Z2Z3Z4 + Z3Z4Z1 + Z4Z1Z2 + Z1Z2Z3), 

an + ao2 + a20 = -^324 + z^zi + Z1Z2 + Z2Z3, a2i + o.\2 = -{zi + Z2 + Z3 + Z4), 0:22 = 1, 

6»'-i e 6' 



P'{z,z) 



P{z,z) 



1^1-9,1- 



(z) 



Z — Z\ Z — Zz Z — Z2 Z — Z4,_ 
OLIO = 6zzZa(z2 - Zi) + 9' ZiZ2{zA - Z3) - Z2Zi{zz + Zi), 

2a2\ + ai2 = -[^zi - Z2) + e'{z3 - Z4) + 2(z2 + Z4) + zi + Z3], 

an + 2q;2o = -[{0 - l)(z3^^4 + Z2Z4 + Z2Z3) + {9' - l){z2Z4, + 2;i2;4 + -21-22) 

-6'(Z3Z4 + Z1Z4 + Z1Z3) - 6''(Z2Z3 + Z1Z3 + Z1Z2)], 



(A.6) 



allowing to express the solution for the T(2;)-inserted correlator as 



1 



Zt{z\ Zi) = uj0fi'{z)uj'-4_0^^_Q,{z)P' {z, z) + -uJe,e'{z)uJi_e,i-0'{z)P{z, z) 
+uje,e'{z)u:i-0^i-0'{z) (^A{zi) + \p"{z, z)^ 

= -iJe,e'{z)ij[_0^i_g,{z)P^{z,z) + ^(jj0,e'{z)uj"_g^-i^_g,{z)P{z, z) 
+{A{zi) + a2o + a2iz + z^)P{z, z), 



(A.7) 



41 



with the primes standing for the derivative 9^, namely, P'{z,'w) = dwP{z.w), P"{z,'w) = d^P{z,w). The equations 
for the coefficients, ai2, a2i, are independently solved, with the solutions given by, ai2 = {6 — — {6 + l)z2 + 
{9' — l)z3 — {9' + l)z4, a2i = —[9zi + 9'z3 + (1 — 9)z2 + (1 — 9')z4]. The condition on the double pole term gives eight 
linear equations for the nine coefficients, a^. However, since the coefficient a2o can always be absorbed inside the 
function A{zi), this coefficient remains arbitrary as long as A{zi) is unspecified. We make the convenient choice [29], 
Q!2o = + 9'){zizs — Z2Z4) + [9 — 9'){z\Zi — Z2Z3) + 2Z2Z4], which then fully determines P{z, w). Evaluating the 
function g{z,w) by taking first the limit, w — > 00, and next the limit, z — > Z2, leads to the differential equation for 
the correlator of interest 

OZ2 I Z23Z24: Z\2 \ZZ2 2^24/ 2^21-2232^24 

where wc use the notation, = z^ — Zj, and have set zj = Z2. Similar equations apply by letting z approach the other 
insertion points. Thanks to the invariance under the SL{2, R) subgroup of the conformal group, one can arbitrarily fix 
three insertion points. With the choice of insertion points along the real axis boundary, z\ = Xi = 0, Z3 = Xs = 1, 24 = 
Xi = X ^ 00, leaving the single real variable, Z2 = X2 = x [—00,00], one can reconstruct the full dependence by 
noting that the Mobius group invariant functions of four variables in C can only depend on the harmonic ratio variable, 
X = ffjf^- (Note that (1 — a;) = ffff^-) The so far unspecified coeflacients of the polynomial P{z,'w) are now given by 

the simpler formulas, Q!2o = X[^{9-9')+x{l- ^)], an = X{-9 + 9' + l + x9'), = f [-a;+ (6*- 6'')(1 - a;)]. It is 
convenient to absorb the known dependence on the variable X by considering the reduced correlator and coefficients 
identified by hats, Z{zi) = Z{x){x- Xf^'^-'^'\ A{zi) = {x-X)A{x), B{zi) = {x-X)B{x). The differential equation 
for the reduced correlator reads then as, 

dx 2(1 (1-x) X x{l-xy ^ 

To proceed further at this point, one must separate out the zero mode part of the coordinate field, associated with the 
classical zero modes, from its quantum or oscillator part. The additive splitting, X{z) = Xf.i{z) + Xg„(z), produces 
the multiplicative splitting of the correlator into classical and quantum factors, Z{x) = Zqu{x)Zf.i{x), [Zci{x) = 
e^"^"''^)] where the classical factor represents the contribution from the world sheet instantons and anti- 
instantons, defined as the Euclidean space coordinate field solutions of the world sheet equations of motion, 
dzXci{z) = 0, dzXci{z) = 0. The general solutions with the appropriate singularities at the insertion points, Zi, are 
described in terms of the meromorphic functions introduced earlier as, dzXci{z) = bue^$'{z), dzXci{z) = cu!i-0^i-0'{z), 
where the complex constants b, c represent free continuous parameters. 

Wc postpone the discussion of the classical partition fmiction until the next subsection, and continue with the study 
of the quantum factor. The coefficients A{zi), B{zi) are determined by imposing the global boundary conditions 
requiring that the quantum components of coordinate fields, unlike the classical ones, are single valued functions on 
the cut complex plane C . For definiteness, we consider the four point amplitude with the configuration of open string 
sectors, < V-0^(^i:)^A){x\)V0^(^A,B){x2)V-ei ^(^B,c){x3)ye' ,{C,d){xa) >, setting the boundary along the real axis of C+, with 
Zi = Xi. The embedding in Tf is described by the four-polygon displayed in Figure 8, with vertices and edges given by 
the images of insertion points, xi, • • • , x^, and arc segments, Sa = (a^i, X2), {x2, X3), (2:3, 2:4), (,X4, xi), along the disk 
boundary. On the four segments covering the periodic real axis, [A,B,C,D] = [{xi,X2), (0:2, xa), (a;3,X4), {xi,Xi)], 
the complex coordinate fields obey the relations, 

e-i'P'sdX^iz) - e''^sdX\z) = 0, e-'-^sBx^iz) - e''^' dX\z) = 0, [S = A, B, C, D] (A.IO) 

with the interbrane angles assigned in the various intervals as, (j>D = 0, (pA = (l> = ""^j 4'b = 0, (j)c ~ 4i' ^ ttO' ■ Since 
the vector space of contours is of dimension 2, it suffices to consider the integrated form of the boundary conditions 
for two independent contours. We choose the contours Ci, C2 along the branes A and B, with angles 4>a = t^9 
and (pB = 0- To establish contact with the formalism used for closed strings, it is convenient to consider instead 
the associated closed contours (or cycles) surrounding the intervals for the open contours, Ci € C+. The cycles 
Ci G C are defined by adding to the contours Ci their refiected image with respect to the real axis boundary, 
Ci = Ci — C'i, [i = 1, 4]. This definition of the fields dispenses one from the need to define the analytic continuation of 
fields on the Riemann sheets of the cut plane C. We now consider the contour integrals on Ci for the sum of the two 
boundary conditions in Eqs.(A.lO), and reshuffle the two pairs of terms so as to express these as cycle integrals on the 
basis of independent cycles, Ci by using the transformation under the reflection, dX — > e^^'^^dX and BX — > e'^'^^^BX. 
The resulting integrated boundary conditions on the coordinate fields 

= / [(5X - e^^-^^aX) + (ax - e^'-^^aX)] = / {dX + BX)= [ (BX + BX), 

JCi JCi-C JCi 
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(A.11) 

arc then used to express the trivial monodromy eonditions on the quantum component of the six-point correlators, 
= dzg{z, w) + dzh{z, w). This is our first encounter with the mixed type correlator, h{z, w). The two integrals 
along the contours Ci, C2 give the relations needed to determine the two unknown coefficients, A{zi), B{zi). We 
note that /^^ ^ JJ', J^^ /J . Substituting Eqs. (A. 5) for g{z,'w), h{z,w) into the integrals over the cycles Ci, C2, 
gives the pair of equations for i = (1, 2), 

A(x) I dzu!e,e'{z) + B{x) / dzu)i-e,i-e'{z) 

JCi JCi 

= -a;^^— ^ / dzioe,e'{z) - {I - 9') [ dz{z - x)u;e,e'{z) = x{l - x)dx [ dzuje,e'{z), 

(A.12) 

where the last step is deduced by integration by parts. The resulting two linear equations for the coefficients 
A{x)^ B{x), obtained with Ci, C2, can be expressed by the matrix equation, 

J dz.,^Hz) J dzc,,.e,^Mz)\ ( Aix)\ ( kdz.eAz)\ ^^,3) 

Ic,dzue,e'{z) J^JzuJi-e,i-e'{z) J \b{x) J ^ ' \J^^dz^e,e'{z) J ^ ' 

These equations are readily solved for A{x), B{x), with the explicit formula for A{x) reading as 

Mx) - xil - X) ""'-''^^-^^^ - id. k ^o,s')(k ^1-0,1-0') I , _ . . 

[A^) = 7^e"(4-f)«-2-«'[(/' ue,e'){f u^i-e,i-0') - { I toe,e'){ f u>,.e,i-e')] (A.14) 

V / J Ci J C2 J C2 J Ci 

where to simplify the notations we have omitted writing the integration measure factors dz or dz in the cycle integrals. 
For the purpose of simplifying the final results, we have set the so far undetermined constant normalization of I{x) 
to a conveniently chosen value depending on the parameter p. Explicit analytic formulas for the cycle integrals will 
be presented shortly. The total derivative structure of the function, A{x)/[x{l — .x)], in Eq. (A.14) transforms the 
differential equation for Zgu{x) in Eq. (A.9) to the simple form, d.ln{Zgu{x)li (x)) = -0{l-0)/x+[0{l-0')+0'{l- 
^)]/[2(l — x)], which is readily solved to give the final form of the quantum partition function, 

Z'^^ix) = - x)-i''''-'-''^-^''^'-'''h-Hx), (A.15) 

where we have reinstated the complex plane label, /. The constant can be determined by making use of the pole 
factorization of the correlator. Note that we have corrected a confusing misprint in Eqs. (3. 21) and (3.37) of the work 
by Biirwick et al., [29]. These would have given a wrong dependence of the exponent of (1 — x), lacking the symmetry 
under 0^0', which our above result exhibits explicitly. The integrals over the two independent cycles have the 
following analytic forms 

/ dzue,e'{z) = 27rie2"(^-^)G2(a;), 

JCi 

[ dzuje,e'{z) = 2isin(7rp^)(l - x)^'^' e'''^^P-^^^+^"^BiHi{l - x), 

JCo 



L 



dzuii-e,i-0'{z) = 27rie2"(^-^)Gi(a;), 



c 

j dzu)i-e,i-e'iz) = -2isin(7rp^)(l - a;)^'-^e*'^Kf-2)^+^'lB2^2(l - x), 

Bi-Bi9,l-9)= B,-Bi0,l-0)= 

Gi{x)=F{0,l-0';l;x), G2ix) = Fil - 0,0';l;x), 

Hi{x) = F{e,l - 0';1 + 9 - 0';x), H2{x) = F{1 - e,e';l - + 0';x), 

Pia, b; c; x) = ^^JJ""^ dtt>>-\l - ty-'-^l - xt)-'' 



r(6)r(c-6) Jo 



(A.16) 
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where B{a, b) and T(z) designate the Euler Beta and Gamma fimctions and F{a, b; c; z) the Riemann Hypergeometric 
function. Although we have kept track in the above formulas of the complex conjugate functions, signalled by the bar 
symbol, this distinction is not relevant here to the extent that the restriction of the various functions to the real axis 
are real. 



b. Classical partition function 

We now return to the calculation of the classical action factor, Z^i = 'J2ci ; where the summation extends over 
the non-trivial monodromies around the cycles Ci of the coordinate field classical solutions, Xci{z, z). For orientation, 
we recall that in the simpler case of 5^ circle compactification, the classical summation extends over the 1-d lattice 
of solutions with winding numbers, ^/2AcX = \/2 J^i dX = 2ttv, [v G Z\. In the torus case, the monodromies 
of the coordinate X^ around the cycles, Ci, C2, belong to the 2-d grand lattice generated by the cycles along the 
fe^Z-DeB— branes, images of C\, C2, possibly shifted by the non-vanishing distance between distant intersection 
points. Denoting the intersection points along the branes D6a and D6b by f^, and in correspondence 

with the images of the insertion points xi, X2 and ^3, one can express the non- trivial monodromy conditions on 
the classical components of coordinate fields, A^X^ = J^^{dzdX^ + dzBX^), for the cycle basis Ci,2 mapped to the 
f-branes A, B, in terms of the 2-d lattice generated by the Z^-brane pair, L^, Lb, as 

V2Ac,X' = 27rvA = 2n{l - e^'^'){ff - + paLa), 

V2Ac,X' = 27rvB = 27r(l - e^'-'){fi - fi + pbLb), (A.17) 

where we have assumed, for simplicity, the brane lattice to be generated by the integers p/i, pB G Z. The relation 
Ac-X = Ac^-c'.X = (1 — e^^^^)X, [i = 1,2] follows from the 27r0 rotation of coordinates along the reflected 
contours. If we had used the open contours in the upper complex plane, C+, as in the presentation by Cvetic and 
Papadimitriou [36], similar formulas to the above would hold, but without the angle dependent factors, (1 — e^"^). 
To relate the integer quantized winding numbers pA, Pb with the coefficients &, c in the instanton and anti-instanton 
solutions, it is convenient to introduce the 2-d bases of holomorphic and antiholomorphic solutions, X^*^ [i = 1,2] 
defined by 

V2dx'^'-\z) = biiJe,e'{z), V2dX^^\z) = b2t09,e'{z), 

V2dX^^\z) = cioji-e,i-e'{z), V2dX^^\z) = C2U>i-0,i-0'{z). (A.18) 

The condition that the X^*) constitute the dual basis to the basis of cycles, Ci, is expressed by the orthonormalization 
relations, 

2iTdij = \/2Ac,X(J^ " ^/ '^^^^^^^(-^) + ^/ dzBx'^^^z), [i,j = 1,2]. (A.19) 

Expressing the general classical solutions as linear superpositions of the two solutions in Eqs. (A.18), 
V2dX = buj0,0>{z) ^ V2{vAdX^^'^ + usax'^)) = {b^vA + b2VB)uJe,e'{z) 

V2dX = cu!i-0,i-0,{z) = V2{vAdX^^^ + vbBX^^^) = {civa + C2VB)u)i-0,i-e'{z), (A.20) 

allows identifying the coefficients b, c as. b = biVA + &2'-'b and c = ciVa + C2Wb, where va, vb are the lattice vectors 
defined in Eqs.(A.17). The four conditions, 2tt = V^Aci^:,X^^^''^'^\ = V2Aci,2-'^(^'^(^\ may now be used to obtain 
four linear equations for the four constant coefficients hi, Ci, whose solution is given by 

bi = ^ [ dzLOi-0^i-0'{z), 62 = -^ / dzu}i-0^i-0'{z), 
ci = -— / dzuj0fi:{z), C2 = — / dzuj0fi:{z), 

•J JC2 JCl 

[J = ( / oj0,0,{z)){ / uJi-0,-i-0'{z)) - ( / u)i-0,-i-0'{z)){ / oj0,0,{z))]. (A.21) 

JCi JC2 JCl JC2 

The final step involves substituting the above solutions into the string classical action, 

S^i = / d^zidXBX + BXdX) 

2-Ea' Jc 
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27ra' 
47rK:,; 



12 



(fz\uj0fi>{z)Y 



d'^z\uji-e,i-0'iz)\" 



1,2] 



+ C1C2 / (fz\LJi-e,i-e'iz)\'' 
Jc 



(A.22) 



and evaluating the complex plane integrals for the holomorphic and antiholomorphic differentials by making use of 
the analytic continuation method of Kawai et al., [116]. Using the analytic formulas for these integrals in Eq. (A. 16) 
along with the previously derived formulas for the coefficients, bi, Ci, in Eq. (A. 21) gives the following results for the 
functions Vij determining the classical action 



^11 = 



4P{x] 



sin npO 



+\B^HiY{2^{B2G2H2) + C00,\G2\ 



V22 



1 



-2Jft[GiG2(B2Gi^2 + B1G2S1) 



4P{x) 

1 £i^gi.pe[5^^^^^(25i(BiGi5i) 



4/2 (x J TT 

-B,GiH,{2?ft{B2G2H2) + C0e^\G2\% 
[cg g, = 7r(cot(7r6') - cot(7r6''))] 



ces'lGif) 



(A.23) 



where I{x) is given by the analytic formula quoted in Eq. (A. 16). To obtain the classical partition factor, Zd, it 
now remains to sum over the 2-d integer lattice generated by the pairs of branes. The comparison with the classical 
partition function for closed strings in orbifolds in Ref. [29], whose notations we have closely followed, indeed shows 
that the open string correlator is the square root truncation of the closed string correlator. With our choice of closed 
contours, the parameter p in the definition of Eq.(A.16) for the function I{x) must be set to unity, p =1. 



c. Useful limiUng formulas 



To ensure that the series for the instanton sums converge at the end points x = 0, 1, thus avoiding the end point 
singularities in the x-integral of string amplitudes, it is necessary to make use of the Poisson resummation formula. 
This is also needed for a proper identification of the particle exchange pole contributions in the field theory limit. We 
record below the limiting forms of some intermediate results at the end points of the interval, x £ (0, 1). The limiting 
formulas at a; — > are given by 



sm{n0) , ,5(6') 
log . 

TT X 



F{e, l-e;l;x)^ 1, F{e, 1 - 6*; 1; 1 - x) 
Bifl'i(x) ^ In -, B2H2{x) In -, 

4 x Ssm{-Kpe)\Ti{5/x) 



4sin (ttp^) ln((5/.x) ' 



ln^(6') = 2V'(1) - - V'(l - 0), \ri5x(e,e') = 2ij;(l) - V(l - 9') - iIj{0), 

i62{e, 6') = 2^{i) - v(i -6)- ^{6'), d{e, e') = [5.{6, ey 

^(1) = -7, ^(1 - z)- il>{z) = TTCotinz), T{z)r{l -z) = 



sin(7rz) 



(A.24) 



where ^(z) = il)^^\z) = denotes the Digamma function and 7 is Euler constant. We also quote the series 

representation for the Hypergeometric function, 



mi-^:l;.).^^ 

n=0 



{e)n{l - 9)n{l - xYl^^5n{9) 



(n!) 



1^2 



M5n{e)) = 2^{n + 1) - + n)-^p{l-e + n), (0)„ 



+ n) 



+ l)---{e + n-l)], (A.25) 
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in order to exhibit the singular logarithmic behavior of the x-integrand near the end point, x = 1. The limiting 
formulas at x ^ 1 are given by 

sin7r(6»-6''),-, .e'-ef ^{0-0') 



sin7r6l' ^ ' \T{9)V{l-9') 

Note that G\{x) and G2{x) differ by the substitution, 9^9'. While the above formulas are specialized to the case, 
9 > 9', the results in the opposite case, 9 < 9', are deduced by using the symmetry under the substitution 9^9'. 



d. Trilinear Yukawa coupling 



The constraints on the four fermion correlator implied by the conformal symmetry can be used to infer an explicit 
formula for the trilinear Yukawa couplings. We consider the four fermion string amplitude discussed in Subsec- 
tion III A. In the limit x ^ I, the four point correlator factorizes into the product of three point correlators as, 
Z{xi) < V-0{O)Ve{x)V-0,{l)Ve,{X) >~< Ve{x)V-e,{l)V-0+g,{z) >< V0-0,{z)V-0{O)V0,{X) >, dominated by the t- 
channel exchange of the mode with vertex operator, V-0+0'{z). Retaining the leading contribution in the limit a; — > 1 
yields the limiting formula for the four fermion string amplitude 



A' 



7r( 2|cot(7r6'^) 



cot(7r6''-f)| 



r(i -0' + e'') 



T(l-9^)T{9'^) 



(A.27) 



To forbid in the classical action factor the singular contribution from uji-0^i-0i(x) at a; — > 1, one must impose the 
restriction, = dX = cuLi\^0^i^ei{x), which implies, = c = CiVa + C2Vb- Using the formulas for ci,2 in Eq. (A. 21) 
reduces the classical action to the form 



lim Scl{,x) = \VA?[Vl^ + 1-1^22 - 2K(Fi2-)] 
x^l C2 C2 

= - Jl^^^n cot(7r^) - cot(7r^')l = 7r—,\MPALA + dA)\ 
z Zwa 



.sm{TT(9 -9'))sin{iT9) 
sin(7r6'') '' 



(A.28) 



where the coefficient of (27ra')~^ in the last step above equals twice the area of the triangle with base, 2n{pALA + dA), 
opposite angle ^{9 — 9') and adjacent angles tt^, tt^'. To illustrate in a concrete way the t-channel factorization 
with respect to the exchange of a scalar boson mode, we consider the orbifold model discussed in Subsection V A. 
The correlator and its pole factorized part, < qit^^^s d1 >~< qldf H >< H^qid% >, couple the conjugate localized 
modes from the open string sectors, q = {a',b) = {b',a), d'^ = (a, c) = {d ,a'), H = {b',c) = {c',b), where the brane 
embedding of the world sheet is given by the four-polygon with sides a',b,a',c'. The Dirac spinor matrix element 
may be simplified by using the Fierz identity: iSi = {q2Lj'^qiL){d'^]^Jij,d'ij^) = 2{q2Ld^B){qiLdAR)* . The trace factor 
dependence on the CP factors can be put in the factorized form , 71 = Trace(A^ l A^d AgjAdc) = 'P|Tr(Aqt A^ct A^)p, 

where T' is a calculable coefficient that we do not attempt to determine here. The comparison with the H mode 
t-channel exchange contribution to the amplitude, using the definition of the trilinear Yukawa coupling, Leff = 
Trace{XgX(ici Xjj)Y_g g, g_gr {q^diiH) + H. c, yields the final formula for the Yukawa coupling constant 



Y-0,0>,0-0' = ^/2^s^Y[ 
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(A.29) 
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The argument in the exponential identifies with the area of the target space embedding given by the triangle with 
adjacent sides lying along the branes a' and c'. The above result coincides with similar ones quoted previously in 
Refs. [36, 56] except for the constant normalization factor. Besides the factor ^/2V^ the ratio of our Yukawa coupling 
constant to that obtained by Cremades et al., [56] amounts to the factor of (27r)^. Identifying the prefactor of the 
exponential with the dependence of the Yukawa coupling constant coming from the canonical normalization of the 
kinetic energy terms in the action, Y-g^g'^-gi oc Y\j[K{~6^)K{6 ^)K{6^ — 6 ^)]~'^ , allows us to identify the Kahler 

potential for the twisted modes as, K cx { -p^i'LJi-j )^ Cg, Cgi. This result agrees with that obtained in Ref. [43]. 

The string amplitude factorization in the limit a; ^ oo can be compared with the u-channel exchange pole contri- 
bution in order to identify the multiplicative normalization constant C„ in Eq. (A. 15). The limiting behavior of the 
string amplitude, restricted for simplicity to the quantum partition function factor, 

-^'-(^--^^)< sin.(^'^-V) ) ne')v\e'!-e')i-u + ir-e''\y ^^-^^^ 

The comparison with the leading u-channel pole exchange contribution determines the multiplicative normalization 
constant introduced in Eq. (A. 15) as, Ca = n/(2sin7r6'^) 2 . One also deduces the expected result for the string mass 
spectrum of modes with the quantum numbers of vector bosons, a'M^ = — 0'\. 



APPENDIX B: Z3 ORBIFOLD MODELS WITH INTERSECTING D6-BRANES 

We here present an encapsulated review of the construction due to Blumenhagen et al., [66] of the orbifold- 
orientifold models with D6-branes. The orbifold T^/ (Z3 + Z^QTZ) uses factorisable 6-d tori symmetric under the 
generator, @ = diag{Q^,Q^,Q^), and the orientifold symmetry, TZ, acting on the complexified orthogonal basis of 
coordinate and spinor fields X' , ip' , by the complex phase rotations. 8^ = e^^" , [v' = (i, i, — |)] and the reflection 
about the imaginary axis, TZ = TZ^ X' —X'. The orbifold Hodge munbers, h^^^^^ = 0, /i^-'^'^^ = 36, entail 
that all the complex structure moduli are frozen, while the 36 complex Kahler moduli decompose into 9 untwisted 
moduli, Tjj, and 27 twisted (blowing-up) moduli labeled by the fixed points. The two inequivalent solutions for the 
2-d symmetric lattices are described by the pairs of cycles, 

A : ei = 1, 62 = e^"/^ B : d = e-^"/^ 62 = e+"/^ (B.l) 

having the same complex structure moduli, Ua = C^b = ff- = 5 +^"^) but differing by the action of TZx which acts as 

A : ei ^ ei, 62 ^ ei - 62, B : Ci ^ 62, 62 Ci. (B.2) 

In the sequel, we consider the alternative solution for the lattice B, using the cycles, ei = 1, 62 = "^e*'^^^) with the 
complex structure moduli, J/b = = | + ^^'^ reflection symmetry realized as, Ci — > Ci, 62 — > Ci — 62. The 

resulting tori lattice solutions, A and B, are described by the following data 

e4 1 , .VS , . -VSr^ , , ^ , 1, .V3 . f 1 ^ 



B: C7« = ^ = -+z— , T^ = 6 + i— ,[ef = 1, e2^ = -+z-^, ff^= 1 ? ] (B.3) 



^2V3' 2^/3'^ ' ' ' 2 ' 



2 3 



where we have displayed the diagonal complex structure and Kahler moduli as a function of the free parameters, 
r and h ~ 6-1-1, associated with the overall radius and the NSNS field VEV; the complexified basis of cycles, 
e{ = 5R(e{) + i3(e{), 63 = 3?(e|) + i9(e|), whose real and imaginary parts describe the orthogonal components of 
the 2-d vielbein vectors; and the lattice basis metric tensor in the matrix representation gab — £a ' = 3fi(e*ef,). The 
lattice basis decomposition of one-cycles, [Ua] = (riaGi + mae2), [na,ma S Z] defines {na,ma) as the 2-d column 
vectors of wrapping numbers. 

The orbifold and orientifold symmetries act from the left on the column vector of wrapping numbers {ua, ma)^ and 
hence on the right on the column vector of cycles (ci, 62)^, with the matrix representatives for the A and B lattice 
solutions given by 
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These matrix representations obey the operator identity, WZQ^ = The wrapping numbers of the one-cycles 

composing the orbits, [a], [a'] are obtained from those of the representative element, a (rio, ma), by left action with 
the matrices 9^, QfRB^ appropriate to the A and B lattices, 

\i"'a/ \"''ag/ \"<'a/ \"<'a/ \"''a'g/ \i'''aj 

where we use the notational convention, = ilTZQ^ a = (a^)' = {a')N-g- 

Our present conventions for the coordinate system conform to those of Blumenhagen et al., [66]. We use sideways 
tilted Tf tori with 06-plane lying along the real Xi coordinate axis, and the relative Z)6a-brane (96-plane angle 

evaluated by means of the formula, tan^^ = ■ No confusion should hopefully arise from the fact that these 

conventions differ from those in the main text, where we used upwards tilted tori, symmetric under the reflection, 
TZy : X^ with 06- plane also lying along the real axis. The two choices, 'R-x,y, of the reflection are related by 

the modular transformation, — > — * — m^, — > n^. 

The 6-d orbifolds /{Zs + Zs^TZ) arc given by direct products of the above 2-d tori solutions. Taking the 
symmetry under permutations of the three complex planes into account, one finds the four inequivalent 6-d lattices: 
AAA, AAB, ABB, BBB. With the reference O6o-plane wrapped around the three-cycle, ttiss, with wrapping 
numbers, (1,0)"^, the two other planes, 06i, O62, obtained by applying the half-angle rotations, 06g = 0~2 • O60, 
form an equilateral triangle whose sides are interchanged under 8. The three-cycles are represented as direct products 
of three one-cycles, e{ 2 in Tf, [I = 1, 2, 3] as illustrated by the examples, ttiss = tti x tts x tts = ej x ef x ef , 7r235 = 
7r2 X TTs X 7r5 = 63 X X ef, with the intersection numbers, [7r2i-i] • [7^21] = 1, [* = 1,2,3]. The total 06-plane 
RR charge is compensated by including stacks of D6-branes wrapped around the factorisable 3-cyclcs. We restrict 
consideration to setups of I?6-branes intersecting the 06-planes at the origin of the coordinate system. A compact 
representation for the wrapped cycles is obtained by considering the orbifold invariant three-cycles belonging to the 
sub-space of dimension 63 = /i^^'^' -|- 2 = 2, with 63 denoting the Betti number. For the symplectic basis of invariant 
three-cycles defined by 

Pi = V7ri35, P2 = V-K2iz, [V = {! + & + e^)] (B.6) 

having the intersection numbers, pi ■ pj = —eij, = 1,2] where ei2 = —£21 = 1, one can define the three-cycles 
eff'ective wrapping numbers, (i^, Zg), for the equivalence class, [!!„], in terms of the decomposition on invariant cycles, 

na^ZaPl+YaP2= ^ + ^2) ■ (B.7) 

g=0,l,2 / 

The above equations may be used to express the eff'ective wrapping numbers as weighted averages of products of the 

one-cycle wrapping numbers. To adapt ourselves with the conventions of Ref. [66], we rather consider the decompo- 
sition on effective wrapping numbers, Ila = ZaPi + YaP2, including the orientifold images, 

^« =lj:(m.+m'X ^« = -^e (n< -n-O^ 

g ^ I I ^ g ^ I 1 ^ 

\n' = in' + raiui), = [nL + mi,U{% (B.8) 



where the wrapping numbers of the mirror branes are calculated from Eqs. (B.5). The mirror branes satisfy the rela- 
tions, Ya' = —Ya, Za' = Z^- The explicit formulas for (Vo, Zg) in terms of the one-cycles wrapping numbers {n\,m^g) 
are listed below for the four inequivalent 6-d lattices invariant under the Z^ orbifold and orientifold identifications. 



• AAA: Z„ = i 



Ya = nllnlml + (n^ + ml)ml) + ml{niml + nl{nl + ml)) 
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AAB 



ABB 



BBB 



mlinlnl - ml{nl + m^)) + nUnlml + nl{2nl + ml)) 



Za = 

Ya = nl{3nlml + (n^ + 2ml)ml) + ml{ml{3nl + ml) + nl{3nl + 2ml)) 



Ya 

Za 
Ya 



mliSnlnl - mlml) + nl{3nl{2nl + ml) + ml{3nl + ml)) 



3(n2(3n>i + « + 2ml)ml) + ml{nl{nl + ml) + ml{2nl + ml))) 



minimi + nl{inl + ml)) + nl{:inl{2nl + ml) + ml{?,nl + ml)) 



^nl{nlml + [nl + ml)ml) + Zml{3nl{nl + ml) + ml{^nl + 2ml)). 



(B.9) 



These formulas reproduce the resuhs previously obtained by Blumenhagen et al., [66], except for the single mismatch 
in the formula of Za for the lattice AAB, where our above term, —^m\mlml, is quoted in Eq.(A.2) of Ref. [66] as, 



-m\mlml. 



The £)6^g-brane stacks at generic angles relative to the 06-planes, ^ /ctt/S, [g = 0,1,2; fc = 0, 1, 2] realize on 
the 4-d world brane the gauge symmetry, Y\ U{Nf^). The states in the non-diagonal sectors, (a, bg), (a', bg), [a ^ bg] 



nlmi) 



which we denote 



carry the bifundamental representations, IabgiNa,Nb) + Ia'bg{Na, Nb), [lab = Yljiniri 
for convenience as, Iabg{Na,Nb) + la'b^iNa, Nb), [lab = ha = -lab = -[tTo] • [TTb]]- From the multiplicities of the 
bifundamental representations, one determines the chiral spectrum multiplicities (distinguished by the suffix label 
x) by summing the intersection numbers algebraically over the orbifold orbits. The chiral spectrum, obtained by 
summing over the distinct pairs of I?6a/i?6{,-brane images in the equivalence classes, ([a], [b]), without distinguishing 
between the sectors related by the orbifold group, {a,ag) ~ {dhidgh), [9ih € Z^] only depends on the effective 
wrapping numbers 



{Na,Nb): lib 



g 



-{YaZb - YbZa), iNa,Nb) : = -J2la'b, = {YaZb + Ife^a). (B.IO) 



The diagonal sectors, (a, Ug), generate the adjoint representation of the gauge group U (Na) with multiplicity /aog Adja- 

The diagonal sector, {a',ag), includes the subset of T-r "'''^ two-index antisymmetric representations, Aa, with the 

remaining subset ^la'agi^ — rt ^i ) realizing symmetric and antisymmetric representations, Aa + Sa- The resulting 

chiral spectrum of antisymmetric and symmetric representations and the complete spectrum of adjoint representations 
can be expressed by the following compact formulas [66] 



^a'a 



g 1. 1. I O, g 

^'^^a'agi'^i^ — T-r / ) = ^a(^a " -7:), 

g ^ 11/ "■ag ^ 

Adja : 4^* = - J E = n l^al'' 



A. 



7=1 



[l^ak = (nf +nX + l^alB = {ni^ + nimi + ^mf )i] 



(B.ll) 



where pb denotes the number of B tori and ji^l^ = (n^ m^) • ■ (ua^^ '^^)^ denotes the squared length squared 

length of the Tf one-cycles wrapped by the _D6,j-branc. For Na = 1, one must set I^/^ = 0, reflecting the absence 
of the antisymmetric representation in this case. Note that the multiplicities are related to the cycle intersections as. 



-[na]-[nb], Ilb = -K']-[nb], I^^l 



X(A,S) 



■i([7ra.]-[7r,]±2[^,,].[7ro6]). 



The RR tadpole cancellation conditions, expressing the net vanishing RR charges carried by the orbifold orbits 
of £'6ju-branes and the 06g-planes, reduce to the unique relation: N^Z^ = 2. The chiral and mixed gauge and 
gravitational anomalies contributed by the massless fermions present m the above representations are given by 

AaiGl) ^ J2 ^"(-^ab + ia'b) + i^a ~ 4)/^'^ + (Na + A)!^'^ 
b^a 

= ^ 2NbYaZb + {Na - 4)Ya + 2NaYa{Za - ^), 
67^:0 
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^^fe(f/(l)„ X Gl) = Na{-I^, + I^,) = 2NaYaZbdab, AabiU{l)a X U (1)1) = 2NaNi,YaZf„ 
Agrav{U{l)a X R^) = Y^[NaNb{-i^, + I^J] + N,{N, - 1)/^'^ + N,{N, + 1)1^'^ = 3N,Y,. 

(B.12) 

The chiral gauge anomalies are found to vanish automatically, Aa{G^) = 0, once the RR tadpole cancellation 
condition is satisfied. The vector space of U{1)^ factors decomposes into two orthogonal sub-spaces of anomaly free 
and anomalous abelian symmetries, corresponding to linear combinations of U{l)fi with vanishing and non vanishing 
mixed gauge and gravitational anomalies from the massless fcrmions. The latter anomalies arc cancelled, as in 
the Green-Schwarz type mechanism, by including the 4-d anomalous counterterm action descending from the 10-d 
Chern-Simons topological action. The dimensional reduction on T^/Z^ with £>6^-branes wrapped around the three- 
cycles n^, produces 4-d couplings between the gauge and gravitational fields and the dual pair of RR scalar fields, 
C2 = /[^^j C5, Co = /j^^] C3, [dCo = -kdC2] with coupling constants determined by the decomposition of 11^ on the 
basis of dual cycles, pi, p2- The resulting anomalous counterterm action 



(B.13) 



with and R denoting the gauge field strength and curvature tensor two-forms, is seen to induce the tree level 
Co — C2 exchange contributions necessary to cancel the abelian anomalies, ^Jff,, Aab and Agrav The unique 
anomalous U{\)x factor here is described by the charge and gauge field strength, Qx — J^fi -^t^^t^Qt^ -^po- = 
A^^F^F^^. Finally, the uncancelled tadpole for the dilaton field, <&, arises through the D6-branes scalar potential, 

vm = e-^{j:^N,Ui\Li\-2). 



APPENDIX C: BARYON NUMBER VIOLATION IN GAUGE UNIFIED THEORIES 

The baryon number violating processes arise in gauge unified field theories [117] from the tree level exchange of 
gauge and scalar bosons gaining mass through the unified gauge symmetry breaking. The massive gauge bosons 
carrying leptoquark quantum numbers with respect to the Standard Model group include the modes, [{y) + ix)] 
[(3, 2)| + (3, 2)_|] in 5(7(5) and 50(10), and [{^',) + (I',)] ^ [(3, 2)i + (3, 2)_i], [X, + X] ~ [(3, 1)| + (3, 1)_|] in 
50(10). The tree level exchange of these modes contribute the effective Lagrangian 

^^^^ = 2^^^^^^ + + ^'-^^ + ^T^^^^"'- + + ^-'-^ + 

0.C = -e^^^{u%ja^.dimWiLl''di^L), O.^ = -e^^^u'^^^i^^^d^^LmRl'^diaR)] (C.l) 

which initiates the quark and lepton subprocesses, u + u ^ d'^ -\- , m -|- u — > u'^ + e+. The change from gauge to 

mass bases uses the transformations of quark fields, {qL)gauge = Vl\qL)mass, {qDgauge = Vlf{ql),nass, {qR)gauge = 
{qRjiyiass, such that thc diagonalization of the quark mass matrices, Xj^qiq^, {X'^)diag = VXV^ , [q = u, d] introduces 
the CKM flavor mixing matrix, V = V^V^^ . In the minimal SU{5) unification case, assuming the mass generation 
from only 5, 5 Higgs boson multiplets, the flavor mixing transformations, d'^ = Vijdj, — > e~^'^'u1, — > e^' = Vijej, 
yield the simplified effective Lagrangian for the first generation of quarks and leptons, 

Leff = ^^[(1 + Kd^Oei + Oe^ + KdO.j^], (C.2) 

where Oe|, Oe^ arise from the couplings, 10^ • loj • lOj • lOj and 10, • lOj • 5j ■ sj, which initiate the proton decay 
processes, p ^ 7r° + and p ^ n° + e^. The current experimental bound, t{p tt^ + e+) > 5. 10^^ yrs is to be 
compared with the predicted proton partial decay width 

T{p ^n' + e+) = 8wC^4Al\(ipf^[l + (1 + cos^ 0cf] 

1 , Mx ,_4, f3p 



5. 1036 yrs ^ 3. 1016 GeV ' ^ -0.015 GeV^ ' 
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K3 = l + —^{D + F)c^2, Arc^ 3.6, U - 139 MeV] (C.3) 

where we use self-explanatory standard notations. The comparison with the experimental lower limit on Mx from 
the nucleon decay rates has long been known to favor the conventional supersymmetric version of grand unification 
over the non-supersymmetric one. However, new possibilities have been opened by the current studies [118], using the 
alternative schemes for the Yukawa coupling constant unification and quark and lepton flavor mixing effects, especially 
those involving orbifold compactification of the extra space dimensions to realize grand unification models [119-121], 
or TeV scale models [122] with variant selection rules on the B, L number non-conservation. 

In supersymmetric grand unification, additional baryon number violating contributions arise from tree level exchange 
of color triplet matter fcrmions between pairs of matter fcrmions and sfermions. The largest contributions arise from 
the B — L conserving F term operators, QQQL, U'^D'^U'^E'^ . Accounting for the quark flavor mixing, leads to several 
quartic couplings between fermions and sfermions of which a few illustrative examples are 



^^^^ = -w^ 



+ H. c. 



1 



H. c. (C.4) 



For minimal SU{5), the trilinear superpotential coupling the matter supermultiplets to the color triplet Higgs 
supermultiplets, He, He, is described in terms of the Yukawa coupling constants of quarks and leptons, A"''^'*', as, 
W = A^ (iQiQj + u1e'^)Hc + {X^fQiLj + X'fj^u'[d'^)Hc, where the choice of one or the other of the down quarks and 

charged leptons Yukawa coupling constants, X^j^, reflects on the ambiguity in describing the uniflcation of Yukawa 
interactions at the scale Mx, and hence of the quark and lepton flavor mixing. For instance, the choice of the 
couplings, {X^jQ^Lj + XfjU^dj)Hc, yields the effective coupling constants for the dimension D = 5 operators, 

Mx ^2Mh, Mx Mh, ^ ^ ' 

Analogous results obtain in the flipped SU{5) x U{1) fi case. The four fermion amplitudes are obtained through the 

one-loop electrowcak gaugino dressing of the fcrmion-sfcrmion couplings using the Yukawa gauge couplings, Leff = 
g2[u*i^{'w^ d) + d^{w~u) -!-•••]+ H. c, and the analogous electrowcak higgsino dressing. Since these calculations are 
standard ones [123], they will not be reviewed here. 

The current experimental bound, r(p + vi) > 1.6 x lO'^^ yrs, should be compared with the partial decay 

width for the dominant decay mode 

r(p ^K+ + Pi) = CKA2m\^^^^^^r\Kg\, smecKsific ei) + f{c, d')] 
+K^ll2l cos6'cK2[/(c, ei) + f{u, d') 



1 Mh, 2. ^2 



(1/3 - 3) X 1034 yrs '2. lO^^ GeV ' -0.015 GeV^ 

W U J W 

{v? - d?){v? - 



[f{u, d) = ,,,,,, log -^ + {u^ d), A'^ 10.5, 



- = ^ + + - = ^32./?' ] (^•«) 

where we use self-explanatory standard notations. The large hierarchy in the mass ratio, Mn^/mz, needed to suppress 
the 2? = 5 operators can be explained by several mechanisms [124-126]. The currently favoured proposals use the 
realization of Higgs bosons as Goldstonc bosons [127], discrete global symmetries not commuting with the unified 
gauge group [128-130] or the flavor physics [131]. The contributions to nucleon decay from the V = b dangerous 
operators are strongly to the threshold corrections to the gauge coupling constants [124]. Thus, in the minimal SU{5) 
unification with the choice of couplings, {XfjQiLj + Xfju':dj)Hc, the experimental limit from the nucleon lifetime 
yields the upper bound, Mh^ > 2. 10^^ GeV, which clashes with the lower bound deduced from the high energy 
extrapolation of the Standard Model gauge coupling constants [125]. By contrast, the alternative choice for the 
couplings, {XfjQiLj + X'ijU'id'^)Hc, suitably adjusted to the fermion masses, suppresses the predictions for the nucleon 
decay rate to values below the experimental upper bound [121]. 
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The interest in the higher order baryon number violation, initiated by dangerous operators of dimension V > 7, 
is motivated by the sensitivity of the nucleon decay processes to mass scales significantly lower than those of grand 
unification and by the variant selection rules. The V — 7 operators involve four matter fermions and gauge or Higgs 
doublet bosons, obeying the selection rules, AJ5 = —AL = —1. Two illustrative examples are given by the operators, 
{Q.RlL){lLdR)(j)^ and {d1dR){e'j^qL)(f>\ with </> denoting the Higgs boson electroweak doublet. 

The baryon number violating six quark operators of dimension D = 9 are especially interesting in view of their 
impact on the N — N oscillation and two nucleon desintegration processes. The local operators in the quark fields of 
fixed chirality and color quantum numbers, x and a, enter in three types involving two pairs of spin-flavor and color 
structures 

('-^l)xi>X2,X3 — (^QiXl^^'""2Xl: )('^Q3X2^^'^"<lX2)('^Q5X3^^'^"6X3)^QlQ2a3a4a50;6' 

(0[2,3])xi,X2,X3 — ('"aixi '^0:2X1 )('"c«3X2^^ '^"4X2 )('^a6X3^^'^«6X3)^afa2a3C«4a5a6' 
^aia2a3a4a5C«6 = l^aiasa^Ca^a^ae + ("5 ^ "e)] + [ai ^ 02], 



-^01020:3040506 ^010205 ^030406 + 



(C.7) 



where Xl^ X2, X3 € [L, R], C denotes the 4-d charge conjugation matrix and the parentheses enclose the pair contrac- 
tions of the anticommuting Weyl spinor fields with respect to the SL{2, C) spin group indices. The corresponding repre- 
sentation using the Dirac spinor fields is given by, (Ol)xi,x2,X3 = -(W-xi%l)(<^-X2'^X2)(<^-X3^X3)^'^. (0[2,3])xi,X2,X3 = 
— {u'L^-^d^-^){u'L^^d^.^){d'L^^d-^^)T^^'''^\ where — x denotes the opposite chirality to x, and we have omitted the color 
indices for notational convenience. In the electroweak symmetry limit, the allowed contributions [132] arise for the 
list of six operators: {Oi)rrr, {02)rrr, {03)rrr, {Pz)lrr, (03)llr, {Oi)llr- {02)llr- 

The current experimental observabihty limit on the N — N oscillation parameter. Fab = G^jy = Tosc — 
10~^'^ sec~^ = 10~^^'^ GeV, is to be compared with the approximate theoretical prediction. Tab = Sm ~ 
fl'£FF|'0Af(O)|^, where the model dependent coefficient Qeff is defined through the eflPective Lagrangian as, Leff = 
QeffO-d^^, and the quark wave function of the nucleon is given by the estimate, |V'jv(0)|"^ =< N\OT>='d\N >~ 
10-5 GeV^. 

The simplest case where a tree level contribution to the A(i3 — L) — —2 processes takes place in the left-right 
symmetric gauge theory [107] with gauge group, G3221 = SU{3)c x SU{2)l x SU{2)r x U{1)b-l, by including scalar 
multiplets with color triplet diquark and dilepton quantum numbers, ~ (3,8,1)2/3, A|j (3,1,3)2/3, A^^ 
(1, 3, 1)_2, A^ ^ (1,1, 3)_2, coupled through the renormalizable interaction superpotential, W = q'^q'^A'j^ + I'^l'^A^j^ + 
(A|j)'^A^. The prediction for the coeSicient, gEFF ~< A^ > /{m\q Y, along with the natural assumption that 
the spontaneous breaking mass scale of the U{1)b-l symmetry group and the color triplet scalar mass are of same 
magnitude, vbl =< A^ >~ '"^a^j yields the estimate for the oscillation rate, Fab ~ '*^b\,\'^n{^)\'^, [vbl =< 
A^ >~ rriA' ]• However, the resulting experimental bound, tia^ ^ "vrl > 10^ GeV, clashes with the range of values, 
Vbl ~ 10^^ — 10^^) GeV, generating an acceptable see-saw mechanism for neutrino masses. 

More favorable conditions are offered in the supersymmetric version [134, 135] of the Pati-Salam gauge theory [133] 
with gauge group, G422 = S'f/(4) x SU{2)l x SU{2)r, by exploiting the possibility to compensate a large mass 
scale Vbl by a suppressed loop induced contribution to the V — 10 operators. This is illustrated by the model 
in Ref. [135] which includes the matter and Higgs boson chiral supermultiplets, / = (4,2,1), — (4,1,2), (pi ~ 
(1,2,2), 015 - (15,2,2_), A'^ -_(10, 1,3), A= - (10,1,3), interacting through the superpotential, W = fp(j)i + 
f.r<Pi5 + .rr^" + A^A'^ + {A^A")^ + ^i^i + 015^15 + (/'iA'^A'=(/)i5. The singlet modes are omitted for simplicity. 
The sec-saw mechanism for the neutrino mass requires the VEV, vrl = 0(10^^) GeV, whereas the loop induced 
contribution to the coefficient, gEFF ~ ^si^^ '^j obtained by assigning the mass value, Mac ^ 1 TeV, leads to an 
N — N oscillation rate. Fab ~ v]^l''^w\'^n{Q)\^, which is compatible with the experimental observability limit. 
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TABLE I: Matter content for the Z3 orbifold model with gauge group SU{3) x SU{2) x U{1)y x U{1)b-l- 



Mode 




Ka Kb 


Sector 
Irrep 


{a',b) {a,a) {o,,c) ib,c) {b',b) (c, c') 
3(3,2)1,1,0 3(3,1)2,0,0 3(3, l)_i,o,i 3(1,2)0,-1,1 3(1,1)0,2,0 3(1,1)0,0,-2 


(a, a) (h,b) 
/„^*(8,1) /,^*(1,3) 



TABLE IL Solutions for the SU{5) unified Zz orbifold model vii the lattices AAA and BBB. The column entries give in 
succession the solution number, the wrapping numbers for the three-cycles a, c in the order, [(n/j,m^), (n^,m^), (njl,mj^)], and 
the branes intersection numbers la'ag lacg Icc'^ la'cg in the four relevant open string sectors, [a ,ag), {a,Cg), (c, Cg), {a',Cg), 
with the last three columns corresponding to the elements of the orbifold group orbit, g = 0,1, 2. 



Solution 






^a'oo ^"=0 4cf| ^a'co 


4'ai 4ci 4-,, 


4'a2 4c2 44 4'C2 


AAA 












I 


-32 - 10 - 11 


-1 - 2 - 1001 


0000 


0337 


300 - 7 


II 


-320 - 101 




8000 


-5330 


0000 


III 


0-10-13-2 




8-30-3 


-5630 


0003 


IV 


01013 - 2 




8-30-3 


-5630 


0003 


V 


011021 




0000 


3-33-2 


0602 


VI 


-32 - 10 - 11 


-1 - 20 - 110 


0800 


0030 


3 - 500 


VII 


-320 - 101 




8000 


-5330 


0000 


VIII 


-321 - 110 




0000 


0337 


300 - 7 


IX 


-2 - 10 - 110 




0000 


3 - 530 


0800 


X 


-131010 




0000 


8030 


-5300 


XI 


0-10-13-2 


-1 - 20 - 110 


8006 


-5-33-6 


0600 


XII 


01013 - 2 




8006 


-5-33-6 


0600 


XIII 


-10 - 10 - 13 


-10102 - 3 


0000 


8 - 5 - 57 


-588 - 7 


XIV 


-100 - 121 




0000 


33 - 50 


0080 


XV 


0-10-13-2 




8-50-7 


-58 - 50 


0087 


XVI 


01013 - 2 




8-50-7 


-58 - 50 


0087 


BBB 












I 


b-1 - 11 - 21 - 2 


-11012 - 1 


0633 


-5-30-3 


8000 


II 


-12012 - 5 




003 - 7 


0000 


3307 


III 


0-11-22-5 




0830 


- 500 


3000 


IV 


Oil - 111 




3330 


0007 


000 - 7 


V 


-25011 - 2 


0-11-12-1 


0630 


0-30-2 


3002 


VI 


-1 - 1 - 1101 




303 - 2 


0602 


- 300 


VII 


-1-1-12-12 




0-33-6 


-5000 


8606 


VIII 


-1-11-21-2 




0-33-6 


-5000 


8606 


IX 


-110 - 111 




3330 


0007 


000 - 7 


X 


0-11-22-5 




003 - 7 


0000 


3307 
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TABLE III: Solutions in the broken SU (5) unified symmetry case with three non-parallel stacks, a, b, c for the Z3 orb- 
ifold model with lattices AAA and BBB. The three column entries give the three-cycles wrapping numbers in the order, 
[(n^, m^), (n^, m^), Jti^)]. The Cases I and II refer to the fully and partially constrained searches. 





(n^, ml) 


("c. mi) 


(ni, m'j 


(n^, ml) 


("c. ml) 


AAA {Case II) 
-32 - 10 - 11 
-131010 

01013 - 2 


-2 - 10 - 110 

-321 - 110 
0-10-13-2 


-1 - 20 - 110 
-1 - 20 - 110 
-10102 - 3 


011021 
-10 - 10 - 13 


-32 - 10 - 11 
0-10-13-2 


-1 - 2 - 1001 
-10102 - 3 


BBB (Case II) 
-1-11-21-2 
-12012 - 5 
-25011 - 2 
-1-1-12-12 
0-11-22-5 


0-11-22-5 

Oil - 111 
0-11-22-5 
0-11-22-5 
-1 - 1 - 1101 


-11012 - 1 
-11012 - 1 
0-11-12-1 
0-11-12-1 
0-11-12-1 


-12012 - 5 
Oil - 111 
-1 - 1 - 1101 

- 11 - 22 - 5 
0-11-22-5 


0-11-22-5 
0-11-22-5 
-25011 - 2 
-25011 - 2 
-110 - 111 


-11012 - 1 
-11012 - 1 
0-11-12-1 

0-11-12-1 
0-11-12-1 


BBB (Case I) 
-1 - 1 - 1101 
-1-11-21-2 


-25011 - 2 
-25011 - 2 


0-11-12-1 
0-11-12-1 


-1-1-12-12 
0-11-22-5 


-25011 - 2 
-25011 - 2 


0-11-12-1 
0-11-12-1 



TABLE IV: Solutions for the Z3 orbifold model realizing the flipped SU{5) unified gauge symmetry model with three stacks, 
a, 6, c. 



(nl, ml) 


("b. ml) 


(nl, ml) 


("a. ml) 


("6. ml) 


(nl, ml) 


ABB (CASE I) 
-211 - 21 - 1 

-111 - 210 
0-1-101-1 


0-1-120-1 
-110101 
01 - 1201 


10 - 1201 
0-11-21-1 
10 - 1201 


-211 - 21 - 1 
0-1-101-1 
10 - 131 - 1 


01 - 1201 
0-1-120-1 
-110101 


10 - 1201 
10 - 1201 
0-11-21-1 


BBB (CASE II) 
-1 - 1 - 1101 


- 11 - 110 


-11 - 121 - 2 


0-11-22-5 


-110110 


-120101 
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FIG. 1: Plot of the x-integral 1(6) defined in Eq. (IV. 15) for the four fermion string amplitude in the large radius limit as a 
function of the brane intersection angles in the equal angle case 6 = 6'. The left hand panel (a) shows, for the case with three 
equal angles, a plot of 1(6) versus the common angle, 6 = 6^ =9^ =9^. The right hand panel (fe) shows, for the case with 
two equal angles, plots of 1(9,9,9^) versus 6^, at three values of the two equal angles, 9^ = 6^ = (0.1,0.3,0.5). Because of the 
symmetry of the integral under 9^ 1 — 9^ , same results apply in the cases involving the angles, 6^ = 9^ = (0.9, 0.7, 0.5). 
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m.a nLb nLc 



FIG. 2: Representative disk couplings for three and four point string amplitudes in the SU{5) x U{1) gauge unified model of 

the Z3 orhifold discussed in Subsection V A, The arc segments on the disk boundary are labelled by the branes, a, 6, ■ • • and 
their mirror images, a' , b' , ■ ■ ■, and the crosses by the open string modes. In the first array (I), the coupling I.a refers to the 
Dirac mass generation for localized fermion and scalar modes in bifundamental and adjoint representations, the couplings I.b 
and I.c to the trilinear Yukawa coupling of down quarks and leptons with the electroweak Higgs boson, and the coupling I.d 
to the quartic coupling, u"^^ u'^v'^^i/'^ . In the second array (11), the couplings Il.a, Il.b represent the baryon number violating 
subprocesses in the SU (5) symmetry limit of the Z3 orbifold model for the configurations 10 • 10^ • 10 • 10^ and 10 • 5 • 10^ • 5^ . In 
the third array (IIIJ, the couplings III. a, 1 1 I.b, 1 1 I.c refer to the contributions to the baryon number violating local operators 
in the quark and lepton fields, Oe<=^ , Oe"^ , Ov'^ , listed in Eq. (C.l) of Appendix C. The particle modes are assigned to brane 
intersection points in the convention where one navigates the disk boundary clockwise. 
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A(x) (Z_3, r = 1 S 2, e_B= 
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FIG. 3: Plot as a function of x of the string amplitude x-integrand 21{9;x) in Eq. (V.8). The regularized integrand is obtained 
by subtracting the leading end point singular terms, x) — I(d; x = Q)— 1{9; x = 1)], associated with the exchange of massless 
gauge boson poles. The present case refers to solution V of the Zs orbifold model in Table II for the equal angle amplitude, 
10 • 10^ • 10 • 10^, at coincident intersection points, tA = = (0,0,0), for the two values of the wrapped three-cycle radius, 
rusr = 1 (upper curve) and mar = 2 (lower curve). 




FIG. 4: The ratio of string to field theory amplitudes, TZs/f, in the subtraction regularization prescription for the Zz orbifold 
model is plotted as a function of the inverse compactification mass parameter nisr for three values of the inverse unification 
mass parameter, s = nis/Mx ~ 1,2,4. The results refer to solution V of lattice AAA in Table II, described by the wrapping 
numbers, ni = nl = (0,1,2), ni = (—1,-1,0), mi = ml = (1,0,1), mi = (—2,0,1), and the intersection angles Oi,^ = 
(0.333,0.667,0.879). The nucleon decay amplitude for the equal angle configuration 6 = 6' with matter mode representations, 
10 • 10^ • 10 • 10^, is evaluated in the broken unified gauge symmetry model. The left hand panel (o) shows the case with 
coincident intersection points, €a = es = (0,0,0) and the right hand panel (b) the case with distant intersection points, 
€A = (0, 0, 0), es = (1/10, 1/10, 1/10). 




FIG. 5: The ratio of string to field theory amplitudes, TZs/f, in the displaced regularization prescription for the Zz orbifold 
model is plotted as a function of the inverse compactification mass parameter, rusr = ms/Mc, for three values of the gauge 
unification mass, s = ms/Mx = 1,2,4. The results for the string amplitude refer to solution V of the lattice AAA, with 
the wrapping numbers ni — nl — (0,1,2), nl = ( — 1,-1,0), mi — ml = (1,0,1), ml = ( — 2,0,1) and intersection angles 
^a'a ~ (0.333,0.667,0.879). The nucleon decay string amplitude in the equal angle configuration 6 = 6' in the matter mode 
representations, 10 • 10^ ■ 10 • 10^, is evaluated in the unified gauge symmetry broken phase with the displaced brane separation 
adjusted to the infrared cutoff mass, Mx = ms/s, by means of the formula, M% = sin'^(7rS^)|d^p. The left hand panel 
(a) shows the case with coincident intersection points, ea = es = (0,0,0), and the right hand panel (b) the case with distant 
intersection points eA = (0, 0, 0),es = (1/10, 1/10, 1/10). 
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FIG. 6: Plots of the ratios of string to field theory amplitudes, Ti-s/f, in the displaced regularization prescription for the Zz 
orbifold model as a function of the 'inverse compactification mass parameter msV = rUs/Mc for solution V of the AAA lattice 
with wrapping numbers Ua — Ub — (0,1,2), nia = rrif, = (1,0,1), ric = (—1,-1,0), rric = (—2,0,1) and distinct intersection 
angles for the 10, 5 modes, d'^,^ = (0.333, 0.667, 0.879), = (0.227, 0.333, 0.560). The nucleon decay string amplitudes for the 
matter representations, 10 ■ 10^ ■ 5 ■ 5^, refer to the unequal angle case, 6 6' . The left hand panel (a) shows results using the 
subtraction regularization prescription with the values of the unification mass scale, s = rris/Mx = 1,2,4, at coincident points 
SA ~ CB ~ {0, 0, 0). The right hand panel (6) shows results using the displacement regularization prescription at s — rUs/Mx ~ 1 
with the intersection point distance parameters set at, eA = sb = (0, 0, 0) and tA = (0, 0, 0), es = (1/10, 1/10, 1/10), respectively. 




FIG. 7: Plots as a function of the radius, msV, of the ratio of string to field theory amplitudes, 7ls/f, for the Z2 x Z2 orbifold 

supersymmetric SU{5) gauge unified model using the subtraction regularization prescription. We show results in panel (a) 
for the case involving equal brane intersection angles, 6^ = 6 ^ , and in panel (6) for that involving unequal angles, 9^ ^ 9 ^ . 
The solution is described by wrapping numbers, Ua = (1,1,1), Ub = (0,1,0), rua = (1,-1,^), mj, = (1,0,-1), and angle 
parameters, iva^ = arctan(m^(72 /^i); o,t the values, 02 = i, cti — i, 03 = {02 — ai) — -pr, yielding the intersection angles. 
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FIG. 8: Embedding of the world sheet disk complex plane, z (left), in the 2-torus complex plane, e Tf (right). The 
four-polygon edges. A, B, C, D are identified with the D-branes and vertices Xi, • • • , X4 with their intersection points Xi. 



